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Introduction 



Entanglement measures the nonclassical correlations between the components of 
a quantum system. The characterization of this intimately quantum property of 
physical systems is one of most challenging issues of quantum information and com- 
putation theory. The reason why the profound comprehension of entanglement is 
considered so important is twofold: on the one hand it is spurred by an obvious in- 
terest for the foundations of quantum mechanics, on the other hand entanglement 
represents one of the most important resources in quantum information processing 
and quantum enabled technologies (Nielsen and Chuang [2000]; Benenti, Casati 
and Strini [2004]). Generation of highly entangled quantum states is therefore one 
of the key elements for the realization of the ideas of quantum information. In 
this scenario a dominant role is played by the investigation of one dimensional spin 
chains, as possible candidates for modeling quantum computers (Lieb, Schultz and 
Mattis [1961]; Pfeuty [1970]; Takahashi [1999]). However, experimental and theo- 
retical difficulties, such as decoherence and imperfections in the quantum hardware, 
impose strong bounds on the realization of large scale systems and at the same 
time have boosted a high interest in finite size systems (Shastry and Sutherland 
[1990]; Schwinger [1998]; Osterloh, Amico and Eckrn [2000]; De Pasquale, Costan- 
tini, Facchi, Florio, Pascazio and Yuasa [2008]; De Pasquale and Facchi [2009]). 
The bipartite entanglement between systems of small dimension (such as a pair of 
qubits) can be given a quantitative characterization in terms of several physically 
equivalent measures, such as entropy and concurrence (Wootters [2001]), as will 
be briefly recalled in Chap. 1. A direct investigation of quantum correlations and 
other basic properties of quantum systems, can be also reached through very simple 
tomographic techniques (Paris and Rehacek [2004]). Among them, scattering has 
always been considered a very powerful way to investigate many physical systems 
in a wide range of fields of physics, from elementary particles to condensed matter 
physics. Some very simple examples of state reconstruction through scattering of a 
probe qubit can be found in (Kawabata [2001]; Hida, Nakazato, Yuasa, and Omar 
[2009]; De Pasquale, Yuasa and H. Nakazato [2009]). The relationship between 
entanglement and tomography is "bidirectional" : if on the one hand tomographic 
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techniques can be very useful in order to unveil quantum correlations between 
subsystems, on the other hand it has been proved that quantum effects enable an 
increase in precision in estimating the parameters of a given target state. The lat- 
ter phenomenon is related to the fact that entangled states can evolve faster than 
untangled configurations, employing the same amount of resources (Giovannetti, 
Lloyd and Maccone [2006]). More involved is the characterization of the global 
features of bipartite entanglement when we deal with systems of larger dimension. 
This represents a very appealing issue from more then one point of view. Indeed 
besides its application in quantum information theory and related fields of inves- 
tigation, such as complexity (Mezard, Parisi and Virasoro [1987]), it represents an 
interesting problem in statistical mechanics. This will be the central theme of this 
thesis. We will tackle this problem by studying a random matrix model that de- 
scribes the statistical properties of the purity of one of the two subsystems, namely 
the local purity of the global state. For the sake of clarity, we will divide this work 
in two parts. 

In the first part, we will focus on the case of pure global states, yielding the 
description of isolated quantum systems. We will study the eigenvalues distribu- 
tion for the reduced density matrix of a subsystem A. Before the advent of the 
field of quantum information, this random matrix model was studied by Lubkin in 
1978 (Lubkin [1978]), who was moved by a more fundamental motivation: "With 
Lucretius (Lucretius [98 a.C. - 55 a.C.]) I find disorder of the universe repugnant. 
. . . my favorite key to understanding quantum mechanics is that subsystems cannot 
be isolated by tracing from an enveloping pure state without generating impurity: 
The probabilities associated with measurements develop because the observer must 
implicitly trace himself away from the observed system" . In this paper, Lubkin in- 
vestigated the following problem: if an L dimensional pure bipartite state is chosen 
at random, what probability distribution will describe the eigenvalues of one of the 
two subsystems? This problem results completely defined if we interpret by ran- 
dom pure states those vectors of the Hilbert space selected by the unique uniform 
Haar measure on the unitary group U{L), that is typical states (with respect to this 
uniform measure). From the computation of the first two cumulants for L ^ 1, 
he concluded that for a typical pure quantum state of a large system, the smaller 
subsystems are very nearly maximally mixed, that is "impure" showing almost no 
signs of the fact that the (initial) global state (the universe) is pure. It was only 
almost ten years later that this topic became of central importance in quantum 
information theory, when Loyd and Pagel, apparently unaware of Lubkin 's work, 
after determining the explicit expression of the distribution of the Schmidt coeffi- 
cients of the reduced density matrix, computed the approximate expression for the 
typical entropy (Lloyd and Pagels [1988]), whose exact formula was conjectured 
by Page (Page [1993]) and later proved in (Foong and Kanno [1994]; Sanchez- Ruiz 
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[1995]; Sen [1996]). For the measure on the Schmidt simplex, see also (Zyczkowski 
and Sommers [2001]). All these results confirmed the conclusions of Lubkin, which 
translated in quantum information language sound: if the global system is divided 
in two subsystems A and B, the amount of quantum information contained in the 
whole system is greater then the sum of the information in the separate parts, that 
is typical (random) states are entangled with high probability. 

Random states play a key role in quantum communication algorithms, such 
as quantum data hiding protocols (Di Vincenzo, Leung and Terhal [2002]) and 
superdense coding (Harrow, Hayden and Leung [2004]). This has spoored the ap- 
pearance of many algorithms for generating random states (Zanardi, Zalka and 
Faoro [2000]), based also on chaotic maps (Bandyopadhyay and Lakshminarayan 
[2002]; Scott [2004]), pseudo-integrable maps (Giraud and Georgeot [2005]) or op- 
erators (Weinstein and Hellberg [2005]) and sequence of two-qubit gates (Znidaric 
[2007]). Ideally, we would like to compute the complete expression for the proba- 
bility distribution of the eigenvalues of the reduced density matrix. Alternatively, 
one should compute all its moments or cumulants. In particular, the exact ex- 
pression for the first three cumulants have been computed in (Scott and Caves 
[2003]), while in (Giraud [2007a]) we find the exact formula for the n-th generic 
moment (and cumulant), and the explicit expression for the first five cumulants. 
The probability distribution for the nonnull eigenvalues of the reduced density 
matrix oi A or B have been computed for the case in which the dimensions of 
the corresponding Hilbert spaces T-La and Hb are dim'H^ = 3 < dimUs — 4 
and dimHA — dim'Hs = 4, in (Giraud [2007b]). In this thesis we will intro- 
duce a different and more general approach for studying the statistics bipartite 
entanglement, relying on the techniques of classical statistical mechanics (Facchi, 
Marzolino, Parisi, Pascazio and Scardicchio [2008]; De Pasquale, Facchi, Parisi, 
Pascazio and Scardicchio [2010]). In particular, in Chap. 2 we will introduce a 
partition function for the canonical ensemble, the role of the energy being played 
by the purity of subsystem A, chosen as a measure of the quantum correlations 
between A and B. The partition function will depend on a Lagrange multiplier 
/3 which corresponds to a fictitious inverse temperature, which fixes the value of 
the average purity. In other words, being (3 the conjugate variable of the purity 
it deforms the Haar measure, localizing it on the corresponding manifold of states 
with a given average entanglement, whose associated uncertainty becomes smaller 
as the dimension of the quantum system increases. We will prove the equivalence, 
for our system, between the canonical ensemble and microcanonical description 
in terms of isoentangled manifolds (Sinolecka, Zyczkowski and Kus [2001]). We 
will therefore explore the entire Hilbert space, from typical to maximally entan- 
gled and separable states, corresponding to different regimes of temperatures, see 
Chaps. 3 and 5. In Chap. 4 we will describe a metastable solution for the system 
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in the region of negative temperatures. In particular, by computing the dominant 
contribution of the partition function thru the method of steepest descent and 
minimizing (or maximizing) the free energy of the system for positive (or negative 
temperatures), we will determine the probability distribution of the Schmidth co- 
efficients corresponding to a given value of average entanglement, and not only for 
(3 — (typical states). Of course, once we know the partition function we can also 
compute the moments of the purity distribution and we will prove the consistency 
of our results with the cumulants' analysis provided in the papers cited above. We 
will unveil the presence of three main regimes in the system, the separable, the 
typical and the maximally entangled phases, separated by first and second order 
phase transitions. Finally in Chap. 6 we will overview our results abandoning the 
temperature and using the purity as our physical variable. One of the main results 
of our analysis will be the computation of the volume of the isopurity manifolds in 
the Hilbert space, that is the probability distribution of the purity. In particular 
this volume will be maximum for typical states, while the iso-entanglcd manifolds 
given by maximally entangled and separable states shrink to a vanishing volume. 
The latter case has been predicted in (Zyczkowski, Horodecki, Sanpera and Lewen- 
stein [1998]) for the more general case of arbitrary mixed states, by proving the 
existence of a topological lower bound for this volume. 

In the second part of the thesis, we will extend the above analysis to the case 
of a large quantum system in a mixed state, with a fixed value of the global purity. 
A statistical mixture describes the lack of complete knowledge of our physical sys- 
tem, when due to the interaction with the environment, it undergoes a dissipative 
evolution. In general, a complete microscopic description of the dynamical evolu- 
tion of a system coupled to the environment (or bath) is a complex many-body 
problem which requires the solution of a potentially infinite number of coupled dy- 
namical equations. According to an open system approach, this issue is tackled by 
retaining only basic information about the environment and describing the system 
dynamics in terms of a master equation (Petruccione and Breuer [2002]; Gardiner 
and Zoller [2004]). The lack of a complete knowledge about the bath leads to 
master equation coefficients (MECs), which may be either unknown or obtained 
from a microscopical derivation carried out within some approximation scheme. 
Some tomographic approaches, based on the idea of studying the evolution of a 
Gaussian probe in order to provide the sought-for bridge between dynamical pa- 
rameters and measurable quantities (the probe's cumulant) have been proposed in 
order to retrieve the above coefficients or checking the approximation schemes for a 
wide class of Gaussian Shape Preserving master equations, both in the Markovian 
and convolution less non-Markovian evolution (Bellomo, De Pasquale, Gualdi and 
Marzolino [2009]; Bellomo, De Pasquale, Gualdi and Marzolino [2010a]; Bellomo, 
De Pasquale, Gualdi and Marzolino [2010b]). If we deal with isolated quantum 
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systems, their density operator reduces to a rank one projection, and the analysis 
of its quantum correlations can be performed with the techniques introduced in 
Chaps. 2-6. The general case of arbitrary mixed states is far more involved than the 
case of pure states studied in the first part, due to more than one reason. The first 
one is related to the definition of a random sampling for the system. While for the 
case of pure states the only natural choice is obtained by requiring invariance under 
the full group of unitary transformations, for the case of mixed states, if we impose 
this unitary invariance we are just stating that the probability density depends only 
on the eigenvalues, while on the other hand we do not have a unique natural mea- 
sure on the simplex of eigenvalues (Wootters [1990]). It follows that, since no single 
distinguished probability measure exists, we can introduce many ensembles for ran- 
dom density operators. One possibility is given by adding an ancillary system of 
the same dimension, then determine a purification of our density operator, and 
consider the measure induced on it, by partial tracing, for instance from the Haar 
measure on the unitary larger group (Braunstein [1996]; Hall [1998]; Zyczkowski 
and Sommers [2001]). This strategy will be followed in Chap. 7. Other sampling 
techniques are based on the Bures distance (Hall [1998]; Slater [1999]; Bengtsson 
and Zyczkowski [2006]) or, alternatively on the multi-partite systems (Collins, Ne- 
chita and Zyczkowski [2010]; Zyczkowski, Penson, Nechita and Collins [2010]). The 
absence of a unique measure for the set of density matrices is an extremely deli- 
cate point. Indeed, there are some properties which strictly depend on the chosen 
measure, such as the volume of separable states, while other statistical proper- 
ties, such as the relation between the probability of entanglement as a function of 
the purity of the global state, are not very sensitive to the measure (Zyczkowski 
[1999]). Another aspect that makes the generalization of the approach introduced 
for the set of pure states very hard, is related to fact that from the definition of 
statistical mixture it follows that the purity of one subsystem (the local purity 
of the global system) cannot be considered a bonafide measure for the quantum 
correlations between subsystems A and B, but should be substituted by its con- 
vex roof (Wootters [2001]; Amico, Fazio, Osterloh and Vedral [2008]; Horodecki, 
Horodecki, Horodecki, and Horodecki [2009]). We will discuss this point at the end 
of Chap. 1. However, the above substitution does not allow for a simple analytic 
treatment of the partition function. We will thus devote the second part of the 
thesis to the study of the statistical distribution of the A-local purity, which can 
be considered as a lower bound for the bipartite entanglement between two parts 
of the global system. By generalizing the partition function introduced for the case 
of pure states, we will compute the canonical moments for the local purity, follow- 
ing different techniques. In particular in Chap. 8, we will introduce the Gaussian 
approximation for the elements of the unitary groups on the set of density ma- 
trices describing our system and on their purification, and compute the average 
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A-purity for /3 = up to order 0{1/N), being dim'HA — N. In Chap. 9 we will 
compute the exact expression for the first moment of the local purity at ^ = 0, 
by exploiting the properties of the twirling transformations (establishing a formal 
connection between our statical problem and the theory of quantum channel), and 
then, thru the solution due to Zuber of some basic integrals over the unitary group, 
we will determine the exact expression for the n-th generic moment, and explicitly 
compute the first two cumulants for the local purity at ^ = 0, that is for typical 
random states. We will therefore generalize the exact expressions computed for the 
pure case in (Scott and Caves [2003]; Giraud [2007a]). Finally, we will compute the 
high temperature expansion of the first moment of the local purity in the canonical 
ensemble. 

As a final remark, notice that the classical statistical mechanics approach we 
will follow in this thesis in order to characterize the behavior of bipartite entan- 
glement for a large quantum system is general, and for instance can be adopted 
also for determining the statistical distribution of the potential of multipartite 
entanglement, as shown in (Facchi, Florio, Marzolino, Parisi and Pascazio [2010]). 



Chapter 1 



An introduction to bipartite 
entanglement 




Figure 1.1: "Entanglement: folding the arms" 



The origin of the word "entanglement" goes back to 1935. It is the English 
translation made by E. Schrodinger of the German word "Vershralung" , which 
is a colloquial expression indicating the action of "folding the arms", Fig. 1.1. 
Schrodinger introduced this name in the dictionary of quantum physics in order to 
describe one of the most striking manifestations of quantum phenomena, that is 
"the best possible knowledge of a whole does not include the best possible knowl- 
edge of its parts, even though they may be entirely separated" (Schrodinger [1935a]; 
Schrodinger [1935b]; Schrodinger [1936]). This is due to presence of correlations of 
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exclusively quantum nature, namely entanglement. In other words, the presence 
of quantum correlations, with no classical counterparts, is responsible for the fact 
that while for a classical bipartite system the knowledge of the total system is 
completely equivalent to the knowledge of both its subsystems, in the quantum 
case this is only a necessary condition. 

The phenomenon of entanglement is profoundly rooted in the mathematical 
structure of quantum mechanics, whose natural framework is given by complex 
Hilbert spaces. In particular, entanglement is a direct consequence of the linearity 
of the Hilbert space (expression of the superposition principle), when applied to 
composite systems, that arc described by the tensor product of the Hilbert spaces 
associated to each subsystem (Horodecki, Horodecki, Horodecki, and Horodecki 
[2009]; Facchi [2009]). 

The peculiarities of this quantum phenomenon make the entanglement an ex- 
tremely important resource in quantum information processing. It is believed to 
be the main responsible for the speed up in quantum computation and quantum 
communication (Nielsen and Chuang [2000]; Benenti, Casati and Strini [2004]). 
Furthermore, many quantum protocols could not be conceived without the exis- 
tence of entangled states, an outstanding example is given by quantum teleporta- 
tion (Bennett, Brassard, Crpeau, Jozsa, Peres, Wootters et al [1993]). Besides its 
important applications in relatively simple systems, that can be described in terms 
of a few effective quantum variables, entanglement has also been widely investi- 
gated in many-body systems (Amico, Fazio, Osterloh and Vedral [2008]; Horodecki, 
Horodecki, Horodecki, and Horodecki [2009]), from many points of view. In partic- 
ular, on the one hand, the methods used in quantum information have proved to be 
extremely useful for the control and the manipulation of this kind of systems, on 
the other hand the interplay between quantum information and quantum statistical 
mechanics has unveiled some peculiar properties of entanglement in quantum crit- 
ical models (Osterloh, Amico, Falci and Fazio [2002]; Osborne and Nielsen [2002]; 
Vidal, Latorre, Rico and Kitaev [2003]). As a consequence of the powerful role 
played by entanglement in the study of fundamental quantum phenomena and its 
applications to innovative technologies, the task of characterizing and quantifying 
entanglement has emerged as a prominent theme of quantum information theory. 

This chapter is devoted to the analysis of some fundamental aspects of bipartite 
entanglement. We will show that there are many basic open questions, revealing 
that we are still far from a complete comprehension of this phenomenon. We will 
start from the definition of entangled state in Sec. 1.1. Then, in Sec. 1.2, we will 
focus on the only class of states for which the basic properties of bipartite entangle- 
ment, together with an appropriate measure, can be considered well understood, 
namely the pure quantum states. Finally, in Sec. 1.3 we will show in what sense 
the study of entanglement and its quantification become extremely involved for an 
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arbitrary quantum state, and report some fundamental results obtained in the last 
decades. 

1.1 Bipartite entanglement 

E. Schodinger defined entanglement as "not one but rather the characteristic trait 
of quantum mechanics" . As already mentioned, this phenomenon is strictly related 
to the mathematical formulation of quantum mechanics, relying on the structure 
of complex Hilbert spaces. 

The smallest non trivial Hilbert space is two dimensional, and describes the 
simplest quantum system, the qubit. This is the fundamental unit of quantum 
information, as the bit is the indivisible unit of classical information. A privileged 
orthonormal basis of this two dimensional Hilbert space is the computational basis, 
{|0),|1)}, which can be seen as the quantum version of the only two possible 
values associated to a single bit, {0, 1}. The way a qubit differs from a bit is 
that a hnear combination of the states |0) and |1) is still a possible state for 
the qubit. This derives from the vectorial structure of the Hilbert space, and 
represents the simplest expression of one of the cornerstone of quantum mechanics, 
the superposition principle. 

Let us now consider a composite bipartite system, and indicate with A and B its 
two parts. In classical physics, the states of the total system belong to the Cartesian 
product of the spaces associated to A and B, and are always given by a convex 
combination of products of states describing each subsystem independently. The 
dimension of the global system is then the sum of the dimensions of its subspaces. 
On the contrary, the Hilbert space T-L where the states of the total system live is 
given by the tensor product of the Hilbert subspaces T-La and Hb associated to 
A and B respectively, H — Ha <8) Hb, and the dimension of H is the product of 
the dimensions of Ha and Hb- If for instance, we refer to a system of n qubits, 
we can study the properties of two subsystems given by Ui and n2 = n — Ui 
qubits. The associated Hilbert spaces have dimensions 2", 2"^ and 2"^^ respectively 
(2" = 2"'i+"2 — 2"i2'*2). Let us now come back to a generic bipartite system, and 
to the difference between the classical and quantum approach. It is again the 
superposition principle which plays the key role. Indeed an arbitrary pure state 

G H cannot in general be expressed by the product of separate states of each 
subsystem, as would happen in the classical case: this is the formal definition of 
the entanglement. In other words, we will define the state |^) entangled, with 
respect of the bipartition [A, B) if and only if 
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for some vectors l'^)^ and lip)^ of each subsystem. If condition (1.1) is not satisfied, 
the state will be said to be separable, or unentangled. Let us remark that the 
separability of the state of the global system strictly depends on the bipartition 
we are considering, as it refers to the quantum correlations between the specific 
subsystems. From Eq. (1.1) it follows that if a pure state l^^) is entangled it is 
given by a superposition of at least two product states: 



This describes interference, i.e. non classical correlation, among probabihty am- 
plitudes for the two subsystems. For example, the probability amplitude of having 

subsystem A in the state and subsystem B in the state \ipj.2)B interferes 

with the probability amplitude of having subsystems A and B in states \'ipki)A 
and \'ipk2)B- ^'^ nonclassical properties of quantum states that many 

schemes of quantum information and quantum technologies (quantum computation 
(Nielsen and Chuang [2000]), quantum teleportation (Bennett, Brassard, Crpeau, 
Jozsa, Peres, Wootters et al [1993]), dense coding (Bennett and Wiesner [1992]) 
and quantum cryptography (Bennett and Brassard [1984]; Ekert [1991]; Deutsch, 
Ekcrt, Macchiavcllo, Popcscu and Sanpera [1996]; Fuchs, Gisin, Griffiths, Niu and 
Peres [1997])) have been realized. 

Until now we have only referred to pure states, which correspond to normalized 
vectors of a given Hilbert space. They describe isolated quantum systems whose 
states arc completely determined as far as the theory allows. However, in practice 
the state of a physical system is not often perfectly determined. For instance in 
a laboratory the systems we measure undergo uncontrolled interactions with the 
environment. An example is given by the emission of atoms from a thermal source: 
we only know the distribution of the kinetic energy of the emitted particles, not 
the kinetic energy of each of them. In such cases, we say that our knowledge of the 
system is incomplete. We only know that the state of the system belongs to a given 
ensemble of pure states associated to a set of probabilities, {pi, . . .pn} satisfying 
the condition of unit total probabihty, ^^Pi — 1. The system is then defined to 
be in a statistical mixture, or equivalently in a mixed state. From a mathematical 
point of view, an arbitrary quantum state is described by a density operator p on 
% (a positive, self-adjoint and unit trace operator), which in case of a pure state 
IV') e % reduces to the projection operator p = A generic mixed state 

can be expressed in an infinite variety of convex combinations of pure states all of 
which have exactly the same consequence for any conceivable observation of the 
system. It naturally follows that entangled mixed states arc no longer equivalent 
to non-product states, as for pure states. A mixed bipartite state is separable if 
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and only if it can be written as a convex combination of product states: 



P= ^ Pi PA,i ® pB,i , Pi > 0, ^ Pi = 1, n > 1, 



(1.3) 




where pA,i and pB,i are in general mixed states of the corresponding subsystems 
A and B. If condition (1.3) is not satisfied, p is said to be an entangled state. 
It can be easily proved that entanglement is responsible for those properties of 
composite systems which do not change under local transformations and classi- 
cal communication (exchange of classical bits between the two subsystems). We 
say that bipartite entanglement is invariant under LOCC, Local Operation and 
Classical Communication. This has led to the introduction of some criterions in 
order to properly identify classes of equally entangled states (Grassl, Rotteler and 
Beth [1998]; Linden, Popescu and Sudbery [1999]; Albcverio and Fei [2001]), and 
to the complementary study of the orbits in the Hilbert space associated to a given 
quantum state by means of local unitary operations (Sinolecka, Zyczkowski and 
Kus [2001]). 

Despite the simplicity of definitions (1.3) and (1.1), it is in general quite hard 
to check whether a state is separable or not, with respect to a given bipartition. 
Another open problem is the quantification of entanglement. In Sec. 1.2 we will 
answer these two points for the case of pure states, and in Sec. 1.3 we will then 
show how this analysis becomes far more involved for the case of mixed states. 



Bipartite pure states represent one of the few cases for which the problem of re- 
vealing and also quantifying quantum correlations has a quite exhaustive answer. 
In the first part of this section we will discuss some fundamental properties of 
bipartite entanglement, by the introduction of one important tool, the Schmidt 
decomposition. In the second part, we will consider the problem of quantification 
of bipartite entanglement through the von Neumann entropy. 

1.2.1 Schmidt decomposition 

Let us consider a pure bipartite state l'^) in ?^ = Ha ® V-b- We introduce 
the following notation for the dimension of the three Hilbert spaces considered: 
dim?^^ = iV, dimT/s = Af and dim'H — L — NM. Without loss of generality we 
also set N < M. If we introduce an orthonormal product basis {|e„)^ km)^}) 
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1 <n < N and 1 <m < M, we can represent the state as 

IV') = ^ Cij\^i)A<»\(^j)B- (1-4) 

l<i<Ar l<j<M 

The coefficients Cij are the elements oi em N x M complex rectangular matrix C, 
whose singular values \/Xk, 1 < k < N {i.e. the square root of the eigenvalues of 
the matrix C^C) determine the Schmidt decomposition of the total state 

l^)= E y^^\k)A®\k)B- (1-5) 

l<k<N 

It is uniquely determined only if the singular values are not degenerate. Indeed, 
the previous expression can be derived from the singular value decomposition of 
the coefficient matrix, namely C — UDV^ , with U and V unitary matrices on Ha 
and T-Lb, respectively, and D the diagonal matrix D = diag{-\/Ai, ■ ■ ■ , V^n}- The 
unitary matrices U and V, determine the orthonormal set 

\k)A^\k)B= Uki\ei)^®Vkj\ej)j^, yke{l,...N}. (1.6) 

l<i<iV l<j<M 

If the singular values of C are not degenerate its singular value decomposition is 
unique, up to multiplication of one or more columns of U by unit phase factors 
and simultaneous multiplication of the corresponding columns of V by the same 
unit phase factors. On the other hand, degenerate singular values do not have 
unique singular vectors. Consequently, if there arc degenerate ^ ^ 'i ^ the 
singular value decomposition of C is not unique. The real numbers Aj are called 
Schmidt coefficients and, due to the normalization condition on the state of the 
global system, (V'lV') = 1) obey Aj = 1 (notice that in the literature instead 
of Aj, 1 < 2 < A^, arc usually defined as Schmidt coefficients). The set of all possible 
vectors A = {Ai, . . . , Aat} forms a (A^ — 1) dimensional simplex A^, known as the 
Schmidt simplex. The Schmidt decomposition is of central importance for the 
characterization and the quantification of the entanglement associated with pure 
states. It enables to investigate many important properties of a quantum composite 
system by looking at the behavior of its subsystems, through the analysis of the 
corresponding reduced density operators. As it will be more clear at the end of 
this section, this is strongly related to the very essence of entanglement. The 
state of each subsystem, A (or B), is described by the density operator pA (or ps) 
defined by partial tracing over the complementary subsystem, pA — Tr b{\iP) 
(or Pb = Tr^i {ip\)), where Tr^ indicates the partial trace over subsystem S (this 
definition of a reduced density operator is also valid for an arbitrary state p) . From 
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Eq. (1.5) we immediately get that pA and pb have the same nonzero eigenvalues, 
the Schmidt coefficients, and their eigenbases coincide with the orthonormal bases 
and 1 < i < A^, called Schmidt bases for "H^ and Hb- The number of 
nonvanishing Aj is equal to the rank of both reduced density operators. It is called 
Schmidt number for the state and together with the Schmidt coefficients, 
is invariant under local unitary operations: they will be the key ingredients for 
quantifying the entanglement between subsystems A and B. Indeed, a bipartite 
pure state is separable with respect to the bipartition {A, B) if and only if there is 
only one nonzero Schmidt coefficient, which must be equal to 1. This is equivalent 
to saying that each subsystem is in a pure state. If it is not the case, the state 
is entangled. In particular, a pure state G "H is called maximally entangled if 
and only if all its Schmidt coefficients are equal to 1/A^, or equivalently subsystem 
A is in a completely mixed state. It is worth noticing that we have introduced 
the expression "maximally entangled state" without referring to any measure of 
entanglement. Indeed these states present the maximal degree of entanglement 
with respect to every (reasonable) measure. They are a very peculiar class of 
quantum states, and can be considered the greatest manifestation of the presence 
of quantum correlations between subsystems. More formally, they give rise to the 
most striking violation of a class of inequalities, Bell inequalities, that must be 
satisfied by classical correlations. If a bipartite system is in a maximally entangled 
state, we have that to a perfect knowledge of the state of the composite system 
(1-0) is a pure state) it corresponds the maximal ignorance of the state of its two 
parts {pa is in a completely mixed state), since the global information is completely 
shared between the two subsystems. This situation never happens in the classical 
case, where the complete knowledge of the state of the global system and of the 
state of both subsystems is equivalent. In quantum mechanics this represents 
just a necessary condition, that becomes sufficient only for the case of separable 
(unentangled) bipartite states, when both the global system and its subsystems A 
and B are in a pure state (Facchi [2009]). This behavior, is the most impressive 
signature of quantum correlations, i.e. entanglement. 

1.2.2 von Neumann entropy and purity 

Bipartite entanglement is a very complex property of quantum states, and the task 
of capturing all its manifestations with an appropriate measure is in general very 
hard. Until now is has not been found a unique answer to this question. Never- 
theless, pure bipartite states represent one of the few cases for which this intricate 
feature simplifies very much. The Schmidt decomposition is an elementary nec- 
essary and sufficient criterion for separability of pure bipartite states. In the last 
section we have seen that the Schmidt coefficients, and obviously their number, are 
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invariant under local unitary operations on the system. Therefore, any reasonable 
measure of entanglement for pure bipartite states must be an appropriate func- 
tion of these quantities. A widely accepted measure of entanglement is the von 
Neumann entropy of the reduced density operators pa and pb describing the two 
subsystems A and B of the global pure state l^^): 

S{pa) = S(pb) = -Tr(pAlog^PA) = - ^ Aaog(A,), (1.7) 

l<i<7V 

Aj, 1 < i < A^, being the Schmidt coefficients of lip), and OlogO defined to be 
zero. The von Neumann entropy goes from for separable states to log2 N for 
maximally entangled states. It is the unique entanglement measure for pure states 
that beside being invariant under LOCC is also continuous and additive if there 
are several copies of the system (5'(|'i/'i) (8) ... (8) |V'n)) — Z]j<n '^(I'^j)))- These 
properties are some of the requirements that, according to an axiomatic approach, 
a good measure of entanglement has to satisfy. We refer to Sec. 1.3.2 for a brief 
review about some peculiar aspects of entanglement quantification. 

Historically, the von Neumann entropy was the first entanglement measure to 
have an operational interpretation: it quantifies the quantum information of the 
system, i.e. the deeply nonclassical interplay between the information provided by 
an entangled state about the whole system and the information regarding its sub- 
systems. In 1995 Schumacher showed that the von Neumann entropy corresponds 
to the minimum number of qubits needed to encode a quantum state produced by 
a statistical source, in an ideal coding scheme (Schumacher [1995]). In this sense it 
is exactly the quantum counterpart of the classical Shannon entropy which counts 
the minimum number of bits (units of classical information) one needs to encode 
the output of a random source using an ideal code. As already pointed out at the 
end of Sec. 1.1, the main difference between classical and quantum information, is 
that if a state is entangled with respect to a given bipartition, it provides more 
information about the global system then about its two parts, whereas this can 
never happen for a classical state. Translated in terms of statistical mechanics, 
the subsystems of an entangled state can exhibit more entropy than global system, 
i.e. they can show more disorder than the system as a whole. This never hap- 
pens in two cases: when we have a classical system (whose statistical disorder is 
quantified by the Shannon entropy) and if we deal with a quantum separable state. 
Once again, entanglement manifests itself as the distinctive character of quantum 
mechanics. 

The definition of the von Neumann entropy as a good candidate for measur- 
ing bipartite entanglement in case of pure states is supported by the so-called 
"uniqueness theorem" (Vidal [2000]; Horodecki, Horodecki and Horodccki [2000]; 
Nielsen [2000]; Donald, Horodecki and Rudolph [2002]) which states that in the 
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asymptotic regime all other entanglement measures, introduced on the basis ei- 
ther of quantum information or of thermodynamical considerations, coincide on 
pure bipartite states and reduce to the von Neumann entropy of the correspond- 
ing reduced density matrices. A more convenient measure of entanglement from a 
computational point of view is the purity ttab- It corresponds (up to a constant) 
to the so-called linear entropy Cab, that is the first-order term of the expansion of 
the von Neumann entropy: 

S(pa) = -Tr(pA log2 pa) ~ -Tr(pA(pA - 1)) = 1 - Tr{p\) = ^^JCab (1-8) 
and the purity of one of the two subsystems is defined as: 

TTAB = IV(pi) = Tr{pl) = J2 (1-9) 

l<i<N 

We can easily see that the purity lies in the compact interval [l/A?", 1], whose 
boundaries from left to right refer to maximally entangled and separable states, 
with respect to the bipartition {A, B) (the corresponding bounds of the linear 
entropy are 1 and 0). In this thesis the purity will be the fundamental tool in order 
to investigate the statistical behavior of bipartite entanglement of pure quantum 
states. 



1.3 Bipartite mixed states 

Beside, and maybe "before" , the problem of properly quantifying entanglement, an 
open question of quantum information theory is to determine whether an arbitrary 
(mixed) quantum state is separable or not with respect to a given bipartition, 
namely the separability problem. We have already stressed that pure bipartite 
states represent a singularity in this complex scenario since the Schmidt decom- 
position on the one hand provides a simple necessary and sufficient criterion for 
checking the separability of pure states, on the other hand it naturally leads to the 
quantification of entanglement in terms of the mixedness of its subsystems. On 
the contrary, the complete characterization of the convex set of separable mixed 
states reveals itself to be extremely involved (Zyczkowski [1999]). We will devote 
the present section to a brief overview of the main separability criteria and mea- 
sures of bipartite entanglement for a generic mixed state (Amico, Fazio, Osterloh 
and Vedral [2008]; Horodecki, Horodecki, Horodecki, and Horodecki [2009]; Bruss 
[2002a]). 
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1.3.1 Separability criteria 



A very strong necessary separability condition for quantum bipartite states is 
given by Peres' positive partial transpose (PPT) criterion (Peres [1996]; Horodecki, 
Horodecki, and Horodecki [1996]). It says that if a quantum state p is separable, 
the state obtained by partial transposition is positive. Partial transposition indi- 
cates the application of the partial transpose operator to only one of two subsys- 
tems, namely {p'^^)m^i,ny = ip)nfi,mu, where the Latin indices refer to subsystem 
A (m, n = 1, . . . A^) and the Greek ones to subsystem B {pi^v = 1, . . . M). The 
same holds for partial transposition with respect to subsystem B. This condition 
which is necessary for arbitrary dimensions of TIa and T-Lb, is also sufficient for low 
dimensional systems, 2x2 and 2x3. The PPT criterion, is an "operational recipe" 
for checking the separability of a given state. Other operational approaches are the 
reduction (Horodecki and Horodecki [1999]) and the majorization criteria (Nielsen 
and Kempe [2001]). From these effective techniques, especially from the PPT con- 
dition, a more general technique in terms of linear maps has been derived. However, 
even if linear maps criteria provide both necessary and sufficient conditions to be 
satisfied by any bipartite separable system, they do not supply any simple proce- 
dure for determining if a bipartite state is separable. In this sense they can be 
defined "non operational" (Bruss [2002a]). The positive maps criterion says that a 
quantum state p is separable with respect to a given bipartition {A, B) if and only 
if for any positive maps and A^ one has that (A^ ® l)p > and (1 ® Ab)p > 0. 
This condition is, automatically satisfied by the PPT criterion, with Ka/b — Ta/b- 
However, the positive maps criterion is nonoperational, since its application would 
require a complete knowledge of the set of all positive maps, and this is itself an 
open problem. As special case of positive maps, there is the criterion of the so-called 
entanglement witnesses. It says that a quantum state p is entangled if and only if 
there exists an Hermitian operator W such that Tr(l^p) < and Tr(l^p*^*') > 
for any separable state p^^^. However, if the expectation value oiW oia. p is non 
negative, this does not guarantee that p is separable. This can be clarified intro- 
ducing the following geometric interpretation: each entanglement witness defines 
an hyperplane in the Hilbert space Hab that separates p from the set of separable 
states. This geometric approach also helps to understand how to optimize W: one 
can perform a parallel transport of the hyper-planes (Lewenstein, Kraus, Cirac 
and Horodecki [2000]) or replace them by curved manifolds (Ghne [2004]) both 
tangent to the set of separable states, see Fig. 1.2. Entanglement witnesses are ex- 
plicitly related to positive maps by the Jamiolkowski isomorphism (Jamiolkowski 
[1972]), according to which one has that given an entanglement witness W there 
exist a positive map A, such that W = (1 ® A)P, where P is the projection op- 
erator on maximally entangled states (P = l/A^AX^fj=i N) Oil))- Observe that 
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notwithstanding these optimization techniques, the problem of finding a complete 
characterization of the convex set of separable states (Zyczkowski [1999]) is still 
open, since it would need in principle an infinite number of entanglement witnesses. 



W 




Figure 1.2: Intuitive picture of an entanglement witness W (continuous straight 
line) and its optimizations (dashed lines). 



1.3.2 Entanglement measures 

Quantifying entanglement is one of the central topics of quantum information the- 
ory. As already mentioned, entanglement represents an extremely complex prop- 
erty of quantum states and for arbitrary states there is not a unique definitive 
measure. In Sec. 1.2.2, we have seen that in the asymptotic regime, when we 
consider the tensor product states of a large number of identical copies of the 
system, due to the uniqueness theorem, the entanglement finds a sort of regular- 
ization, at least on pure states. Following an axiomatic approach, there exists a 
set of postulates that a good measure of entanglement should fulfill: monotonicity 
under LOCC, vanishing on separable states, continuity, additivity, subadditivity 
and convexitiy (Horodecki, Horodecki, Horodecki, and Horodecki [2009]; Amico, 
Fazio, Osterloh and Vedral [2008]; Bruss [2002a]). However, there are contrast- 
ing believes about the necessity of all these requirements for properly quantifying 
entanglement, and furthermore it is still not known if there exists a measure sat- 
isfying all of them (Vedral, Plenio, Rippin and Knight [1997]; Vidal [2000]). A 
natural question may also be whether, despite the absence of a unique measure of 
entanglement, we can order all the states of a given bipartite system with respect 
to the degree of entanglement (Virmani and Plenio [2000]). The answer is negative. 
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since there exist incomparable states even in the case of pure states. This lack of a 
single order can be explained, at least qualitatively, in terms of the complexity of 
entanglement, in the sense that there are many manifestations of entanglement and 
it may happen that one state shows more entanglement of one type than another 
state, and the situation is reversed for another type of entanglement (Miranowicz 
[2004]; Verstraete, Porras and Cirac [2004]). 

Two very important measures of quantum correlations are the entanglement of 
distillation Ed, defined as the maximal number of singlets that can be produced 
from a given bipartite quantum state by means of any LOCC operation, and the 
entanglement cost Eq, which is complementary to E^, and measures the minimal 
number of singlets needed to produce a given quantum state by LOCC operations. 
We recall that the singlet state is a particular maximally entangled state for a qubit, 
whose expression in the computational basis is l/\/^(|01) — |10)) (for simphcity, 
we indicate the tensor product \a)^ ® as \af3)). It has been proved that the 
entanglement of distillation and the entanglement cost are, respectively, the lower 
and upper bounds of any entanglement measure satisfying appropriate postulates 
in the asymptotic regime (Horodecki, Horodecki and Horodecki [2000]). Both Ed 
and Ec fulfill the additivity axiom, but their continuity has not been completely 
proved. There are also indications that the entanglement of distillation is not 
convex (Shor, Smolin and Terhal [2001]). 

A technique which enables to extend entanglement measures defined on pure 
states to arbitrary states is the convex-roof method. It says that if is a mea- 
sure defined on pure states, its natural extension to a mixed state p is E{p) — 
min^-PiE{\i/ji))j where the minimum is taken over all possible convex combina- 
tions of pure states, p = YlkPk iV'fc) (V'fclj being E an increasing function of bipartite 
entanglement. In this way it is defined another important measure of entanglement, 
the entanglement of formation: 



being p^^k = TTB{\ipk) {^kD- It is the convex-roof extension of the von Neumann 
entropy (Wootters [2001]). This measure is continuous and convex but its full 
additivity for bipartite systems has not been established yet (Vidal, Diir and Cirac 
[2002]). Finally, by analogy with pure states, a measure of entanglement more 
convenient for analytical treatment is the convex-roof extension of the purity of 
one of the two subsystems: 




(1.10) 




(1.11) 
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Note that in this case we have to maximize over the mixture of pure states corre- 
sponding to p, since the purity is a decreasing function of the bipartite entangle- 
ment. 

The main topic of this thesis will be the study of the statistical properties of 
bipartite entanglement, shared by large bipartite quantum systems. A complete 
overview will be discussed for the case of bipartite pure states, as anticipated at the 
end of Sec. 1.2.2, by considering the purity as the tool for measuring entanglement. 
The generalization of this analysis to the case of bipartite mixed state, as might 
be expected, will reveal far more complicated then the previous one. 



Part I 
PURE STATES 
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Summary 



In the first part of tliis tliesis wc will characterize the statistical properties of bipar- 
tite entanglement for a large quantum system in a pure state. This problem will be 
tackled by studying a random matrix model that describes the distribution of an 
entanglement measure, the purity of one of the two subsystems (Facchi, Marzolino, 
Parisi, Pascazio and Scardicchio [2008]; De Pasquale, Facchi, Parisi, Pascazio and 
Scardicchio [2010]). In particular, we will introduce a partition function for the 
canonical ensemble as a function of a fictitious temperature. The purity will play 
the role of the energy of the system, that is different temperatures will correspond 
to different degrees of entanglement. The outcome of this statistical approach 
to bipartite entanglement will be twofold. On the one hand we will distinguish 
three main entanglement phases (or regions): maximally entangled, typical and 
separable phase. These regions are separated by a first and a second order phase 
transition. On the other hand, by the explicit computation of the entropy of the 
system, we will evaluate the volume of the manifolds in the Hilbert space with 
constant purity (isopurity manifolds), answering to the question of determining 
how probable is finding an entangled state, for a large quantum system. 

This part of the thesis is organized as follows. In Chap. 2 we will introduce 
the notation and set the bases of our classical statistical mechanics approach to 
the problem. In Chap. 3 we will study the case of positive temperatures, where at 
very low temperatures we will find very entangled states. Negative temperatures 
will be investigated in Chaps. 4 and 5 where we will show the existence of two 
branches, a stable one associated to a partial factorization of the state, and a 
metastable branch which contains the 2D quantum gravity point. We will also 
investigate finite size corrections. Finally, in Chap. 6 we will overview the results 
of the previous chapters, by reinterpreting them directly from the point of view of 
quantum information. 
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A statistical approach to bipartite 
entanglement 

In this chapter we will introduce the bases of our approach to the study of bipartite 
entanglement for a large quantum system. The first requirement will be setting 
the statistical ensemble of pure states, and this will led to the introduction of 
a random matrix model for the system in Sec. 2.1. The next step will be the 
definition, in Sec. 2.2, of the fundamental tool for this statistical approach: the 
partition function of the system for the canonical ensemble, where the role of the 
energy will be played by the entanglement measure we will choose, the purity of 
one of the two subsystems. Since the exact computation of the partition function 
results very complicated due to the positivity constraint on the eigenvalues, we will 
recur to the method of steepest descent and impose a set of saddle point equations 
for the Schmidt coefficients. 

2.1 Random matrix model and Haar measure 

Consider a bipartite system {A,B). In Sec. 1.1 we have seen that the states of 
the total system live in the tensor product Hilbert space H = T-La <S)'Hb- Without 
loss of generality we set dim'Hyi = N < dim'HB = M. We also assume that the 
system is in a pure state e "H. A statistical analysis of bipartite entanglement 
is equivalent to answer the following question: what is the probability of finding 
a state with a given value of entanglement in the Hilbert space T-Ll In order to 
answer this question, we have to solve two preliminary points: 

• introducing a proper measure of entanglement 

• defining a probability measure on T-L, according to which the states are sam- 
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pled. 

The answer to the first point is almost straightforward. In Chap. 1 we have seen 
that for the case of pure states the von Neumann entropy of one of the two subsys- 
tems, S{pa) = S{pb) = — Pa log Pa, is a good measure of bipartite entanglement. 
However, at the end of Sec. 1.2.2 we have also noticed that the von Neumann 
entropy can be replaced for computational convenience by a more simple function, 
the purity, which derives from its linearization: 

TV 

TTAB = Tr(p^) = Tr(p|) = ^ ttab e [1/N, 1]. (2.1) 

i=i 

The minimum is attained when all the eigenvalues Xj, the Schmidt coefficients of 
the global state \ip), are equal to 1/A^ (subsystem A in a. completely mixed state 
and maximal entanglement between the two bipartitions) , while the maximum is 
attained when one eigenvalue is 1 and all others are (this detects a factorized 
(unentangled) state). Let us now switch to the second point. The idea of taking a 
quantum state randomly is equivalent to assuming the minimal knowledge about 
the system. Random states can be considered typical in the sense that they rep- 
resent the states to which an arbitrary evolving quantum state can be compared. 
For the case of pure states there is a unique unbiased measure, as it is enough to 
identify the minimal knowledge with the maximal symmetry for the system. Thus 
the sampling criterion has to be invariant under the full group of unitary trans- 
formations (Hall [1998]), and corresponds to the unique (left and right) invariant 
Haar (probability) measure djiuiU) on the unitary group U e U{l-L) ~ U{L): 

dpniU) = dpH{V), yU,Ve U{n). (2.2) 

We will define the typical vector states l^') (Lubkin [1978]; Page [1993]) as pure 
random states given by the action of a random unitary matrix U G IA{'H) on an ar- 
bitrary reference state G "H, = U \ipo), see Fig. 2.1. The final state {ip) will 
therefore be independent on l^^o)- Henceforth, all measures will be tacitly assumed 
to be probability measures (normalized to 1), if not stated otherwise. For exam- 
ple, ijLh{'U{H)) — 1. Nevertheless, it is important to add that the Harr measure is 
not the only possible sampling choice. Other distributions can also be considered 
but they encode additional information on the system (in this sense, Haar is the 
most neutral). These alternative distributions could be treated in our approach 
by constraining the system by means of Lagrange multiplier. For an approximate 
realization of this Haar measure by means of short quantum circuits see (Arrow 
and Low [2009]), where it is proved that one can extract Haar-distributed random 
states by applying only a polynomial number of random gates. 
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Figure 2.1: Typical states \ip) sampled according to the unique left and right 
invariant Haar measure on the unitary group W('H), from an arbitrary state \ipo)'- 
|^) = f/|^o),withf/GW(H). 

2.1.1 Induced measure on subsystems 

In order to describe the statistical properties of the purity of a bipartite system we 
need to determine the induced measure on the reduced density operator describing 
the state of one of the two subsystems, for instance pA- Representing pA in terms 
of its spectral decomposition pA = J2f ^tPi^ with {Pn}n=i,...,N a complete set of 
orthogonal projections defined up to a unitary rotation, we have that the set of 
the Hermitian reduced density matrices Sa can be written as a Cartesian prod- 
uct (Zyczkowski, Horodecki, Sanpera and Lewenstein [1998]; Zyczkowski [1999]): 

Sa = Pa^ Aa, (2.3) 

where Pa is the family of the complete sets of orthonormal projections on "H and 
is the Schmidt simplex. It follows that the measure on the reduced density 
operators pa is a product measure 

dp{SA) = dp{pA) = dfi{UA) X rf(7(A^), (2.4) 

where dp{UA) is defined on the unitary group U{1-La) — U{N), i.e. on the eigen- 
vectors of pai and dij{KA) on its eigenvalues. By construction the two measures 
are both induced by the Haar measure over the larger unitary group W('H) and are 
given by the Haar measure on U{1-La), dp{UA) = dpuiUA) and, for the eigenvalues, 
by (Lubkin [1978]; Page [1993]): 

dti{kA)=Ci,,M n (^-^-^3? n aS'^-^^^Ii- ^H^"^^- (2-5) 

l<i<j<N l<e<N \ l<k<N J 
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The square of the Vandermonde determinant {Y[i<i<j<Ni'^i ~ -^i)) derives from 
integrating out the eigenvectors, and the normahzation constant can be written in 
terms of the Euler gamma function r(x): 

m 

^"'^ - n„<,<«-,r(Ar-i + i)r(M-j)- 

It is worth observing that the measure over the eigenvalues in not uniform. In 
particular the Vandermonde determinant introduces a repulsion between pairs of 
Schmidt coefficients. 



2.2 The model 

In Sec. 1.1 we mentioned some studies related to the idea of determining iso- 
entanglement manifolds, through the introduction of invariants under LOCC (Grassl, 
Rotteler and Beth [1998]; Linden, Popescu and Sudbery [1999]; Albeverio and Fei 
[2001]), and by means of the complementary study of the orbits in the Hilbert space 
associated to a given quantum state (Sinolccka, Zyczkowski and Kus [2001]). The 
statistical approach presented in this thesis is slightly different, in the sense that we 
will describe manifolds with constant average purity. However, at the end of this 
chapter (Sec. 2.3) we will prove the equivalence between this approach and the one 
based on the microcanonical ensemble, in which the purity of distinct manifolds 
is fixed. Let us start with the introduction of a partition function (Facchi, Mar- 
zoUno, Parisi, Pascazio and Scardicchio [2008]) from which all the thermodynamic 
quantities, for example the entropy or the free energy, can be computed: 

ZabW) = J diiipA) exp {-^N'-t^ab) ■ (2.7) 

Recall that the measure d/j{pA) is a product measure, see Eq. (2.4). The parameters 
a and /3 in definition (2.7) arc extremely relevant for our analysis. The coefficient /3 
can be interpreted as the inverse of a (fictitious) temperature T selecting different 
regions of entanglement, ^ — 1/T. More formally, it is a Lagrange multiplier which 
enlightens different regions of the Hilbert space with a given value of average purity, 
(ttab)- In particular we have that: 

• for /3 = we obtain typical states (sampled by dp{pA)), 

• for /3 > we expect to find more entangled states, which in the limit ^ — >■ oo 
should become maximally entangled {{t^ab) = 1/-^)) 
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• for /3 < we are going towards the set of separable states, reached when 
P ^-oo {{-kab) = 1). 

In this scenario, the purity corresponds to the energy of the system, and varies 
in the interval [1/iV, 1] according to the different profiles of the distribution of the 
Schmidt coefficients. The parameter a is a scaling coefficient. It is a positive 
integer (either 2 or 3, as wc shall see) and refers to the scaling properties of the 
system. Its value needs to be chosen in order to yield the correct thermodynamic 
limit, that is in order to keep the exponent of the partition function extensive. 
Since the number of degrees of freedom of is A?"^ — 1 ~ N"^ (for a large system) , 
a needs to satisfy the constraint: 

N'^iTiAB) = O (7V2) . (2.8) 

Around maximally entangled states (for /3 > 0) we have {tvab) = O {1/ N) so a — 3, 
while around separable states (for /3 < 0) we have (tiab) = O (1) and hence a = 2. 

In the following wc will assume = M, since this docs not change the qualita- 
tive picture, the extension to N ^ M being straightforward but computationally 
cumbersome. 

2.2.1 Saddle point equations 

Prom the cycUc property of the trace, we have that the purity depends only on 
the eigenvalues of p^, thus the integral over the unitary group deriving from the 
measure on pA factorizes in the partition function, and wc have only to integrate 
on the (positive) Schmidt coefficients, modulo the constant factor Cjv,iv: 

Zab^ I d^X n (A.-A.f^A- 5^ A,')e-'^^"2:.,.,.Af ^3.9) 

JXi>0 l<i<j<N \ l<i<N J 

Recall they we are assuming all the measures to be probability measures, thus 
liH{U{l-iA)) = 1- By introducing a Lagrange multiplier for the delta function in 
Eq. (2.9), we get: 

J -00 27r Jxi>o 

(2.10) 

The argument of the partition function can be interpreted as the Boltzmann factor 
of a gas of point charges (Coulomb gas) at positions Aj's on the positive half- 
line (Dyson [1962]). The potential energy of the gas of eigenvalues is given by a 
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harmonic potential and a two dimensional electrostatic repulsion between pairs of 
charges. The analogous of the integral in Eq. (2.10) is known for the case in which 
the integration limits are — oo < Aj < +00, as Selberg's integral (Mehta [2004]). 
The constraint of the positivity of the eigenvalues makes the computation of this 
integral far more complicated. The exact solution can always be found by means of 
the orthogonal polynomials method, but the expressions for Zab grows enormously 
in complexity with increasing N, see (Giraud [2007a]; Facchi, Marzohno, Parisi, 
Pascazio and Scardicchio [2008]) for the first few moments. In the limit of large N, 
we will overcome this problem by looking for the stationary point of the exponent 
of the partition function Zab'- 

l<i<7V l<i<j<N \ l<i<N ) 

The function V represents the potential associated to the Coulomb gas of eigen- 
values. By deriving with respect to both the Aj's and ^ we get + 1 saddle 
point equations for the system: 



-2^7V«A, + 2 E ^^-^^'^ = 0' Vie{l,...,7V}, (2.12) 

E A, = 1. (2.13) 

l<i<Af 

In particular, the domain of integration for ^ lies on the real axis, but we will 
see that the saddle point for ^ lies on the imaginary axis. Thus the contour 
needs to be deformed to pass through this point along the line of steepest descent, 
see appendix A. The minimum or the maximum of V, for positive or negative 
temperatures respectively, in the thermodynamic limit, will correspond to the free 
energy of the system. In the following chapters we will separately analyze the 
range of positive and negative temperatures, and unveil the presence of some critical 
temperatures for the system, associated to first and second order phase transitions. 



2.3 Equivalence between canonical and microcanon- 
ical ensemble 

Before considering different ranges of temperatures, we will devote this section to 
the proof of the equivalence, for our system, between the microcanonical and the 
canonical ensembles. In the first case we represent the total Hilbcrt space % as 
given by isopurity manifolds, that is with a fixed value of the purity. On the other 
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hand in the canonical ensemble we fix only the average purity for each submanifold. 
What follows will be independent on the specific value of the temperature 1//3. 

In order to compare the two approaches, we compute the relative fiuctuations 
of the rescaled average purity: 



where the average purity (vr^^) = N°'{nAB) (j^'ab = A^"vrAB)and its mean square 
fluctuation cx^, in the canonical ensemble are given by (Huang [1987]): 



'AS 



^"^-^ = ^B = ^'-''^ 



= {^ab}-{^ab} = (2.16) 



'AB 



We recall that the coefficient a has to be chosen in order to yield an extensive 
quantity at the exponent of the partition function, namely of (9(iV^) (see Eq. 
(2.8)). It can be easily shown that Jvab = -^n'^g- Fig. 2.2: the indices i and 
j, labeling the purity and its average, identify different isopurity manifolds in Ti. 




Figure 2.2: Microcanonical (a) and canonical (b) ensembles. In the microcanonical 
ensemble the Hibert space "H of a bipartite system is represented in terms of iso- 
purity manifolds, that is manifolds with a fixed value of the purity habi, the index i 
identifying the corresponding submanifold. With the canonical ensemble we release 
the assumption of fixing the exact value of the purity on each submanifold, and fix 
only its average value {'KAB)i by taking into account its fiuctuations (7^^^^. 
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The proof of the equivalence of the two ensembles is thus straightforward: 

- ^^'a. - — — " ^ U J ' ^ ^ ^ 

namely the fluctuations of the average purity can be neglected for a large quantum 
system, N » 1. 



Chapter 3 

Positive temperatures 



In this chapter wc will study how bipartite entanglement between A and B can 
affect the probability distribution of eigenvalues of one of the two subsystems, in 
the range of positive temperatures. In particular, we will show the emergence of 
two classes of symmetries. In Sec. 3.1, by the introduction of a natural scaling for 
the eigenvalues of subsystem A in the limit of a large quantum system, we will 
translate the problem of finding the most probable set of eigenvalues Aj, into the 
task of getting the corresponding density distribution. In Sec. 3.2 the saddle point 
equations will be then reduced to an integral equation for this density function, 
know as Tricomi equation. In Sec. 3.3 we will find a two parameter continuous 
family of solutions, among which we will select the one yielding the higher Boltz- 
mann factor for the Coulomb gas of eigenvalues, i.e. the most probable density 
distribution for every fixed value of the temperature. We will see that there exist 
two classes of symmetries for the density distribution, lying in nearby intervals of 
positive temperatures, that can be interpreted as two phases for the system: the 
maximally entangled and the typical phases. The two regimes are separated by a 
second order phase transition for the system, whose presence will be unveiled in 
Sec. 3.4 by studying the entropy of the system. We will conclude with a reinterpre- 
tation, in Sec. 3.5, of the results obtained in this chapter, directly from the point 
of view of entanglement. 

3.1 Natural scaling and density of eigenvalues 

In Sec. 2.2 wc argued that, in order to yield the correct thermodynamic limit, the 
exponent of the partition function 




(3.1) 
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has to be an extensive quantity, that is of the same order A^^ of the volume of 
states Pa- See Eq. (2.8). As far as we increase (3 typical states enhance their 
degree of entanglement and become maximally entangled in the limit (3 — ?► +oo. 
It then follows that for this range of temperatures the correct scaling exponent is 
a = 3. In order to determine the thermodynamic properties of the system, we solve 
the saddle point equations (2. 12)- (2. 13) in the continuous limit by introducing the 
natural scaling 

= ^A(t,), < = ^ < 1 with Ati = -ti = J^- (3-2) 

Observe that AjAj = A(tj)Atj, where Ai = {i + 1) —i = 1. In Fig. 3.1 we show the 
effects of this scaling from the point of view of the eigenvalues and their labeling 
indices: they become continuous functions in the limit N ^ oo. 

A,. = 1 



J I L 



1 2 i JV 1/N f\ 

AU = l/N 







1/^ 1 1 2 " TV +00 x{t,) 

(a) (b) 



Figure 3.1: Natural scaling, (a) Contraction of indices i from {1,...A^} to the 
interval [1/A^, . . . , 1] associated to ti. (b) Expansion of the eigenvalues Aj from the 
interval [0, 1] to the interval [0, A^] for X{ti). 

In this scaled setup, the properties of the system are determined by the density 
distribution p(A) of the Schmidt coefficients. Obviously, it will depend on (3, that 
is on the level of entanglement between A and B. In the limit A^ — )■ oo we define 
in general: 

p(A) = / dtSiX-Xit)). (3.3) 
Jo 

The probability distribution of the Schmidt coefficients needs to satisfy two con- 
sistency conditions: 

• nonnegativity: p(A) > 0, for all A 



normalization: 



dXp{X) = 1. 



(3.4) 
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As will be discussed in the next section, the above constraints on p(A) will select 
the physical solutions of the saddle point equations, for different values of /3. Notice 
that the standard definition (3.3) for p(A) naturally satisfies the above constraints. 



3.2 Tricomi's equation 



The set of + 1 saddle point equations (2.12)-(2.13) for the Schmidt coefficients 
Xi and the Lagrange multiplier ^, reduces to a couple of integral equations for the 
density p(A) and ^. In both cases, they depend on the inverse temperature (3. 

Let us consider the sum over index j in Eq. (2.12). According to the natural 
scaling it becomes: 



E 



NAt, 



mi 



JV 



N^P 

N^P 
N^P 



dy 



X{x) - X{y) 
dy 

X{x) - X{y) 



dX'5{X' - X{y)) 



dX' 



P(A') 



X{x) - X' 



(3.5) 



where P indicates the Chauchy principal value: 



dX 



> P(A') 



A 



= lim 
A' £->o+ 



A+£ 



A- A' 



(3.6) 



In Eq. (3.5) we have also introduced the identity dX'5{X' — A) = 1. With this 
technique, we get that the saddle point equations, with o; = 3 and /3 > 0, become: 



/3A + P 



dX 



, P(A') 



A -A' 



^2 ^ °' 



^0 



dX p{X)X 



(3.7) 
(3.8) 



Eq. (3.7) is an integral equation for the density of eigenvalues p(A). It is a singular 
Fredholm equation of the first kind, known as Tricomi's equation (Tricomi [1957]). 
According to Tricomi (Tricomi [1957]) the solution, p(A) lies in a compact interval 
[a,b], {0 < a < b). We expect that as this integral equation depends on the 
temperature, also its solution will depend on (3. Let us set 



X — 171 + xS, 4>{x) = p{X)S, 



(3.9) 
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where 



a + b 



b — a 



m — 



2 , ^ = 0<<5<m. (3.10) 

The last expression maps the domain [a,b] of p(A) into the interval [—1, 1] associ- 
ated to 4>{x). Observe that by definition 0(a;) fulfills the normalization condition 
(3.4) 



dx (j){x) 



1. 



(3.11) 



In the new rescaling Eq. (3.7) becomes 

1 



TT 



-1 



y-x 



with 



9{x) 



1 

TT 



i^- + 136m + 135'^x 1 . 



(3.12) 



(3.13) 



We recall that a Lagrange multiplier introduced in the partition function due 
to the unit trace condition on the reduced density operator p^, see Eq. (2.10). 
According to Tricomi's theorem, the normalized solution of this integral equation 
is 

9{y) 1 



The constraint (3.8) 



x^ y — X 



+ 



TT 



X^ 



J A 0(a;) dx — 1 



(3.14) 



(3.15) 



fixes the Lagrange multiplier ^ 



(52 



+ m 2/3 - 



(3.16) 



Notice that ^ is an imaginary number, as anticipated at the end of Sec. 2.2.1. Thus 
the contour of integration for ^ in the partition function needs to be deformed with 

respect to real axes, in order to pass through this point along the line of steepest 
descent. Finally, by substituting the last expression in Eq. (3.13) and performing 
the principal value integral in (3.14) we get 



(f){x) 



' PS"" 2(1 -m) 2 
1 + ^ + — -X - (36^x^ 

Z 



(3.17) 
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The solution of the Tricomi equation depends on the fictitious (inverse) tempera- 
ture p. The physical solutions must have a density (plx) that is nonnegative for all 
X G [—1, 1]. The numerator of (3.17) is non negative for x G [a;_, with 



1 / 1 — m 



x± 



±v^), (3.18) 



where 

(3,., 

Notice that A > for every m and 6. The condition [—1, 1] C imposes 
that the points {6, m) should be restricted to a (possibly cut) eye-shape domain 
given by: 

max{5,r^((5,/3)} < m < r+(5,/3), (3.20) 

where (3.10) expresses the positivity of eigenvalues and m = Tf[5,l3) are the level 
curves corresponding to x± = ±1: 

r^(5,/3) = l±^(l-^). (3.21) 

See Fig. 3.2(a). 




Figure 3.2: (a) Solution domain for /3 = 0.5. The domain is the region with 
full color, whereas the curved and straight dashed lines are Tf{S,l3) and m = 5, 
respectively. The line m = 6 corresponds to the positive eigenvalues condition, (b) 
Solution domains for different temperatures: for each value of /3 (indicated) the 
relative full line encloses the region of the parameter space where the eigenvalue 
density is positive. For each temperature, the red dot indicates the right corner of 
the eye-shaped domain. 
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Observe that the level curves TfiS, 13) are symmetric with respect to the hne m — 1 
and intersect at 5 = and at 5 = •\/2//3. The right corner of the eye is at 



{5, m) = 




(3.22) 



and belongs to the boundary as long as (3 >2. For ^ < 2 the eye is cut by the hne 
m — S, as shown in Fig. 3.2(b). In other words the domain is symmetric as long as 
P >2, whereas for lower values of (3 this symmetry is lost. Let us remark that all 
points inside each domain correspond to solutions of the saddle point equations, 
that is, we have a two parameter continuous family of solutions. 

In the next section we will show how to select the density of eigenvalues which 
maximizes the partition function. It corresponds to the most probable distribution 
of the Schmidt coefficients, for a given temperature. In particular, we will distin- 
guish two classes of distributions. They will correspond to different ranges of /3, 
associated to the two symmetries we have found for the physical domains of the 
solution of the Tricomi equation: < /3 < 2 and /3 >2. 



3.3 The solution 

We will now compute the density of eigenvalues which gives the maximum contri- 
bution to the partition function. It is the solution corresponding to the maximum 
Boltzmann factor for the gas of eigenvalues, or equivalently to the minimum of PV 
(or V since /3 > 0). The potential computed for this minimizing solution will be 
the free energy of the system. 



3.3.1 Thermodynamics of the system and minimization of 
the free energy 

By applying the natural scaling (2.8), the potential for the gas of eigenvalues (2.11) 
with a — 3 becomes 

l<i<N ^ l<i<N i<j<N 

E E in'v 

^ l<i<N i<j<N 

1 1 , .^/IniV 
= u- -s+ - \tlN + 



f3 (3 \ N 

= / + ^lniV + o(^). (3.23) 
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where 



l<i<N 

^ = ]^ E ^^\m-m\ (3-25) 

l<i<j<N 

f = lim (f^-hnN) . (3.26) 



In the hmit — >■ oo we have: 

-1 



u = / da;A>(a;), (3.27) 



dy4>{x)<i>{y)ln{S\x-y\). (3.28) 



Observe that /, u and s satisfy the thermodynamic relation 

Pf = /3u- s, (3.29) 

and they will be the free energy, the internal energy and the entropy densities of the 
Coulomb gas of eigenvalues when computed on their most probable distribution. 
In order to compute s we integrate Eq. (3.12) 

r dx P C dy = TT r dz giz) (3.30) 

J-i 7-1 y-x j_i 

from which it follows 

j dx (j){x) J dy (l){y)]n\y — x\ — J dx (p{x) hi{x + 1) 

/I PX 
dx(f){x) J dy g{y). (3.31) 

Thus by inserting in Eqs. (3.27)-(3.29) the physical solution (3.17) of the Tricomi 
equation and observing that 



\y '"'^^-h2, f iy "■"'j^' = 1. (3.32) 

1 Try 1 — J-i Try 1 — 7/2 

dy == = ----ln2 (3.33) 

1 7ryl-y2 4 I 
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we get 



u{S, m, (3) 
m, (3) 



2(1 - 



m] 



m] 



52 



52 (55^ 
' 2" ~ ~8~' 

^- hln-. 

16 2' 

2n -m)^ 
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(5^ 
2 



16 



In 



(3.34) 
(3.35) 
(3.36) 



In Sec. 3.2 we have determined the domain of p(A) in the space (m,(5), for dif- 
ferent ranges of the inverse temperature /3. As aheady observed all the points of 
this domain yield a solution of the saddle point equations (2.12)-(2.13). In order 
to find the most probable distribution of the Schmidt coefficients we now have to 
minimize /?/ (or /) in the physical domain (3.20). The contour plots of the free 
energy density are shown in Fig. 3.3. 




i = i 




(b) 



1 .0 1 .5 
5 



m 1.0 



/3 = 3 



0.0^ 
0.0 



1.0 
5 



Figure 3.3: Contour plots of /3/ in regions of the parameter space such that 0(a;) > 
0, for (a) /3 = 0, (b) /3 = 1 and (c) /3 = 3. Darker regions have lower free energy. 



Notice that /, as well as u and s, is symmetric with respect to the line m = 1. 
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This Z2 symmetry will play a major role in the following. The only stationary 
point (a saddle point) of (3f is at the right corner of the eye (3.22), see Figs. 3.2 
and 3.3. Thus, the absolute minimum is on the boundary. The behavior of the free 
energy at the boundaries of the allowed domain is shown in Fig. 3.4 for different 
temperatures. 




Figure 3.4: Free energy density 13 f on the boundary of the region of the domain 
where 4>{x) > 0, for different temperatures (indicated). Dashed lines: (3f on the 
lower boundary of the eye-shaped domain. Full lines: Pf on the upper boundary. 
The sought minima of the free energy can be inferred from the graphs and coincide 
with the dots in Figs. 3.2 and 3.3. 

We deal with two regimes of solutions corresponding to different ranges of tem- 
perature. Recall that the temperature in our model is the Lagrange multiplier 
associated to isopurity manifolds of the Hilbert space. Different symmetries of the 
physical domain will correspond to different symmetries in the distribution of the 
Schmidt coefficients. More precisely we have: 



for /3 > where /3+ = 2, when the domain of the physical solution of the 
Tricomi equation is eye-shaped, the absolute minimum is at the corner of the 
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eye, Eq. (3.22) 



171=1, S = 




(3.37) 



• for < /3 < the absolute minimum is at the right upper corner of the 
cut eye-shaped allowed region, 5 = r^{6,p), namely at 

m = 5, /3=^-|. (3.38) 

See the red dots in Figs. 3.2(b) and 3.3. For the distribution of eigenvalues corre- 
sponding to these minima, / coincides with the free energy of the system, while u 
and s can be identified with the internal energy density and the entropy density, 
respectively. 



3.3.2 Distribution of the Schmidt coefficients 

We will study the behavior of our system starting from high values of /3, that is 
low values of internal energy u (purity). In the previous sections we unveiled the 
presence of two regimes. In this section we will show that they correspond to a 
different symmetry in the distribution of the eigenvalues. 



Wigner semicircle distribution: /3 > - towairds mctximally entangled 
states 

For /3 > /3+ = 2, we have found that the minimum of the free energy for the gas of 
eigenvalues corresponds to the corner of the eye (Eq. (3.37)) From the expression 
(3.17) for the physical solution of the saddle equations, we get the semicircle law 
(see Fig. 3.5) 

2 / 

(j){x) = -Vl -^2, (3.39) 

TT 

whence, by (3.9), 

p(A) = -^/X^Vb^, (3.40) 

TT 

where 

a=l-5=l-J%, 6=1 + 5=1+,//^. (3.41) 



This distribution is displayed in Fig. 3.5. Observe that as ^ becomes larger the 
distribution of the Schmidt coefficients becomes increasingly peaked around 1. 
This means that all the eigenvalues tend to 1/A^. See the natural scaling (3.2): 
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(a) 1 



-0.5 0.5 1 

X 




Figure 3.5: (a) Density of the eigenvalues for (3 > (3+ = 
eigenvalues for /? = 2, 4 and 10. In the temperature range 
is given by the semicircle law. 



2. (b) Density of the 
' e [^+, oo[ the solution 



for temperatures T = tending to zero quantum states become maximally 
entangled. The internal energy and the entropy of the system for this solution can 
be easily obtained by plugging the condition on the parameters m and 5 (3.37) 
into the expressions (3.34)-(3.35) for u and s respectively: 

- ' (3.42) 



u 



s — 



1 + 



1 + 



1 , ^ 

h In - = 

4 2 



1 
1 

~4 



ln(2/3), 



and thus 



ln^ = /3+^ + ^ln(2/3). 



(3.43) 



(3.44) 



Wishart distribution: < /3 < /3+ - typical states 

At higher temperatures < /3 < /3+ the solution acquires a different shape, be- 
coming asymmetric with respect to the center of its support. Analogously to the 
previous case, by plugging the relations (3.38) between m, 5 and P at the minimum 
of /3f into the physical solution of the saddle point equations (3.17), we get the 
distribution of the Schmidt coefficients: 

By (3.9) it yields 
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with b — 25. See Fig. 3.6. This is the Wishart distribution. 




X A 

Figure 3.6: (a) Density of the eigenvalues for (3 = 1. (b) Density of eigenvalues 
for ^ = 0, ^ = 2/3, and p = /3g = —2/27 (dashed). In the range of temperatures 
/3 e {/3g,/3+), with /3+ = 2, the solution is given by the Wishart distribution. We 
refer to the next chapter for the range /3g < /3 < 0. 

The change from the semicircle to the Wishart distribution is accompanied by a 
phase transition (the first of a series), as we shall see in Sec. 3.4. In Fig. 3.7 we 
plot the relation (3.38) between the half width of the distribution 6 = b/2 and /3. 
It runs monotonically from (3 = 2 when S = 1 to (3 = when 5 = 2. 




Figure 3.7: Plot of Eq. (3.38) for positive (sohd hne) and negative (dashed line) 
temperatures. The minimum (3g = —2/27 is attained at S = 3. 

Moreover, the inverse temperature (3 reaches a minimum equal to 



(3.47) 
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at 5 — 3. Therefore, the above solution can be analytically continued down to 
/3g, which is shghtly negative, but not below. We will study the solution for neg- 
ative temperatures in the next chapter. Let us now compute the thermodynamic 
quantities for the Wishart distribution. The internal energy u (proportional to the 
average purity) is obtained by plugging (3.38) into the expression (3.34) 

3 P 

u = ^S-j. (3.48) 

In particular, at ^ = (5 = 2) we get u — 2, at /3+ — 2 {6 — 1) we get u — 5/4, 
and at pg = -2/27 {5 = 3) we get u = 9/4. See Fig. 3.8. Observe that u{(3) is 
continuous at From Eqs. (3.38) and (3.35) we can also compute the entropy 

9/4 
2 

u 

5/4 



PgO 2 4 

Figure 3.8: -u as a function of the (inverse) temperature. Notice that m = 2 at 
(3 = (typical states). In the /9 — )■ oo limit we find the minimum u = 1. We 
will see that (3g = -2/27, u = 9/4 (left point) and (3+ = 2,u = 5/4 (right point) 
are critical points for the system. However the phase for < ^ < is unstable 
(dashed line) towards another phase as soon as ^ < (in this scahng of /3) . 

and the free energy for the Wishart distribution 

9 5 3,5 

' = -4 + ^-5^ + ^"2' 
9 9 11 5 

= ^-^ + T-^"2' 

in terms oi d = S{/3) e]l,3]. 

Summarizing, for (3 — >■ +00 all eigenvalues of pa are equal to 1/A^, that is sub- 
system ^4 is in a maximally mixed state. As the inverse temperature is decreased, 
for /3 > all eigenvalues remain of O (1/iV) their distribution being character- 
ized by the Wigner semicircle law, while for < /3 < /3+ the eigenvalues, still of 




(3.49) 
(3.50) 
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O (l/N), follow the Wishart distribution which diverges at the origin, that is some 
eigenvalues of pA vanish. Let us remark that the distributions of eigenvalues we 
have found are consistent with the scaling coefficient a = 3. Indeed, for any value 
of the inverse temperature /3, the domain of p(A) involves the rescaled eigenvalues 
A which are of that is the eigenvalues Aj are of 0{1/N). In other words, 

the range of positive temperatures refers to those submanifolds of the total Hilbert 
space T-L given by pure states showing a high level of entanglement between 
subsystems A and B. Observe that the internal energy is of that is the purity 
is of 0{1/N), as will be discussed at the end of this chapter. 

3.3.3 Analysis of moments and cumulants 

Notice that from the relation 

2AB(/3) = e-« (3.51) 

we have that the free energy is proportional to the generating function of the 
connected correlations of tt^b- More precisely, the derivatives of Zab{P) evaluated 
for /3 = 0, yield the moments oii^AB with respect to the measure dii{pA) introduced 
in Sec. 2.1.1, i.e. 



'-1Y d^ZABil^) 



(3.52) 

/3=0 



These functions fully determine the statistical distribution of tt^b on the set of 
pure bipartite states and in the high temperature regime provide an expansion of 
Zab{I^)- In (Facchi, Marzolino, Parisi, Pascazio and Scardicchio [2008]) we can 
find the expression for the n— th cumulant /C„ of the purity tiab (in one to one 
correspondence with the first n moments) in the limit of large N and ^ — )■ 0, for 
the case of balanced bipartitions, dim I-La — dim I-Lb — N: 



(-l)na"(/3/) 



7V3n Qfjn 



<^"-^)' (3.53) 



/3-)-0 ^ ^ 



and also the cumulants for the more general case of unbalanced bipartitions. For 
instance, the first and second cumulants are given by 
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with N = dim Ha and M = dimUB- They reduce to (3.53) for = M and 
n — 1,2. The exact expression (that is vahd also for small A?") at ^ = of the 
first 5 cumulants can be found in (Giraud [2007a]). In particular for the case of 
balanced bipartitions, dim?^^ — dim?^B = A^, the first two cumulants are given 

by 

= if^ (3-56) 

r = 2(iv^ - 1)^ 

' (l + Ar2)2(2 + Ar2)(3 + iV2)' ^ ■ ) 

while for more general case of dim Ha 7^ dim T-Lb 'we have 

2(Ar2-i)(M2-l) 



(1 + NMf{2 + ArM)(3 + ATM) ' 

(3.59) 

We recall the relationship between the cumulants and moments of a generic prob- 
ability distribution 

^n^Mn- Yl (^Ii)^mMn-m, (3.60) 

l<m<n-l ^ ^ 

and in particular 

/Ci = Ml (3.61) 
/C2 = A^2--M?- (3.62) 
More generally, we can also define the moments of tvab ^or /3 7^ as 

where diJ,i3{ip) is the canonical measure 

dl^nipA) = dp{pA) — 7^ . (3.64) 

The latter is a deformation of the Haar measure dfi{pA) obtained by including 
a non uniform weight which explicitly depends upon the purity ttab, thru (3. In 
particular, as /3 increases, dpp{pA) enhances the role of the states with lower values 
of ttab (i-e. larger values of bipartite entanglement between the two subsystems 
A and B), to the extent that for ^ — >■ +00 only the maximally entangled states 
contribute to (3.63). 
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3.4 Second order phase transitions 

We are now ready to unveil the presence of the first critical point for the system 
at /3+ = 2. In the previous section, we have found that the densities of eigenvalues 
corresponding to the two contiguous intervals of temperature < (3 < (3+ and 
(3 > (3+ show a different symmetry, which correspond to the Wishart and the 
Wigner semicircle distributions. They are associated to two distinct entanglement 
phases for the system: the maximally entangled and the typical phases. In the 
former all the eigenvalues are of 0{1/N), whereas in the second regime some 
eigenvalues go to zero. A second order phase transition at /3+ separates the two 
phases. It is due to the restoration of the Z2 symmetry P ("parity"), present in 
the thermodynamic quantities, see Eqs. (3.34), (3.35) and (3.36). It corresponds 
to the reflection symmetry of the semicircle distribution p(A) around the center of 
its support {m = 6 = b/2 ioT 13 < (3+ and m = 1 for /3 > /3+), that is not present 
in the Wishart distribution. 

The direct proof of the presence of a second order phase transition in the system 
at p — f3^ can be obtained from the expression of the entropy density s — f3{u — f). 
See Eq. (3.43) for l3 > l3+ and Eq. (3.49) for < 13+ . Prom both equations we get 
the same values of the entropy at /3 = /3+ (5=1) 

s = -1/4 -In 2. (3.65) 

On the other hand the first derivative of s with respect to 5 is discontinuous at 
5 — 1. However, also ^ as a function of 5 has a discontinuous first derivative at 
5 — 1, see Eqs. (3.37) and (3.38). By reminding that 

ds ds /d(3 
dp ^ ds/dS' 

d?s d^s / fdpV dsd^p/fd[3V 

(3.66) 



dd'^l \dS J dS dS^/ \dS 
one easily obtains that the discontinuities compensate. For /3 — )■ /5+ we have 

. ^ -i - ln2 - ^ + 0(13 - ^+)^^^, (3.67) 

where 6 is the step function. Thus the system undergoes a second order phase 
transition at (3+ that, however, is fairly mild. Indeed, the entropy s, together with 
its first derivative, is continuous, although the second derivative shows a finite 
discontinuity of 1/8 at ^ = /3+. See Fig. 3.9. In other words, the heat capacity of 
the system, corresponding to ds/dj3, does not diverge, although its first derivative 
is discontinuous. 
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Figure 3.9: Entropy (a) and its second derivative with respect to /3 (b). The entropy 
is continuous in /?_)_ while its second derivative presents a finite discontinuity. 



Observe, also, that the entropy is unbounded from below when /? — > +00, see 
Eq. (3.43). The interpretation of this result is quite straightforward: the minimum 
value of {71 ab) is reached on a submanifold (isomorphic to SU{N)/7jn (Sinolecka, 
Zyczkowski and Kus [2001])) of dimension N"^ — 1, as opposed to the typical vectors 
which form a manifold of dimension 2N'^ — — 1 in the Hilbert space Ti. Since 
this manifold has zero volume in the original Hilbert space, the entropy, being the 
logarithm of this volume, diverges. 



3.5 Entropy density vs average entanglement 

We conclude the analysis of positive temperatures, reinterpretating our results from 
the point of view of quantum information. To this end we will write the entropy 
density s as a function of the internal energy density u, or average entanglement 
{t^ab)- In Sec. 3.3.2 we computed these thermodynamic quantities as functions of 
the parameter 5, both for the semicircle law and for the Wishart distribution. If 
we invert the expressions for the internal energy, Eqs. (3.42) and (3.48), we get 

r 2v^7^, i<u<\, 

S = I (3.68) 
[3- V9-4m, l<u<2 
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which plugged into Eqs. (3.43) and (3.49) give 



s{u) 



4' 



1 < U < 



4' 



Inff 



(3.69) 



2 ; 4 



<u<2. 



This function is plotted in Fig. 3.10. It counts the number of vectors with a 
given (average) degree of bipartite entanglement between subsystems A and B: 
s = {lnV)/N'^ where V is the volume occupied by each submanifold. 




Figure 3.10: Entropy density s versus internal energy density u 
Eq. (3.69). 



N{ttab)- See 



In particular by expanding (3.69) in the critical region, that is for u ^ 5/4, we get 



4 U- 



+ - 



32 



-6 [ u 

2 V 4 



u 



+ 



u — 




(3.70) 



Therefore, the second order phase transition at u = 5/4 is signaled by a dis- 
continuity in the third derivative of s{u), see Fig. 3.11, whereas at this point 
the entropy shows a discontinuity of the second derivative with respect to /3, see 
Fig. 3.9. In fact, in general, if s,m and T are entropy, energy and temperature, 
respectively, and C = du/dT is the specific heat, one gets ds/du = (3 = \/T and 
d'^s/du^ = -1/{T'^C) = -{l/T^){ds/dT)-^. Discontinuities of the n-th derivative 
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Figure 3.11: Entropy (a) and its third derivative with respect to u (b) in the critical 
region u — )■ 5/4. The entropy is continuous at m = 5/4 while its third derivative 
displays a finite discontinuity. 

of s{T) translate, therefore, into discontinuities of the (n+ l)-th derivative of s{u). 
Let us finally discuss the significance of these results. 

In the range of positive temperatures T — 1/p the eigenvalues of the reduced den- 
sity matrix of our dimensional system are always of O {1/N). As a consequence, 
the value of energy (average purity) in Eq. (2.1) 



is always small: there is therefore a lot of entanglement in our system. Notice that 
the purity becomes an intensive quantity if multiplied by iV: u = N((Kab)- There 
are however, important differences as purity changes (it is important to keep in 
mind that in the statistical mechanical approach pursued here, the Lagrange mul- 
tiplier 13 fixes the value of the energy /purity). As pointed out in Sec. 3.3.2, when 
1/N < (tiab) < 5/4A^ (where /3 > /3+) the eigenvalues are distributed according to 
the semicircle law (Fig. 3.5), while for 5/4A^ < (ttab) < 2/iV (where < /3 < /3_) 
they follow the Wishart distribution (Fig. 3.6). The two regimes are separated by 
a second order phase transition. Notice that the value (ttab) — 2/N, which cor- 
responds to infinite temperature /9 = 0, refers to vectors in the Hilbert space that 
are typical according to the Haar measure. One is therefore tempted to extend 
these results to negative temperatures (Facchi, Marzolino, Parisi, Pascazio and 
Scardicchio [2008]) and indeed this can be done up to {hab) = 9/4A^, correspond- 
ing to the slightly negative temperature — —2/27. However, as we have seen, a 
mathematical difficulty emerges and no smooth continuation of this solution seems 
possible beyond 



( i:A?)4/wA^o(i) 



(3.71) 
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In the next chapter we will see that two branches exist for negative /3: one 
containing a phase transition ai f3 — f3g and in which purity is always of O (l/N) 
and one in which purity is of O (1). The latter becomes stable as soon /? becomes 
negative through a first order phase transition. We remind that larger values of 
purity, towards the regime ttab = (1) yield separable (factorized) states. We are 
therefore going to look at the behavior of our quantum system towards separability 
(regime of small entanglement). 



Chapter 4 



Negative temperatures: 
metastable branch 



By analytic continuation, the solution at positive temperatures of the previous 
chapter can be turned into a sohition for negative /3, satisfying the constraints of 
positivity and normalization. As this is an analytic continuation of the solution 
obtained for ^ > 0, we are not assured that this is indeed a stable branch. In 
this chaper we will study this analytic continuation, but we anticipate that this is 
metastable as soon as /3 becomes negative (namely for /3 < — 2.455/iV). We call 
this solution the metastable branch, postponing the proof of its metastability to 
Chap. 5. The metastable branch is of interest in itself because, as we shall see, it 
entails two phase transitions. The first one is at /? = — 2/27 = l3g and corresponds 
to the so-called 2D quantum gravity free energy (see (Di Francesco, Ginsparg and 
Zinn- Justin [1995])), provided an appropriate double-scaling limit (jointly /3 — )■ /3g 
and N — )■ oo) is performed. The second phase transition, at /3 = — /3+ = —2, is 
due to the restoration of the Z2 symmetry that was broken at the phase transition 
at ^ = ^+ = 2, described in the previous chapter. 

This chapter is organized as follows. In Sees. 4.1 and 4.2 we will determine the 
physical solution of the Tricomi equation for the range of negative temperatures 
and select the one which maximizes the potential energy density of the system. 
Thus, in Sec. 4.3, we will unveil the presence of the two mentioned second order 
phase transitions for the system, and discuss them by studying the behavior of 
the entropy density. Finally, we will devote Sec. 4.4 to the overview and the 
reinterpretation of our results in terms of the average entanglement. 
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4.1 Physical domain and maximization of the free 
energy 

By keeping the same value of the scahng coefficient a — 3 used for the positive 
temperatures, the saddle point equations for the potential energy of the Coulomb 
gas, lead in the continuous limit to the Tricomi equation (3.12)-(3.13), now with 
/3 < 0. We recall that its solution is given by Eq. (3.17) that is 



If. 2(1 -m) 2' 



(4.1; 



being defined for x e [—1,1]. The density distribution is nonnegative for x ^ 
(x_,a;+), where (t>{x±) — and 

with 

-"'^^,,.^+^V (,3) 



See Eqs. (3.18) and (3.19). For /3 < we have A > when 

r2-(<5,/3)<m<r+((5,/3) (4.4) 
where m = Pf (5, (5) are the level curves associated to A = 



r±(5,/3) = l±5y-/3(^l + ^). (4.5) 

As displayed in Fig. 4.1(a) at /3 = —5, they are symmetric with respect to the line 
m — 1 and intersect at 5 = and 5 — •>/— 2/^. On the other hand the condition 
(— 1, 1) n [x-, x+] — 0, implies 

rr(^,/3)<m<r+(5,/3) (4.6) 

being m = ^"^{5, (5) the level curves (3.21), x± = ±1, already computed for /3 > 0. 
For /3 < they are still symmetric with respect to the line m = 1 but intersect 
only in 5 = 0. The two couples of level curves, Vf and intersect at points 



(..™)H,/-|.l±^f|^). (4.7) 
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Figure 4.1: Metastable branch. Domain of existence of the solution (m,(5) for 
negative temperatures, (a) Solution domain at /3 = —5. (b) Solution domains at 
different temperatures: for each value of /3 (indicated) the relative full line encloses 
the region of the parameter space such that the eigenvalue density is positive, and 
the red dot indicates its right corner. 



Therefore, conditions (4.6) and m < (5, expressing the positivity of the eigenvalues, 
determine the physical domain of the solution 0(x) of the Tricomi equation, namely 
the (possibly cut) eye-shaped domain 

max {(5, /i_(/3, 5)} < m < /i+(/3, 5), (4.8) 

where 



rf(5,/3), 0<(5<y/-^, 
/i±(5,/3) = <' (4.9) 



The right corner of the eye is given by 



2^ 



(5,m)= (4.10) 

and belongs to the boundary as long as /3 < — /3+ = —2. This equation is the 
analogue of Eq. (3.22) for positive temperatures. For /3 > —4 the eye is cut 
by the line m = 6. See Fig. 4.1(b). Notice that, similarly to what found for 
the case of positive temperatures, the physical domain does not display the same 
symmetries for the entire range of negative temperatures. We have already noticed 
that all points inside the physical domain correspond to solutions of the saddle 
point equations, for every fixed temperature. The solution corresponding to the 
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most probable distribution of the Schmidt coefficients can be found by maximizing 
the free energy density / (or — /?/, /3 < 0) of our Coulomb gas. For the explicit 
expression of the free energy we still refer to Eq. (3.36). The contour plots of the 
free energy density are shown in Fig. 4.2. 




1 2 3 

S 



Figure 4.2: Metastable branch. Contour plots of — / in regions of the parameter 
space such that 0(a;) > 0, at (a) (3 = f3g, (h) f3 = —1 and (c) (3 = —5. Darker 
regions have lower free energy. 

The behavior of Pf at the boundaries of the allowed domain is shown in Fig. 4.3 for 
different temperatures. Notice that the free energy density displays a Z2 symmetry 
with respect to the line m = 1, as already found for the range of positive /3. See 
Fig. 3.3. The symmetries displayed by the physical domains for different ranges of 
/3 will naturally emerge in the distribution of the Schmidt coefficients. Apart from 
/3 = Pg, the free energy density / has no stationary points for /3 < 0. Thus, the 
most probable distribution will correspond to an absolute minimum of /?/ on the 
boundary. 

In particular we have that: 

• for P < — /5+ = —2 the right corner of the eye (4.10) is the global minimum, 
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Figure 4.3: Metastable branch. Free energy on the boundary of the region of the 
domain where 0(a;) > 0, for different temperatures (indicated). Dashed hues: free 
energy /3/ on the lower boundary of the eye-shaped domain; full lines: free energy 
on the upper boundary. The sought minima of the free energy can be inferred from 
the graphs and coincide with the dots in Figs. 4.1 and 4.2. 

• for —(3+ < (3 < the absolute minimum of — / is at the right upper corner 
of the allowed region, namely at 

m = 6, 6 = h+{P,6). (4.11) 

The division of the range of negative temperatures in two parts is the first signature 
that the point /3 = — will correspond to a change of symmetry in the distribution 
of the Schmidt coefficients, as can be expected from what seen for /3 = in the 
previous chapter. Condition (4.11) implies a further arrangement in subintervals 
of the range — /3+ < (3 < 0. Recall that the boundary h+{6, f3) consists of two 
level curves, r^(5, /3) and r^{6,P), corresponding to nearby intervals of inverse 
temperature. See Eq. (4.9). Therefore Eq. (4.11) entails 

• ioT (3g</3 <0 

m = 6, P = ^,~l (4.12) 
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for -l3+ < l3 < l3g 



m 



6, 6-1 



/3 1 + 



,5^ 



(4.13) 



Notice that Eq. (4.12) coincides with the condition (3.38) found for < /3 < 
and leads to the Wishart distribution. We conclude that Eq. (4.12) is the analytic 
continuation of the curve (3.38) which runs monotonically from /3 = when 6 = 2 
to its minimum f3g = —2/27 at 5 = 3. See Fig. 3.7: the continuous and the dashed 
lines refer to < /3 < /3+ and f3g < (3 < 0, respectively. On the other hand, 
Eq.(4.13) 



6 



/3 1 + /3 



52 



(4.14) 



admits two real solutions: 



i, + js^{2 + V2-6)i6-2 + V2), 3<6<2 + V2 



/3 



1 i_ 

'<52 S3 



(4.15) 



{2 + V2-6){6-2 + V2), 1 <6 <2 + V2. 



They run from /3 = Pg, when 6 = 3, (with derivative zero) up to (3 = -3/2 + y/2 
when 6 = 2 + a/2 (with derivative —00) and then from j3 = —3/2 + \/2, when 
6 = 2 + a/2, (with derivative +00) up to /3 = — when 6 = 1. See Fig. 4.4. 



/5 




3 2 + V2 



Figure 4.4: Metastable branch. (3 vs 6. See Eqs.(4.12) and (4.15). 
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4.2 Distribution of the Schmidt coefficients and 
thermodynamics of the metastable branch 

In this section we will compute the density function for the Schmidt coefficients. 
They will display different symmetries for different ranges of /3 (or S), reflecting 
the symmetry of their domains discussed in the previous section. We will explore 
the whole range of negative temperatures, from /3 = (5 = 2) to /3 ^ — oo (5 = 0). 
We will then compute the entropy, the internal energy and the free energy of the 
Coulomb gas of eigenvalues. 



Wishart distribution: l3g < ^ < - typical states and quantum gravity 

When (3g < f3 < the probability distribution of the Schmidt coefficients is ob- 
tained by plugging (4.12) into the generic expression (4.1) for the solution of the 
Tricomi equation 



0(^) = ;|y^(l + (2-5)a;), 2<S<3. (4.16) 

It is the Wishart distribution (3.45) we have found for ^ < ^+ in Chap. 4. At /3g, 
where the inverse temperature 

" = ^-1 ("^) 

is minimum, for S = 3, (see Fig. 3.7) the density function becomes 



or equivalently 



It displays a different shape with respect to the other distributions in this range 
of temperatures. Indeed for (3g < f3 < (3^ the derivative at the right edge of the 
distribution diverges, whereas at /3g it vanishes. See Fig. 4.5. In Sec. 4.3 we will 
show that the Coulomb gas of eigenvalues undergoes a second order phase transition 
at /3g. The internal energy, entropy and free energy for the system corresponding 
to this family of solutions can be determined by inserting (4.17) with m = S in 
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Figure 4.5: Metastable branch, (a) Density of the eigenvalues at /3 = —1/27 (solid) 
and (3g (dashed), (b) Density of eigenvalues at (3 = —1/54, the solution is given by 
the Wishart distribution. 



Eqs. (3.34)-(3.36) 

u = (4.20) 
9 5 3,5 

11 9 9 5 

These expressions determine the thermodynamic properties of the Coulomb gas, 
and will be employed in the analysis of the phase transition for the system. 



Asymmetric Eircsine distribution: 
mally entangled states 



-2 = < ^ < /3„ 



towards mctxi- 



For this range of temperatures the relation between j3 and 5 is more involved. 
Indeed we find two famihes of solutions. By (4.15) and (4.1) we have that: when 
— 3/2 + \/2 < /3 < /3g the distribution of the Schmidt coefficients is 



rl 
2 



+2(1 - 6)x +(6- + 45-2 



X 



(4.23) 



with 



+ ^^(2 + V2- S){S -2 + V2), 3<6 = m<2 + V2, (4.24) 
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while for -2 < ^ < -3/2 + ^/2 we have 



+2(1 - S)x +(^S + V-S^ + 45-2) 



X 



with 



1 



1 
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(2 + a/2-5)(5-2 + \/2), 1<5<2 + a/2. 



(4.25) 



(4.26) 



The density functions corresponding the these two nearby subintervals of inverse 
temperatures are shown in Fig. 4.6. Notice that this eigenvalue density diverges 
both at the left edge x — —1 and at the right edge x — +1. It can be seen as an 
asymmetric arcsine distribution, as it is not symmetric with respect to the center 
of its domain, x — 0. 





(b) 



X 



A 



Figure 4.6: Metastable branch, (a) Eigenvalue density for /3 = /3g — 0.01 G [—3/2 + 



\/2 < 



< 



(see Eq. (4.23)). (b) Eigenvalue density for /3 = -1/4 e [-2 < 



/3 < -3/2 + ^/2] (see Eq. (4.25), with X^m + x5, Eq. (3.9). In both cases, the 
distribution of the Schmidt coefficients diverges at the edges and is not symmetric 
with respect to the center of its domain. 

By plugging (4.24) and (4.26), with m — 5 in the Eqs. (3.34)-(3.36) we compute 
the internal energy, entropy and free energy for the system. Beyond Pg we get for 
-3/2 + ^/2<P<Pg{3<6<2 + ^/2) 



u — 


2S - -S' - 
8 




s — 


AS 


15 

" w 




5 25 


17 

8^ " 
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8^ 



2 



+ 45 _ 2 - In 



2' 



(4.27) 
(4.28) 
(4.29) 
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and for -2 < /? < -3/2 + \/2(l<(5<2 + \/2), 

q r 

= 25- -5' + -V-52 + 45-2, (4.30) 
8 8 

. = -2+S-^-^V-.^ + 4.-2 + ln^, (4.31) 

/5/ = + V-5^ + 45-2-ln^. (4.32) 



2 4(5 8(52 ' \^8(5 (5V 2 

Arcsine distribution: fi < — /3+ = —2 - almost maximally entangled states 

At /3 = — /3+ (5 = 1) the Z2 symmetry is restored. We have seen that for /3 < — 
the global minimum of the physical domain is given by Eq. (4.10) 



m = l, <^=y-|- (4.33) 

Observe the analogy with Eq. (3.37) for the Wigner distribution, p > j3+. We 
immediately have that for all inverse temperatures lower than — /5+ the solution of 
the Tricomi equation is given by 

2r2 

<P{x) = - , < 5 < 1 (4.34) 



X 



P(A) = , P < -2, (4.35) 



that is 



7r^l-f(l-A)2 

being A = m -\- 5x (Eq. (3.9)). See Fig. 4.7. The eigenvalues distribution is 
symmetric with respect to the center of its domain. In the next section, we will 
show that at /3 = — ^+ = —2 the system undergoes a second order phase transition 
related to the Z2 symmetry, that mirrors the critical point at /3+ = 2. The function 
(4.34) or (4.35) is the arcsine distribution. Furthermore, while (t){x) does not 
depend on the temperature, p(A) shrinks around A = 1 as ^ lowers, i.e., recalling 
the natural scaling (3.2), Ai = A2 = . . . = Aat = 1/N for ^ — )■ —00. In other words, 
the metastable branch for the system leads to maximally entangled states, for very 
high — /3's. However, as already anticipated at the beginning of this chapter, we 
will show that this solution, characterized by the same scaling coefficient a = 3 
of positive temperatures, belongs a metastable branch for the system. On the 
contrary, the stable branch will approach the region of states showing a significant 
separability between systems A and B. Furthermore, we will find that the volume 
occupied by states with (tt^b) = 1 (/3 ^ —00) diverges. 
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Figure 4.7: Metastable branch, (a) Eigenvalue density for (3 < — /3+ = —2. (b) 
Eigenvalue density for /3 — —2,— 5,— 10. For /3 < —/3+ the Z2 symmetry of the 
eigenvalues distribution with respect to m = 1 is restored. 



The thermodynamic functions u (internal energy density) and s (entropy den- 
sity) computed for the arcsine distribution are 

u = l + ^^' = l-^> (4-36) 
s = ln^-^ = -^ln(-2/3)-i, (4.37) 

The next section will be devoted to the analysis of the two critical points for the 
metastable branch, and — 



4.3 Second order phase transitions for the 
metastable branch 

The metastable branch for /3 < 0, is characterized by the presence of two second 
order phase transitions for the system. Similarly to what seen in the previous 
chapter for the range of positive temperatures, they correspond to a change in the 
symmetry of the distribution of eigenvalues. The thermodynamic quantity we will 
consider in order to characterize these critical points is the density entropy. 

The first phase transition of the system, as /3 lowers from to —00, appears 
at the gravity temperaure, Pg, that is S — 3. This point is at the convergence of 
two solutions, the Wishart and the asymmetric arcsine distribution, given by Eqs. 
(4.16) and (4.23) respectively. We have also observed that the Wishart distribution 
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at (3g shows a vanishing derivative at the right edge. See Figs. 3.6(b) and 4.5(a). 
In order to characterize this phase transition we expand Eqs. (4.17) and (4.24) for 
5 — >■ 3 and invert the series: 



.| + i(2-3»(*-3))(*-3r-^(^-3)' 



- 3)) (6 -3)' + {{S - if) , 

which, by setting x = ^/9\P — /3g\/2, returns 

253^ e{/3 - /3g)^ + O {x^) . 



(4.39) 



^ / /35 
+ ^ -253+ + 



(4.40) 



Therefore the entropy in the critical region (3 ^ Pg is given by 

-(/3 - /?,) + 3 (5 + (x/2 - 5)^(/5 - /5,)) 1/3 - ^gf/' 



11 , /3\ 3 
' - -12+^H2j+i' 

+ {{P - Pg)') . (4.41) 

The entropy and its first derivative are continuous at f3 — f3g, while its second 
derivative diverges at this point. See Fig. 4.8. Thus the critical exponent for the 
first derivative of the heat capacity is —1/2. 
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Figure 4.8: Entropy density (a) and its second derivative (b) versus ^. The entropy 
is continuous at /3g while its second derivative diverges. 



By expanding the free energy for ^ ^ we have 

3 3 ' 

M - 



^log^ + ^(/3-/3.)-^(/3-/3.)^ 



81v^ 



(/5 - figfl^ + O ((/5 - figf) 



(4.42) 
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If one relaxes the unit trace condition (3.8), our model coincides with one that 
has been studied in the context of random matrix theories (Morris [1991]). The 
objects generated in this way, thru an appropriate double scaling limit /? — >■ and 
N — )■ +00 correspond to chequered polygonations of 2D surfaces, a theory of pure 
gravity. Therefore, the constraint Tr p^i = 1 is irrelevant for the critical exponents 
in this region. This is why we call this solution the gravity distribution for the 
gas of eigenvalues. The second critical point for the system is given by ^ = — ^_|_ 
{S = 1), and corresponds to the restoration of the Z2 symmetry of the density 
function of the eigenvalues, namely from the asymmetric to the symmetric arcsine 
distribution. We already noticed that it represents the specular counterpart of the 
phase transition at /3+ = 2. Similarly to the gravity point, we expand Eqs. (4.26) 
and (4.33) for 5 ^ 1 

^ = -2 + 4(5 - 1) - (6 - e{6 - 1)){6 -If + O {{5 - if) , (4.43) 

and for /? — >■ — /?+ — —2 we have 

<^ = 1 + ^ + ^ (3 - \e{P + 2)) (/? + 2f + ((/? + 2)=^) . (4.44) 

We finally expand the entropy for /3 — >■ — /3+ 

5 = -^-ln2 + ^ + ^^-^(/3 + 2))(/3 + 2)2 + 0((/3 + 2)3). (4.45) 

The entropy together with its first derivative is continuous, while its second deriva- 
tive displays a finite jump (5/32) at /3 = — /3+ = —2. See Fig. 4.9. 




Figure 4.9: Entropy density (a) and its second derivative with respect to 8 (b). 
The entropy is continuous at /3 — — /3+ = —2 while its second derivative presents a 
finite discontinuity. 
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Let us briefly overview the behavior of the system for the scahng a — 3. We plot 
the interesting behavior of the minimum eigenvalue Amin — a — m — 5 (see Eqs. 
(3.9) and (3.10)) plotted in Fig. 4.10. 




Figure 4.10: Minimum eigenvalue a — m — 5 versus p. Solid line: stable branch. 
Dotted line: metastable branch. 

For — < /3 < a coincides with the origin (left border of the solution do- 
main). This variable signals both second order phase transitions at —2 = —/3+ 
and at 2 = The Z2 symmetry is broken for —/3+ < /3 < /3+. Notice, however, 
that the gravity critical point at (3g = —2/27 remains undetected by a, while it 
is detected by the vertical derivative displayed by the average and the width of 
the solution domain at Pg. See Fig. 4.11. However, since (3g and — /3+ lie on an 




Figure 4.11: Average m (a) and width of the sohition domain 6 (b) (Eq. (3.10)) as 
a function of /3. Solid line: stable branch. Dotted line: metastable branch. 

analytic continuation of the solution obtained for /3 > 0, we arc not assured that 
this is a stable branch. In the next chapter we will show that a first order phase 
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transition occurs at /? ~ — 2.4:55 /N in this scaling. In Fig. 4.12 we overview the 
results presented for the scaling a = 3 by plotting the internal energy u and the 
entropy density s as a function of the inverse temperature, —oo < f3 < oo. For 
/3 > we have a stable branch for the system (solid line), whereas for /3 < we 
will see that the solution is metastable. For the sake of future convenience let us 
record that at /?+ = 2 the internal energy density reads u = 5/4, at f3g = —2/27 
u — 9/4, and at —/3+ — —2 we have u — 7/4. 
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Figure 4.12: Internal energy density u (a) and entropy density s (b) versus /?. Solid 
line: stable branch. Dotted line: metastable branch. 



4.4 Entropy density as a function of average en- 
tanglement 

Let us now reinterpret the behavior of the system directly from the point of view 
of quantum information. In this section, we will study the behavior of the entropy 
as a function of the internal energy, that is the average purity (see Eq. (3.71)). Its 
expression for 3 < 5 < 2 + V2 and 1 < 5 < 2 + -\/2 is extremely involved. On the 
contrary for the other two ranges of 5 we have: 

(2^{u-l)/3, l<u<l 
5 = I (4.46) 
[3 + V9-4h, 2<u<l, 
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which plugged into Eqs. (4.21) and (4.37) give 

f|ln(ii-l)-|ln3-|, l<u<l, 

s(u) = 



(4.47) 



As can be expected, by analytic continuation we find that s (it) for 2 < it < 9/4 
has the same expression obtained for the interval 5/4 < u < 2 (see Eq. (3.69)). 
Let us now compute s{u) in the critical region unveiled in the last section. At fig 
{5 = 3) the average purity is (ttab) = u/N = 9/4:N. By expanding it for 5 ^ 3 we 
have 

9 



^=4 4 



\ (l - l^iS - 3)) {6 -Sf + O {{6 - 3f) (4.48) 



X 



and inverting the series we find 

5 = 3+ (2 + {2V2-2)e(u-^y. 

+0(u - ^ (-16x^ + n6V2x^ - 2272x^) + O {x^) , (4.49) 

where x = ^/\u - 9/4|. Prom Eqs. (4.21) and (4.28) we get that for u 9/4 the 
entropy is given by 

11 3 2 2 4 4 
s — hln X X 

12 2 27 81 

I 256 + (1504^2 - 256)/m - - ) ]x^' + (x^) . (4.50) 
3645 \ \ A J j 

Notice that even if we are exploring a metastable solution for the system, and 
thus the Maxwell relations may not be satisfied, at /3 = /3g we have ds/du — 
f3g. Furthermore, the second order phase transition at u — 9/4 is unveiled by 
a divergence for u (9/4)"*" of the third derivative of the entropy, d^s/du^ (see 
Fig. 4.13), while we have seen that its second derivative with respect to /3 diverges 
with different rates on both sides of Pg. This behavior can be justified by observing 
that the solution m — 5 — 3 found by minimizing the free energy in the physical 
domain, is also a stationary point for the system. This can be inferred from the 
plot in Fig. 4.3(a) of the free energy on the boundary of the domain, as m and S 
tend to 3 (see Fig. 4.2(a)). 
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Let us now consider the other critical point at p — {S = 1). We have 

+0{{6-lf) (4.51) 

4 / . / 7\\ / 7\^ / / 7^ 



+_^8 + 99^«-jjj^»-jj |. (4.53) 

See Fig. 4.14. Observe that in this case, the second derivative of the entropy with 
respect to the internal energy is discontinuous, as found by deriving s with respect 
to /?. 

Note incidentally that the point p — — 3/2 + y/2 (5 — 2 + y/2), connecting real 
solutions of Eq. (4.13), is not a critical point for the system. Indeed for 5 — > 2 + 
we get 

("4-^)' - ''^(';;;g;'' +o(p+v2-.n (4.54) 

{s-sof = ^(2 + V2-5)(3\/2 -4) + O ((2 + ^2-5)2^,(4.55) 
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Figure 4.14: Entropy (a) and its second derivative (b) versus p. 



where 

It then follows that the entropy has a linear behavior for /3 = —3/2 + \/2, that is 
(^-^o)= Q + (^P+l-V2^+o[{2 + V2-6r). (4.57) 

Indeed the divergences in the derivative of the entropy and in the derivative of the 
inverse temperature with respect to S at this point compensate. In Fig. 4.15 we 
show the behavior of the entropy for the metastable branch. 
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U 

Figure 4.15: Entropy density versus average entanglement u = N{ttab), for the 
metastable branch. The four solutions are plotted with different colors. In partic- 
ular both the red and the blue lines refer to the asymmetric arcsine distribution, 
and correspond to the real solutions of Eq. (4.14). 
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Let us briefly comment on the fact that the metastable branch which emanates from 
the analytic continuation of the solution at positive f3 described in the previous 
chapter has not led us towards separable states. See also Fig. 4.16 where we plot 
the entropy as a function of the internal energy. We recall that in this scaling 
{a = 3) the average purity is (tiab) = u/N, namely of 0{1/N). The eigenvalues 
remain of C (1/A^) (and so does purity) even though the temperature can be (very) 
negative (as /3 crosses 0). In order to find separable states we will have to look at 
the stable branch in the next chapter. We conclude that until now we have studied 
the region of the Hilbert space given by highly entangled states. 
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Figure 4.16: Entropy density versus average entanglement u — N{ttab) for the scal- 
ing a = 3: positive temperatures and metastable branch of negative temperatures. 
See Figs. 3.10 and 4.15. 



Chapter 5 



Negative temperatures: stable 
branch of separable states 



In this chapter we will search the stable solution of the system at negative temper- 
atures. We will find that for ;5 — > — oo the most probable solution for the Coulomb 
gas of eigenvalues involves the evaporation of one Schmidt coefficient from O {1/N) 
to 1. In other words, it will appear a new phase for the system related to the set 
of separable states. This solution will be associated to the appearance of a new 
critical point for the system, whose nature will be different from the phase transi- 
tions discussed in the previous chapters. In particular, by maximizing the potential 
energy of the Coulomb gas, we will find that at /3 = —2/N a competition starts 
between two solutions: one with all the eigenvalues of order O (l/N), and one with 
an isolated eigenvalue of O {1) and all the others of 0{1/N), namely the typical 
and separable phases. At /3_/N ~ — 2.455/A^ we will unveil the coexistence of 
these two phases, and the system undergoes a first order phase transition. Finally, 
hi /3 < P^/N the stable phase for Coulomb gas will be given by separable states. 
However, in order to yield the correct thermodynamic limit in this statistical anal- 
ysis, we will need to properly modify the scaling coefficient with respect to the 
stable solution of positive temperatures (and the metastable branch for negative 
temperatures). This will introduce a new scale in the inverse temperatures and, in 
particular, the phase transition in the new scaling will be at ~ —2.445. 

In Sec. 5.1 we will discuss the new scaling for the system, which will affect the 
saddle point equations, whose physical solution will be obtained in Sec. 5.2. In 
Sec. 5.3 we will look for the most probable distribution of the Schmidt coefficients 
by maximizing the free energy of the Coulomb gas, and then unveil the presence 
of the new separable phase. In Sec. 5.4 we will discuss the thermodynamics of the 
system and identify the first order phase transition from typical to separable states. 
We will conclude this chapter with a brief overview about the finite size corrections. 
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in Sec. 5.5: on the one hand we will explain our numerical results, on the other 
hand we will provide a mapping between the results shown in Chaps. 3 and 4 and 
the stable solution for negative temperatures described in the next sections. 

5.1 Natural scaling for /3 < 

From definition (2.7) of the partition function, one expects that for any N, as 
(3 — )■ — oo, the system approaches the region of the phase space associated to 
separable states: here the purity is C (1) and thus the scaling coefficient in Eqs. 
(2.7)-(2.8) is a = 2. From the point of view of the parameter /3, this can be 
interpreted as follows: by adopting the scaling A^^ for the exponent of the partition 
function, we will explore the region (3 = (l/N) of the scaling A^'^ introduced for 
positive temperatures. Notice that the critical point Pg = —2/27 for the solution 
at negative temperatures now reads (3 = —{2/27)N and escapes to — oo in the 
thermodynamic limit, as pictorially shown in Fig. 5.1. 




^ ► a = 2 



Figure 5.1: By setting the scaling coefficient a = 2 for /3 < we zoom on the region 
(3 = O (l/N) of the scaling (i.e. a = 3) introduced for positive temperatures. 
The critical point (3g = —2/27 now reads (3 = —{2/27)N and escapes to — oo. 

We will show that the solution (3.46), according to which all the eigenvalues are 
O [l/N), becomes metastable in the region of negative temperatures, and the most 
probable distribution of the eigenvalues is such that one eigenvalue is O (1): this 
solution in the limit P — )■ — oo will correspond to the case of separable states. 
By following an approach similar to that adopted for positive temperatures, we 
will first look for the set of eigenvalues {Ai, . . . , Xjy} satisfying the saddle point 
equations (2.12)-(2.13) with a = 2, and get, as in Sec. 3.3, a continuous family of 
solutions. We will select among them the set maximizing {/3 < 0) the free energy 
(2.11), with a = 2: 

fN= ^'-Wb ^ 1HA,-A.|. (5.1) 

l<j<Af ^ l<i<3<N 



|__ 

-OO iVA 
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Since we are approaching the hmit ^ — )■ — oo the states occupying the largest 
volume in phase space are separable. We define Ajv = as the maximum eigenvalue 
of order of unity, whereas the other eigenvalues are O (l/N): 

X^ = I^ = 0{1), Yl '^^ = 1-/^- (5-2) 

l<i<N-l 

From this it follows that we need to introduce the natural scaling only for the sea 
of the first N — 1 eigenvalues: 

A, = (1-/^)]^' °<*^^]V^-^' Vie{l,...,7V-i}. (5.3) 

In the limit iV ^ oo we will describe the behavior of the rescaled eigenvalues in 
terms of a continuous, non negative and normalized density function as discussed 
in section Sec. 3.1 (see Eqs. (3.3) and (3.4)). 



5.2 Saddle point equations 

In order to solve the saddle point equations for the expression (2.12)-(2.13), we 
will separately consider the contribution of the first A'" — 1 eigenvalues and of the 
isolated eigenvalue: 

^ ^ ^ l<j<N-l,j^i * ^ ^ 

(5.4) 

l<j<N-l 

In the N ^ oo hmit, by neglecting contributions of O (l/N), we get: 

2/x/3 + i^ =0, (5.8) 
Xp{X)dX = 1 (5.9) 



/ 

^0 



where p is the density of the first — 1 cigcnvahics and (5.8) is the condition 
deriving from the saddle point equation (5.5) associated to /i. By the same change 
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of variables introduced in Sec. 3.2, see Eqs. (3.9) and (3.10) 

\ = m + x6, (5.10) 

with m = [a + b)/2, 6 = {b — a)/2 and A G [a,b], the solution of the integral 
equation (5.7) can be expressed in terms of 4>{x) = p{X)S. It can be easily shown 
that the solution of this Tricomi equation is 

1 2x(m-l) 



(j){x) 



1 



6 



(5.11) 



where the normalization constraint (5.9) on the rescaled eigenvalues fixes the La- 
grange multiplier 



. 4(m- r 



(5.12) 



The region of the parameter space (m, 6) such that the density of eigenvalues is 
nonnegative reads 

' (5.13) 



max <5, 1 [><m<l, 

1 ' 2/ - - 2 

Notice that the physical domain for the sea of the first — 1 eigenvalues has the 
same expression of the domain found for the range of positive temperatures (3.20) 
at /3 = 0, namely Tf{6, 0) = 1 ± S/2 (see Fig. 5.2, which is the analog of Fig. 3.2). 




Figure 5.2: Physical domain for the sea of the first — 1 eigenvalues, for /3 < 0. 

This is consistent with the change in the temperature scaling from to iV^: we 
are "zooming" into the region near /3 — )■ 0~ of the range of temperatures ana- 
lyzed in Chap. 3. Furthermore, as could be expected from what seen for positive 
temperatures, the solution of the saddle point equations is a two parameter con- 
tinuous family of solutions. In the next section we will determine the distribution 
of eigenvalues that maximizes the free energy of the Coulomb gas of eigenvalues. 
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5.3 Thermodynamics of the Coulomb gas 

Will now set the basis for the thermodynamic characterization of the Coulomb gas 
of the Schmidt coefficient, when the scaling coefficient is a = 2. From Eqs. (5.1) 
and (5.2) we get 

= E ln|A.-A,| + o(l) (5.14) 

l<i<j<N ^ ^ 

By applying the natural scaling we have (5.3) 

In = /^'-^^^(l-/^)-]^ E ln|A(t.)-A(t,)| + ^lniV + o(^^^^ 



N 



1 /In TV' 

= f + -\nN + 0(—], (5.15) 



= M--s+-lniV + C 



where 



(3 (3 \ N 

1, ,^/lnA 



w = (5.16) 



ln(l-/.) + A \n\X{t,) - X{t,)\ (5.17) 

l<i<j<Ar-l 



Finally, in the limit N ^ oo 

u ^ (5.19) 
s = ln(l-/x)-/3/red(m,5,/3) (5.20) 
^/ = f3u-s^piJ,^-Hl-ij,)+Pf,,^{S,m,P), (5.21) 

being /red the reduced free energy of the sea of eigenvalues 

fred{S,m,/3) = J dx4>{x) J dy4>{y)lii{6\x -y\) 

- \J'-). (5.22) 



P \2 
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As shown for the range of positive temperatures and for the metastable branch, /, 
u and s will be the free energy, the internal energy and the entropy densities of our 
Coulomb gas when computed at its most probable distribution, i.e. on the solution 
{n,(f){x)} which maximizes the exponent of the partition function or equivalently 
minimizes /?/ (/3 < 0). 

5.3.1 Minimization problem: separable states 

It is easy to see that (3f^ed{m,6), has no stationary points, but only a global 
minimum /3/red = 1/2 at (5, m,/3) = (2,2), indicated by the red dot in Fig. 5.3(a). 
This point yields the Wishart distribution (3.45) or (3.46) found at /3 = for the 
case of positive temperatures: 

where one should remember that now the A's are also scaled by 1 — /i, see Eq. (5.3). 
We stress that this result is valid for all /3 < 0. In order to check this solution we 
have to compute the free energy on the boundary of this domain, see Fig. 5.3(b) 
(which is the analog of Fig. 3.4(a)). 




2 3 



U.U I ' ' ' ' • ' ' ' ' • ' ' ' ' • ' ' ' ' • r- 
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s 

Figure 5.3: (a) Contour plot of the reduced free energy density /3/red(^, '^) of the 
sea for negative temperatures, (b) Reduced free energy /3/red on the boundary of 
the triangular domain in Fig. 5.3(a) for the case of negative temperatures. Solid 
line: upper boundary; dashed line: lower boundary. 
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Thus by minimizing the potential energy for the sea of eigenvalues we get 

u ^ f/ (5.24) 
s = ln(l (5.25) 

/?/ = pi,^-ln{l-i,)+^-. (5.26) 

A new stationary solution, in which the largest isolated eigenvalue /i becomes O (1), 
can be found by minimizing f3f and yields 

which is defined for /3 < —2. This expression can be also obtained directly by plug- 
ging Eq. (5.12) into the saddle point equation (5.8) corresponding to the isolated 
eigenvalue /i. This eigenvalue, of O (1), evaporates from the sea of eigenvalues of 
O (l/N), as pictorially represented in Fig. 5.4. 




Figure 5.4: Evaporation of the eigenvalue /i — O {1) from the sea of eigenvalues 
0{1/N). 

The isolated eigenvalue moves at a speed —d^/d^ — l/(2-\//3^ + 2^^), which di- 
verges at ^ = —2: another symptom of criticality. However, this new solution, 
when it appears at /3 = —2, is not the global minimum of /?/: it becomes stable at 
a lower value of /3. This can be immediately seen from the plot, in Fig. 5.5, of the 
free energy (5.26) as a function of the detached eigenvalue ji for different values of 
p. When f5 > —2 there are no stationary points, but only a global minimum at 
/J, — 0. According to the solution corresponding to this point, n is still in the sea 
of the eigenvalues O {1/N) and the stable density function is given by the Wishart 
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Figure 5.5: Reduced free energy as a function of /i for different values of /3(< 0). 
Notice the birth of a stationary point for /3 = —2 (/x = 1/2) which becomes the 
global minimum for /3 < /3_ {fi > fi^). 



distribution (4.16) with the free energy (4.22) (remember that, in the zoomed scale 
considered in this chapter, f3g corresponds to the very large inverse temperature 
NPg). At (3 = —2 it appears a stationary point for the free energy density corre- 
sponding to fi = 1/2 = C(l), see (5.27). Notice, however, that (3f at this point 
remains larger than its value at the global minimum, until /3 reaches Finally, 
for /3 < /3_ the global minimum of /?/ moves to the right, to the solution containing 
H = O (1). The inverse temperature /3- represents a critical point for the system, 
which we will soon see, undergoes a first order phase transition. The critical tem- 
perature /3_ is the solution of the transcendental equation /(/3_,0) = /(/3_,/i_), 
that is 

^" --ln(l-^_), (5.28) 



2 1 



which yields 



0.71533, (3- 



2/i_(l - /i_ 



-2.45541. 



(5.29) 



The maximum eigenvalue is then a discontinuous function of the temperature at 

/3 = /3_ 

'O, 0</3</3_, 
/i(/3) = { (5.30) 

and in the limit (3 — )■ — oo, /i approaches 1: the state becomes separable. 
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See Fig. 5.6. We stress that Eqs. (5.23) and (5.27) are the distributions of eigen- 
values that gives the largest contribution to the partition function, that is the the 
most probable distribution for the isolated eigenvalue and for the sea of eigenval- 
ues, respectively. 




Figure 5.6: Maximum eigenvalue at negative temperatures. The two solutions are 
exchanged at /3_ —2.45541. The solution with fi = O (1) appears at /3 = —2 
( fi = 0.5) but is not the stable solution. Full line: solution of minimal /?/, i.e. 
maximal free energy; dashed thick line: solution of higher /?/. 

Thus, by inserting Eq. (5.27) in Eqs. (5.24)-(5.26), we determine the internal energy 
u, the entropy s and the free energy / densities for the system: 



for /3_ < /3 < 



u 




1 

2 
1 



(5.31) 
(5.32) 

(5.33) 



for /3 < /3_ 



u 



1 
2 

In 



1 1 

2/3 2 

1 1 

2 ~ 2 

2 2 



1 + 



/3' 



1 + 



/3 



1 

2' 



In 



1 + 



/3 



(5.34) 
(5.35) 
(5.36) 
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Summarizing, for /3_ < /3 < the solution of saddle point equations maximizing the 
free energy of the system is such that all eigenvalues are O {1/N), and through an 
appropriate scaling reduces to the Wishart distribution described in Sec. 4.2. We 
remark that, since the scaling is a = 2, we are neglecting the contributions deriving 
from the sea of eigenvalues, since they scale as O {In N/N) (see Eq. (5.15)). At 
/3 = —2 a new metastable solution appears with one eigenvalue of O (1), and for 
/3 < /3_ it becomes the stable solution, whereas the distribution of the eigenvalues 
described in Sec. 4.2 becomes metastable. It is important to emphasize that in the 
scaling a = 3, adopted in Chaps. 3 and 4, the interval of negative temperatures 
where the Wishart distribution refers to all the eigenvalues, reads P-/N < P < 0, 
and shrinks to zero in the thermodynamic limit (that is why the solution found in 
Chap. 4 is metastable). 



5.4 First order phase transition 

We are now ready to unveil the presence of a first order phase transition in the 
system. The branch (5.34)- (5.36) is stable for /3 < /3_ while it becomes metastable 
for 13- < (3 < —2. On the other hand, the solution /U = 0, corresponding to (5.31)- 
(5.33) has a lower value of [3f for /3_ < /3 < 0, and a higher one for (3 < (3^. See 
Fig. 5.7(a). At /3_ there is a first order phase transition, signaled by a discontinuity 
in the first derivative of the free energy, as shown in Fig. 5.7(b). 
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Figure 5.7: Free energy (a) and first derivative (b) with respect to /3. The two 
solutions are exchanged at /3_ —2.45541, where there is a first order phase 
transition. Full line: solution of maximal free energy; dashed line: solution of 
lower free energy. 
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Indeed for (3 ^ (3- we have 
1 



/3/ 




P) + ((/3 - /3_)3) 



with m(/3_) = /i^ ~ 0.5117 and 



du 



/3- 




-0.2755, 



(5.37) 



(5.38) 



see Eq. (5.34). We have seen that at /3_ the free energy of the global minimum at 
/X = is equal to the free energy of the stationary point at = Thus at this 
temperature our statistical model is characterized by the coexistence of two phases, 
corresponding to typical and separable states. In Fig. 5.8 we plot the entropy and 
the internal energy densities versus (3. 
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Figure 5.8: Internal energy (a) and entropy density derivative (b). At fixed inverse 
temperature /3_ the internal energy goes up from mt = to ~ 0.5117 whereas 
the entropy goes down from st = —0.5 to ss — —1.756. 



Observe that at fixed temperature, /3 = /3-, the internal energy of the system goes 
from mt = up to us = fi'^ — 0.5117, while the entropy goes from st = —1/2 down 
to Ss = — l/2 + ln(l — /i„) ^ —1.75643: the system moves from typical to separable 
states. We have that u G [ut, us] and s G [ss, st] are given by convex combinations 
of the boundaries of the associated intervals and describe the coexistence of the 
two phases. More precisely, if our statistical ensemble consists of n identical copies 
of the system, we say that during this first order phase transition a certain amount 
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riT of them is in the typical phase, that is with the largest Schmidt coefficient null 
on average (i.e. of O (l/N) in the rescaling a = 3), and a complementary amount 
ns = n — riT is in the seperable phase, the average of ^ being Notice that 
there can be some (few) copies of the system whose free energy is not maximal, in 
this sense we should have written n ~ nr + ns- However, these fluctuation can be 
considered negligible for our scopes. Thus, as the energy of the system increases 
we have a gradual transition from the typical to the separable phase. In particular 
for a generic state X, when we have both phases, the internal energy u and the 
entropy s are 



0.5' 





utUt + usns 

u = 

riT + ns 


and 


stUt + ssns 

s = . 

riT + ns 


(5.39) 
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Figure 5.9: First order phase transition, at fixed temperature /3 = In (a) and 
(c) the system is the typical and separable phase, respectively. In (b) the two 
phases coexist. 



Since this first order phase transition refers to a reversible isothermal process, we 
have 

/ ds = 13^ du (5.40) 



86 



Negative temperatures: stable branch of separable states 



thus 



We also have 



s = St + (3-{ut ~ u) = —- + (3-u. 



As 
Am 



/3- 



(5.41) 



(5.42) 



with As = St — Ss and Am = ut — us- The variation Am of the internal energy can 
be interpreted as the specific latent heat of the evaporation of the largest eigenvalue 
from the sea of the eigenvalues, from 0{1/N) up to 0{1). Summarizing, the 
entropy density as a function of the internal energy density for the stable solution 
at negative temperatures reads 



I3.U 



2' 



S[U 



ln(l - v^) 



1 

2' 



< M < /i^ 



/i^ < M < 1. 



(5.43) 



It is continuous together with its first derivative at n = /x^, while its second 
derivative shows a finite gap. See Figs. 5.10 and 5.11. 
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Figure 5.10: Entropy density s versus internal energy density u = (itab)- See Eq. 
(5.46). 



At the critical point, for u ^ fi^, we have 



s{u) = Ss + P-{u - Us) 
2 



-2 + /3_ 



'2 + /3_ 



/3- 



(u - usfeiu -us) + {{u - Usf) , 

(5.44) 
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where 



d/3 
du 



'2 + /3_ 



/3- 



-3.629, (5.45) 



see Eq. (5.34). 



-3 



cS?-4 



(a) 



U 




u 



Figure 5.11: Entropy density (a) and its second derivative (b) versus u. 



A few words of interpretation are necessary. As we have seen, it is the stable branch 
of the solution that lead us to separable states at negative temperatures. The 
analytic continuation of the stable solution for positive temperatures emanates a 
metastable branch in which all eigenvalues remain O (l/N). By contrast, the new 
stable solution consists in a sea of — 1 eigenvalues 0{1/N) plus one isolated 
eigenvalue O (1). Let us now translate these results in terms of purity (we stress 
again that ^ is a Lagrange multiplier that fixes the value of the purity of the 
reduced density matrix of our N dimensional system). The average entanglement 
of the states described by the separable phase is given by 



l<j<N 



1 



/ij 



A^ 



46) 



Recall the scaling (1 — //)/( A" — 1) for the sea of eigenvalues. It is the analog of 
(3.71) for the scaling a = 3. Assume that we pick a given isopurity manifold in the 
original Hilbert space, defined by a given finite value of the average purity (hab)- 
If we randomly select a vector belonging to this isopurity manifold, its reduced 
density matrix (for the fixed bipartition) will have one finite eigenvalue ~ ^/ttab 
and many small eigenvalues 0{1/N) (yielding a correction 0{1/N) to purity). 
In this sense, the quantum state is largely separable. The probability of finding 
in the above mentioned manifold a vector whose reduced density matrix has, say. 
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two (or more) finite eigenvalues /^i and 112 (such that /xf + ~ (t^ab), modulo 
corrections O (l/N)) is vanishingly small. By contrast, remember (from the results 
of Chap. 3) that if the isopurity manifold is characterized by a very small value 
O {1/N) of purity, the eigenvalues of a randomly chosen vector on the manifold are 
all O {1/N) (being distributed according to the semicircle or Wishart, depending 
on the precise value of purity) . This is the significance of the statistical mechanical 
approach adopted in this thesis. We will come back to this point in Chap. 6. 

5.5 Finite-size corrections 

The results of the previous section refer to the N ^ 00 hmit. In order to under- 
stand how finite-A^ corrections affect our conclusions we have numerically max- 
imized the free energy for various temperatures. The two phases of the system 
discussed in the previous sections correspond to the two solutions obtained by 
minimizing /S/jv (5.1) on the N dimensional simplex of the normalized eigenval- 
ues. Indeed, we have numerically proved that pfj^{p) presents two local minima 
at negative temperatures: for /^i^'' < /3 < the minimum giving the lower value 
of /3/iv(/3) corresponds to the distribution of eigenvalues (4.16), found in the last 
chapter. The other minimum is reached when the highest eigenvalue is O {1). The 
point (3 = /3^^ is a crossing point for these two solutions, and for /3 < /3^^ these 
two solutions are inverted, see Figs. 5.6 and 5.7. Summarizing, there exists a neg- 
ative temperature at which the system undergoes a first order phase transition, 
from typical to separable states. The first thing to notice is that qualitatively the 
phase transition remains of first order even for finite A^. The second is that the 
finite- corrections are quite relevant for the location of the phase transition and 
the value of the maximum eigenvalue as a function of (3. For example, for N — 30, 
the negative critical temperature = —1.935 instead of —2.455. This is evinced 
from Fig. 5.12, which is the finite-size version of Fig. 5.7. This can be understood, 
as the corrections to f{fi) around /i = are quite large. In the limit ^ = 1/N there 
is a hard wall for the maximum eigenvalue 11, as the condition Aj = 1 cannot be 
satisfied \l ji <1/N. It is therefore hkely that all sorts of large corrections occur 
as jJL tends to 1/N, probably yielding an effective size to the corrections which is a 
lower power oi 1/N (or even possibly l/lnN). The limits — >■ and A?" — >■ 00 do 
not commute. 

To further explore this effect we have minimized 13 with respect to Ai, . . . , Xn-i, 
for fixed values of the largest eigenvalue = /j, and for different temperatures. 
The results for N — SO are shown in Fig. 5.13. One can see that between /3 — —1.8 
and p — —2.5 there is a competition between two well defined local minima, cor- 
responding to the two solutions discussed above. 
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0.16 

0.104 

/5/iv-lniV 

0.05 




/3 /3 

Figure 5.12: Finite- version of Fig. 5.7. Free energy and maximum eigenvalue in 
the saddle point approximation as function of /3 at = 30. The local minimum 
is in blue, the global one in red. The two minima swap stability at /3 = —1.935. 
Notice the birth of the new local minimum at /3 = —1.8 (for N = oo this takes place 
at /3 = —2) and the exchange of stability at /3 = —1.93 (for N = oo, /3 = —2.45). 



N = 30 




0.2 0.4 0.6 0.8 1 



Figure 5.13: Finite- iV version of Fig. 5.5. /S/at- In as a function of the maximum 
eigenvalue /x, obtained by numerical minimization over the remaining — 1 eigen- 
values for various /3. Observe the formation of a new minimum and the exchange 
of stability, although the critical values of /3 at which these phenomena occur differ 
from the theoretical ones, due to large finite A^ corrections. However, it is clear 
that at small /i, 1/A^ corrections tend to increase the value of /3/jv, making the 
critical value /9_ move towards 0, as observed in the numerics. 
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At p — P_ — —1.935 their free energies are equal. For higher /5 the global 
minimum corresponds to the solution (4.16), whereas on the other side of P^^'' 
the solution with // = O (1) minimizes Pf^. We have seen that for /3 > the 
stable solution has no detached eigenvalues. By taking into account the scaling 
P ^ P/N we get that the solution is given by the very first part of the gravity 
branch (4.22). In particular, recalling the natural scaling (5.3) for eigenvalues, the 
maximum eigenvalue is given by b/N = {m + 5)/N = 26 /N. On the other hand 
for P < P- the maximum eigenvalue is given by (5.27). Therefore, we get 

US{P/N), P-<P<0, 
Pi = I (5.47) 

where 6{P) can be retrieved from (4.17). The numerical results for = 40 arc com- 
pared with the expressions in Eq. (5.47) in Fig. 5.14. The agreement is excellent. 

N = 40 

ii 1 



0.5 



-4 



-3 



-2 



/3 



-1 



Figure 5.14: Maximal eigenvalue in the saddle point approximation as function of 
/3. The points are the result of a numerical evaluation for iV = 40, while the full 
line is the expression in Eq. (5.47). 

The corresponding free energy follows from (5.36) and (4.22) with the appropriate 
scaling 




mm 



/3- < /3 < 0, 
/5</3_. 



(5.48) 
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In the limit N ^ oo we get (3f = 1/2, as expected from Eq. (5.33). Notice that in 
order to have a correct finite size scahng of the critical temperature ^i^^ one should 
take into account O (l/N) corrections to the expression of /3f, see Eq. (5.14). This 
would enable a more exhaustive mapping between the two branches. This analysis 
goes beyond our scopes. In the next chapter we will overview the results of our 
statistical approach to bipartite entanglement, trying to connect what found in the 
two scaling regimes a — 3 and a — 2. 



Chapter 6 
Overview 



In this chapter we will summarize the main results obtained in the first part of 
this thesis, dealing with the analysis of bipartite large quantum systems in pure 
states. We will focus on those quantities that are more directly related to physical 
intuition. In the classical statistical mechanics approach we have adopted, the 
temperature plays the usual role of a Lagrange multiplier, whose task is to fix the 
value of energy (purity in our case). A given value of /3 determines a set of vectors 
in the projective Hilbert space whose reduced density matrices have a given purity 
(isopurity manifold of quantum states). The distribution of the eigenvalues of 
(the reduced density matrices associated to) these vectors has been investigated in 
the last three chapters and yields information on the separability (entanglement) 
of these quantum states. The distribution of eigenvalues is the most probable one 
(Mehta [2004]), in the same way as the Maxwell distribution of molecular velocities 
is the most probable one at a given temperature. In this chapter we will abandon 
the temperature and fully adopt purity as our physical parameter. In Sec. 6.1 
we will determine the probability distribution of the purity which is proportional 
to the volume of the isopurity manifolds. Then in Sec. 6.2 we will overview the 
different shapes of the distribution of eigenvalues for different phases of the system, 
corresponding to different ranges of tiab from 1/N to 1. 



6.1 Entropy and volume of isopurity manifolds 

Entropy counts the number of states with a given value of purity and is propor- 
tional to the logarithm of the volume in the projective Hilbert space. The explicit 
expressions of the entropy density s as a function of the purity ttab of the state 
vectors in that volume, can be read directly from Eqs. (3.69) and (5.43) for the 
stable solution, by taking into account the correct scaling when the system moves 
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across different regions of the Hilbert space. In order to discuss this dehcate and 
important point let briefly recall the expression of the partition function (2.7): 



where N is the dimension of subsystem A in the state pa- The power A^"" can be 
read as N°' — N^N"'"'^, where A^^ is the number of degrees of freedom of pA- The 
scaling coefficient a depends on the degree of entanglement of the global system, 
or equivalently on the temperature /3. In the region of positive temperatures, 
the purity of the system scales as 0{1/N), and becomes an intensive quantity 
if multiplied by N'^~^ — N, that is a = 3. The internal energy density is then 
given by u — Nttab- On the other hand, when we approach the region of separable 
states, (3 — > — oo, the purity is O (1), from which it follows that a — 2 and u — ttab- 
Thus, the entropy density (Eqs. (3.69) and (5.43)) reads 




(6.1) 
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s{71ab) = < 
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with //^ ~ 0.5117 and /3_ ~ —2.455, given by Eqs. (5.28) and (5.29). By exponen- 
tiating expression (6.2) we get the volume V = exp {N^s) (i.e. the probability) of 
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the isopurity manifolds 



1 / -rr ^5 
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X exp 



(6.3) 



2^ 
N 



e (1 - ^vfl^) 



2 < TTAS < 1. 



We have three main regions of entanglement, corresponding to the maximally en- 
tangled phase < t^ab ^ 5/4), the typical phase {5 /AN < ttab < A*?.) and 
the separable phase < ttab < 1)- Notice that ttab — 2/N refers to typical 
states. f3 = 0. at the boundary between the solution for the positive and the stable 
solution at negative temperatures. We have already noticed in Sec. 5.5 that the 
thermodynamic potentials show large corrections in the region /j, = O (1/iV). In 
particular, if on the one hand we have neglected contributions of O (l/N) to the 
potential energy density (5.14) in the region of negative temperatures {a — 2), on 
the other hand the average purity of the system in the region of positive tempera- 
tures (a = 3) is also of O (l/N) (see Eq. (3.71)). In other words when we derive 
the expression of the entropy as a function of the internal energy we deal with the 
same order of fluctuations in u and s. It follows that in order to yield a uniform 
connection between the two branches at ttab — 2/N, we should have considered 
higher order corrections in l/N. However, this goes beyond our scopes. In order 
to plot s{u) for finite N, notwithstanding the presence of these fluctuations, we 
have artificially modified the entropy as 



s{N ex.p{x)), 



X < \n{2/N) 



s{x) 



(6.4) 



s((exp(x) -2/N), x> \n{2/N). 



This function is continuous in x = In for finite N, see Fig. 6.1. The presence of 
discontinuities in some derivatives of the entropy detects the two phase transitions 

for the stable solution. Notice that, as already pointed out, discontinuities of the 
n-th derivative of s{T) translate in discontinuities of the {n + l)-th derivative of 
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Figure 6.1: Volume V = exp(A^^s) of the isopurity manifolds versus their purity 
tvab for N = 50. We indicate in blue, green and red maximally entangled, typical 
and separable phases, respectively. The black line corresponds to the first order 
phase transition from the typical to the separable phase. The dashed line refers to 
the metastable branch, see Figs. 4.15 and 4.16. 

s{u). At ttab = 5 /AN there is a second order phase transition signaled by a 
discontinuity in the third derivative, see Sec. 3.5 (Eq. (3.70)). The first order 
phase transition, which takes place between ttab = 2/N and ttab = /U^ ~ 0.512 is 
signaled by discontinuities in the second derivative of the entropy at those points. 
Actually, differently from tiab = f^-, analyzed in Sec. 5.4 (Eq. (5.44)), the same 
check can not be explicitly performed for hab = 2/iV, due to the large fluctuations 
in this region, mentioned above. 

Observe that entropy is unbounded from below: at both endpoints of the range 
of purity, hab = 1/A^ (maximally entangled states) and ttab = 1 (separable states), 
when the isopurity manifold shrinks to a vanishing volume in the original Hilbert 
space, the entropy, being the logarithm of this volume, diverges, and the number 
of vector states goes to zero (compared to the number of typical vector states) 



6.2 Distribution of eigenvalues 

In the previous chapters we have stressed that the presence of phase transitions 
in the system can be easily read out from the behavior of the distribution of the 
Schmidt coefficients, i.e. the eigenvalues of the reduced density matrix pA of one 
subsystem. From Eq. (3.68) we get the expression of the minimum eigenvalue 



96 



Overview 



a — m — 5 as a function of ttab 



^ (1 - 2VNnAB - 1) , ^<7rAB< 



A. 



.0, 



4Ar' 



(6.5) 



< TTAB < 1- 



See Fig. 6.2. We recall that m = 1 for 1/N < ttab < 5/4A^ and m = 5 for 
5 /AN < TiAB < 1, Eqs. (3.37) and (3.38). The second order phase transition at 
TiAB — 5/4 A^, associated to a Z2 symmetry breaking, is detected by a vanishing 
gap. 




Figure 6.2: Minimum eigenvalue as a function of ttab- At ttab = 5/4N the gap 
vanishes. 
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Figure 6.3: Maximum eigenvalue as a function of tvab- The dashed line in (a) 
represents the first part of the metastable branch, from ttab = 2/iV to ttab = 9/4:N 
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On the other hand, the maximum eigenvalue Amax coincides with the upper edge of 
the sea of eigenvalues b — m + 5, as given by (3.68), until it evaporates according 
to Eq. (5.34). Thus, 



Amax — ^ 



^ (1 + 2VN7rAB - 1) , ^ < TTAS < 

|(3-2^|-7V7r^s), 



5 

4Ar' 



AN 



N 



as shown in Fig. 6.3. 
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Figure 6.4: Summary (stable branch): profiles of the eigenvalue density for < 

T^AB <^,^<T^AB< jj, and I < TTylB < 1. 

Notice that due to the double scaling in the maximum eigenvalue and the purity 
in the two nearby intervals, 1/N < ttab < S/A?" and 2/N < ttab < 1, we pre- 
ferred to separate the plot of Amax vs ttab in two distinct figures. Fig. 6.3(a) and 
(b). In Chaps. 3 and 4 we have seen that in the different phases the distribu- 
tion of the eigenvalues of pa has very different profiles. See Fig. 6.4. While for 
1/N < ttab < the eigenvalues (all 0{1/N)) follow Wigner's semicircle law 
(maximally entangled phase), they become distributed according to Wishart for 
larger purities, ^ < ttab < jf (typical phase), across the second order phase tran- 
sition. This is a first signature of separability: some eigenvalues vanish and the 
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Schmidt rank decreases. For even larger values of purity, jj < ttab < 1, across the 
first order phase transition, one eigenvalue evaporates, leaving the sea of the other 
eigenvalues O (l/N) and becoming O (1) (separable phase). This is the signature 
of factorization, fully attained when the eigenvalue becomes 1 at hab = 1- 

In the previous chapters, we have associated to each phase transition in the 
system, a critical exponent, related to the behavior of the entropy as a function of 
the inverse temperature /3. See the table in Fig. 6.5. 









order 
of derivative 


jump 


critical 
exponent 


stable 
branch 


0+=2N 


5 

4iV 


2 


1 

8 





P-(l-M-) 







/?- 




mestable 
branch 


27 


9 

In 


2 


15-3\/2 


1 

~2 


-/3+ = -27V 


7 

iiV 


2 


5 

32 






Figure 6.5: Table of the critical points for the system, for the stable and metastable 
branches. In particular the third and the forth columns refer to the behavior of 
the entropy density s at the critical point while the last column refers to derivative 
with respect to /3 of the heat capacity. Notice that the first order phase transition 
at /3_ ~ —2.455 involves a broad interval of purity, from 2/N to /x^ ~ 0.5117, 
associated to the coexistence of typical and separable phases. 

For the first order phase transition, we have replaced the critical exponent with 
the thermodynamic relation (5.42). Note that in this case, at fixed temperature 
/3_ the purity, together with the entropy of the system, varies in a finite interval 
associated to the coexistence of the typical and separable phases. 

We conclude the first part of this Thesis with an overview in Fig. 6.6 of the 
phenomenology of the phase transitions in bipartite entanglement. Observe that 
while the second order phase transitions correspond to a given value of the purity, 
the first order phase transition from typical to separable states lies on a broad 
interval of O (1), namely 2/N < ttab < A*?. — 0.5117. We have also included the 
two metastable branches. A metastable branch is born at ttab — 2/N: it starts as 
Wishart, undergoes a second order phase transition at hab = 9/4 (2D gravity), 
where a singularity is developed at its right edge, then its support starts decreasing. 
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undergoes a second order phase transition at tcab = 7/4N (Z2 symmetry restora- 
tion) and eventually becomes sharply peaked (with two singularities). The other 
metastable branch is born at vr^^ = 1/4 and corresponds to the separable phase: 
it becomes stable at ttab = A*?.- 




0{1/N)^ I ►O(l) 

1st order phase tr. 

Figure 6.6: Overview of the "evolution" of the eigenvalue densities as a function of 
purity (1/A^ < ttab < 1)- The straight segment represents the stable branch. The 
three phases of maximally entangled, typical and separable states are indicated 
in blue, red and green, respectively, as done in Fig. 6.1 for the volumes of states. 
The curved lines refer to the two metastable branches. The diamonds indicate the 
three second order phase transitions. 
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Summary 



In the first part of the thesis we have introduced a statistical approach to the 
study of the behavior of bipartite entanglement of a large quantum system in a 
pure state. In particular we have chosen as entanglement measure the purity of one 
of the two subsystems, namely the local purity of the global state. Nevertheless, 
pure states represent only a special set of states as they refer to isolated quantum 
systems. In practice physical systems are described by a statistical mixture of 
pure states due to the interaction with the environment. In this second part of the 
thesis we intend to generalize our statistical analysis to the case of arbitrary mixed 
global states. Let us first observe that in order to study the properties of bipartite 
entanglement for a large quantum system in a mixed state, we should substitute 
the purity of one of the two subsystems with its convex-roof counterpart, as seen 
in Chap. 1. This considerably enhances the difficulty of the analytical treatment of 
our problem and will not be discussed in this thesis. On the other hand, a related 
interesting topic is the analysis of the statistical distribution of the local mixedness 
of the global density matrix, as an indicator of the average lower bound for the 
bipartite entanglement (see Eq. (1.11)). This problem will be tackled by studying 
the moments of the local purity of the system. 

This part of the thesis consists of three chapters. In Chap. 7, by following a 
standard approach in statistical mechanics, we will determine the formal expression 
of the moments of the purity at a given temperature. In fact, for the case of 
arbitrary mixed global states, due to the nonfactorization of the integral over the 
unitary group, we cannot directly compute the partition function by applying the 
method of the steepest descent, as was done for the case of pure states. In Chap. 8 
we will introduce a Gaussian approximation for the elements of the unitary group 
acting on the global and local systems, in order to determine the first moment of 
the purity at infinite temperature (/3 = 0), at 0{1/N), being the dimension 
of the subsystem whose purity distribution we aim to study. This will lead to 
a first generahzation of the results achieved for the case of pure states in (Page 
[1993]; Lloyd and Pagels [1988]; Facchi, Marzolino, Parisi, Pascazio and Scardicchio 
[2008]). Finally in Chap. 9, by exploiting the symmetry properties of the twirling 
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maps (VoUebrech and Werner [2001]) and thru the solution of some basic integrals 
over the unitary group due to Zuber (Zuber [2008]), we will compute the exact 
expression for the first two cumulants of the local purity at /3 = 0. Our results, 
when reduced to the case of pure states, will be in agreement with (Scott and 
Caves [2003]; Giraud [2007a]). We will conclude our analysis by determining the 
high temperature expansion of the first moment of the local purity. 



Chapter 7 

Statistical distribution of local 
purity 

We will devote this chapter to the introduction of the physical and mathematical 
setup that generalizes to the case of mixed quantum states the statistical approach 
introduced in Part 1 for the study of the local mixedness of pure states. In Sec. 7.1 
we will introduce the concept of local purity for a generic mixed bipartite state. 
It will play, in Sec. 7.2, the role of energy of the system in the definition of the 
partition function for a generic set of states with a given value of the global purity. 
Finally in Sec. 7.3 we will write the formal expression for the moments of the 
local purity, with respect to the canonical ensemble, in order to characterize its 
statistical distribution. 

7.1 Entanglement and local mixedness 

Consider a bipartite system given by two subsystems A and B, described by the 
Hilbert space "H = Ha ^'Hb, with dimTiA = N, dimTiB = M and dimTi = L = 
NM. Without loss of generality we will assume that N < M. In this second 
part of the thesis we will slightly change the notation introduced in the first part. 
In particular we will define the total Hilbert space as Hx, being X = AB our 
bipartite system, and will denote with Aj, 1 < i < L, the spectrum of the global 
states in "Hx, while the eigenvalues of each subsystems will be indicated by Xa^i 
and As J, for 1 < i < and 1 < j < M. (Recall that in the first part of the thesis 
the symbol \, 1 < i < N , referred to the eigenvalues of the states belonging to 
Ha). 

In the previous chapters we have explored the total Hilbert space l-ix from 
the point of view of bipartite entanglement between A and B, that is in terms of 
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Statistical distribution of local purity 



manifolds with a fixed value of the purity of one of the two subsystems (isopurity 
manifolds) . The main result of our analysis was the computation of the probability 
of finding a random pure state (of a large quantum system) with a given value 
of average entanglement between A and B, if sampled according to the unique 
invariant Haar measure. 

As mentioned in Chap. 1, in general the states of X are represented by the set 
of nonnegative unit trace operators (density matrices) p which act on the Hilbert 
space Hx, and reduce to (projections onto) vectors in the Hilbert space Hx if 
and only if they have only one nonnuU eigenvalue (equal to 1). We indicate with 
&{T^x) the set of all the states of X and with 

nx{p)^Ttp' e[l/L,l], (7.1) 

the purity of an arbitrary state p. The latter quantity provides a characterization 
of the global mixedness of the system and induces a partition of &{T-Lx) into a 
collection of distinct subsets 

Gxinx) = {p e e{nx) ■■ Trp^ = x}. (7.2) 

The minimum value of x = 1/L is attained when X is in the completely mixed 
state I/L, whereas the maximum a; = 1 is attained over the set &i{T-Lx) which 
includes all pure states IV')^^ . For each p e &{'Hx) the quantyties 

7rA(p)=Trpi and 7rB(p) = IVp^ (7.3) 

represent its A-local and B-local purity functions, pa = Tr^p and ps — Tr^p being 

the reduced density matrices of the subsystems A and B, respectively. Notice that 
the local purities of p can be different in general, and vary in different ranges: 

7rA(p) e [1/7V,1] and 7rB(p) e [1/M, 1]. (7.4) 

Clarifying the connection between the global and local purities of a state is an im- 
portant problem of quantum information theory (Bose and Vedral [2000]; Adesso, 
Serafini and Illuminati [2004]; Mazzola, Bellomo, Lo Franco and Compagno [2010]) 
which is closely related to the characterization of bipartite entanglement (Giihne 
and Toth [2009]; Bruss [2002] b). In particular we have seen that for the pure states 
p — \'il^)x (■01 of belonging to the set &i{Hx), due to the Schmidt decomposition 
(1.5), the A and the S-local purities coincide 

ttabW = Mmx (V'l) = M\i^)x (V'l) (7.5) 

and provide a measure of the bipartite entanglement between A and B: the smaller 
t^ab{'4') is, the larger is the entanglement contained in \'4>)x (V'l, see Sec. 1.2.2. The 
double index of "Kab refers to the equality between the local purities when the 
global state is pure (see Eq. (7.5)). 
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7.2 The statistical approach: partition function 

In the first part of tlie tliesis, we liave cliaracterized tlie bipartite entanglement 
between A and B by studying tlie statistical distribution of TTABi'ip) on ©i('Hx)- 
In this section we will generalize this statistical approach, based on the canonical 
ensemble, to the case of generic states, being the pure states only a special class of 
them. 

Let us briefly recall what seen for the x — 1. We have introduced the partition 
function 

Zab{P)= J dfi{ij)e-f'^^^^^\ (7.6) 

where the local purity TTAsii^) of plays the role of the internal energy density of 
the system and ^ is a Lagrange multiplier which fixes the value of the purity/energy 
selecting an isopurity manifold (Sinolecka, Zyczkowski and Kus [2001]). See Eq. 
(2.7) for comparison: here the coefficient iV", related to the scaling properties of 
the system, has been reabsorbed in the definition of the fictitious temperature. 
We also have that, d/i(V') is a (normalized) measure on the space of pure states 
&i(T-tx)- From the structure of complex vector space of Hx, wc have that the 
natural choice for dulip) has been induced by the Haar measure duniU) on the 
unitary group U{'Hx) — i^{L), through the mapping \ip)x = Ux \ipo)x with 
\ipo)x an arbitrary reference unit vector of T-Lx- Noticing then that for every 
Ux e ^{L) the reduced density matrix pA — Trs(|V')x (V'l) can be written as 
Pa = Ua^aUa, with Ua G V({N) and = diag{Ayi,i, ^a,2, ■ ■ ■ ^a,n}, the measure 
induced over the density matrices pA by duniUx), factorizes into the product of 
a measure over the unitary group d^H{U a) (related to the eigenvectors of pa) and 
a measure (ia(A^) over the (iV — 1) dimensional simplex of its eigenvalues A^j, 
1 < i < (Zyczkowski, Horodecki, Sanpera and Lewenstein [1998]; Zyczkowski 
[1999]). Thus the partition function becomes 

Zab^ = j dpH{UA) J da{KA) e'^^^^^-) = j da{KA) e'^^^^^^), (7.7) 

where 

da{KA) ^ Cn,M ^{l- ^ \A,i^ ^{^A,j)>^Aj^ Y\. - >^A,myd^XA , 

l<i<N l<j<N l<l<m<N 

(7.8) 

is the measure on the spectrum, with 9{x) the unit step function and Cjv,m a 
normalization coefficient, see Eqs. (2.5) and (2.6). 

A natural question is what happens when the system X is in a mixed state 
p of purity x ^ 1, rather than in a pure state. A generalization of Eq. (7.6) is 
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obtained by replacing tiab{'4^) with (say) the A-local purity 7rA(p) of Eq. (7.3) and 
the measure dii{ip) with a proper measure dfj,x{p) on the set Sxi'Hx)- This yields 
the following definition of the partition function of the A-local purity 



Za{x, = / dfi,{p) '^^("^ = / dti{p) 5(Trp2 _ x) e"/^ '^^(^^ , (7.9) 



where dp^p) is a probability measure on the set of mixed states (see below), and 



is a normalization factor. An analogous expression for the fi-local purity partition 
function Zb{x,I3) is obtained by simply replacing TTAip) with 7rs(p) in Eq. (7.9). 
Notice however that for x ^ 1, Zb{x, j3) will in general differ from Za{x, /3). More- 
over, it is worth stressing that the function Za{x, (3) provides only statistical infor- 
mation on the local mixedness of X, but not directly on its bipartite entanglement 
properties: this is due to the fact that for generic mixed states p of X the local 
purities /T^(p) and /Ts(p) are not entanglement measures, see Sec. 1.3.2. We have 
also seen that a generalization of Eq. (7.6) that retains the ability of characterizing 
the statistical properties of the bipartite entanglement of X for x 7^ 1 could in 
principle be constructed by replacing 7rA(V') with its convex-roof 



where the maximum is taken over all ensembles {pj, \ipj)AB}j which provide a 
convex decomposition of the mixed state p (e.g. see (Bennett, Di Vincenzo, 
Smohn and Wootters [1996]; Cornelio and de Oliveira [2010])). Eq. (7.11) coincides 
with Eq. (1.11) being Ep — tta- The quantity 7rA{p) is a proper measure of the 
bipartite entanglement, but the resulting partition function does not allow for a 
simple analytic treatment and will not be discussed in this thesis. The average local 
purity we will determine with our statistical approach can be seen as the average 
lower bound for the bipartite entanglement of the global system, with respect to 
the bipartition AB. 

Finally, another delicate point is that, differently from what seen for the case of 
pure states, there is no unique measure on mixed states (Zyczkowski and Sommers 
[2007]). Indeed, as previously discussed for pA, the Hermitian matrix p can be 
always diagonalized by a unitary rotation, and as consequence we can write the 
measure dp{p) as the product of a measure on the (L — 1) dimensional symplex 
of the eigenvalues and a measure of on the unitary group W(L) related to the 
eigenvectors. However, if on the one hand it would be natural to take the Haar 




(7.10) 




(7.11) 
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measure on U{L) so that = diJ,{UxpUx), on the other hand the measure on 

the eigenvalues can be chosen in different possible ways (Slater [1999]). Thus, we 

need to properly specify the choice of dfi{p), since the consistency requirement that 
for pure states Eq. (7.9) should reduce back to Eq. (7.6), i.e. Z^(l,/3) = Zab{P), 
does not eliminate such ambiguity. In order to overcome this point we will use a 
balanced purification strategy. In other words we introduce a composite Hilbert 
space 

-Hxx' ^V-x® nx', with -Hx ^ ^x' (7.12) 

that is T-Lx and T-Lx' are isomorphic Hilbert spaces. In this dimensional Hilbert 

space, each p of X can be represented by those pure states 1^)^^' which provide 
a purification for such density matrix, i.e. which satisfy the identity 

p^Trx'mxx'm- (7.13) 

Thanks to this identification we can now induce a measure on &{1-Lx) by sampling 
the pure states on Hxx' according to the unique, unitarily invariant Haar mea- 
sure which, as usual, is induced by the Haar measure on the unitary group U{L^) 
through the mapping \'^)xx' ~ ^xx' |^o)xx'' where |^o)xx' arbitrary ref- 

erence vector and Uxx' G U{L'^). Thus, with this choice the partition function 
becomes 

ZA{x,f3) = C,JdpHiUxx')S{x-Tr{Trx'mxx'{'^\?) 

(7.14) 

where in computing ttaIp) we have used the fact that the reduced density matrix 
of subsystem A is obtained from \'^)xx' W tracing over X' and B, i.e. 

PA = Trjsip) = TrjsiTrx' |*)xx' (7-15) 

Analogously to what seen for the pure case, a; = 1, by writing p = Ux^xUx with 
Ax = diag(Ai, A2, Al), we finally get 

Za{x,P)^C,J dpniUx) J da{Ax) s(x-TrAl) ^-^^^^Biu^K^u.^)? ^ 

(7.16) 

where dpniUx) is the Haar measure on U{L) and 

da{Kx)^CL5(l- ^) n ^(^^) n (Ai - A.O^d^A, (7.17) 

l<i<L l<i<L l<i<j<L 
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with ^ 

° r(L + i)n,<.<,W 

Notice that, from the properties of the Euler gamma function F, the normahzation 
constant Cl concides with Cn,m in Eq. (2.6) if we replace and M with L, namely 
Cl is equivalent to Cn,m in Eq. (7.8) for iV = M = L. The above definition is 
equivalent to have identified the measure dn^ip) of Eq. (7.9) with 

d^ia.(p) = C,diiH(Ux)da{Ax)5(^x - IVA^) . (7.19) 

In the case of pure states, i.e. x = 1, when the density operator of the system 
reduces to p = the matrix Ax becomes a rank one projection. Thus 

expression (7.16) reduces to (7.7): 

Za{1,I3) = Zab{I3) . (7.20) 

7.3 Asymptotic behavior and analysis of moments 

In the first part of this thesis we have apphed the method of the steepest descent in 
order to determine the leading contribution to the partition function at x — 1. This 

analysis could be performed because of the factorization of the integration over the 
unitary group. It led to a complete characterization (in the thermodynamic limit) 
of the statistical distribution of the local purity of X, tiabW = ''^aH') = 7Tb{'4') 
(and hence of the bipartite entanglement of the system) computed with respect 
to a canonical ensemble for different regions of the Hilbert space, selected by the 
inverse temperature /5. In other words, in the limit of large N, we have determined 
the most probable distribution of the Schmidt coefficients which maximizes the 
distribution (3.64) 

dpp{pA) = dp{pA) — 7^ , (7.21) 

i.e. typical states with respect to the above canonical measure. In this context we 
have noticed that the system undergoes two main phase transitions (for the sta- 
ble solution), related to difi^crcnt distributions of the eigenvalues A^ of the typical 
states: a second order phase transition, mentioned above, associated to a Z2 sym- 
metry breaking, and related to the vanishing of some eigenvalues of p^, followed by 
a first order phase transition, associated to the evaporation of the largest eigenvalue 
from the sea of the others. By tuning the inverse temperature (or equivalently by 
deforming the Haar measure with the non uniform weight e~^'^^^^''^y Zab{I^))i we 
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have determined the statistical behavior of the system at different levels of entan- 
glement. As noticed in Sec. 3.3.3, this is equivalent to determine all the moments 
of the purity at different values of /3. 

For X ^ 1, the integration over the unitary group U{L) in Eq. (7.16) does not 
factorize, making the computation of the partition function far more complicated 
than for the case of pure states (7.7). The only notable exception is the case of 
maximally mixed states (i.e. x — 1/L), when the Dirac delta of Eq. (7.16) selects 
a unique diagonal matrix Ax (the totally mixed state of X). This makes the 
exponent equal to e~^/^ for all Ux and yields the following exact expression 

Z4l/L,/3)=e-^/^ (7.22) 

However, even if for x ^ 1, 1/-^ the situation is much more complicated, for small 
/5, the evaluation of the moments M.^{x, (3) associated with Za{x, /3) admits exact 
analytical treatment, as we will show in the next two chapters. The moments of 
the purity for a generic set of states Sxi'K-x) are formally defined as 

M^ix, 13) = J df^M ^A-{P) = ^^^^^ , (7-23) 

and represent the average value of t:a^{p) over the canonical measure 

dpxAp) = dnx{p) ^7 — ^ , (7.24) 

with dfij.{p) as in Eq. (7.19). For pure states {x = 1) they coincide with the 
moments A^„(/3) defined in Eq. (3.63). In the case of a totally mixed state {x = 
1/L) Eq. (7.22) yields instead values which are independent from the temperature 
/3, namely 

M^{l/L, /3) = Mt{l/L, 0) = N-\ (7.25) 
For intermediate values of x and for sufficiently small /? we can write 

M^{x, /3) ~ Miix, 0) -(3 [Mtii^, 0)^Mt{x, 0) Mf,{x, 0)] , (3^ 0, (7.26) 

which has been obtained by expanding Eq. (7.23) up to the first order contribution 
in (3. Incidentally, notice that in agreement with Eq. (7.25), the /3-corrections of 
Eq. (7.26) vanish when x = 1/L. The above expression shows that, at least in the 
high temperature regime, we can focus on the unbiased moments A4^{x, 0). In the 
remaining part of the thesis we will compute the moments of the purity and infer, 
at the end of Chap. 9, the asymptotic behavior of the average local purity in the 
mean field approximation. 
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In this chapter we will focus on the first moment of the local purity 7r_4(V^) for a 
generic state p G &x{T^x)- In particular, we will consider two distinct approaches 
based on the Gaussian approximation for the elements of the unitary matrices 
belonging to hl{L) and U{L'^), where L is the dimension of the global system. 

In Sec. 8.1 we will use the fact that the integral over the eigenvalues of the 
density operator p, with purity x, and the integral over the unitary group associ- 
ated to its eigenvectors can be evaluated independently, the former by introducing 
the empirical density distribution of the eigenvalues, the latter by approximating 
the elements of the unitary matrices acting on T-ix as independent and identical 
complex Gaussian random variables. In Sec. 8.2, we will follow a different strat- 
egy which applies the same Gaussian approximation to the enlarged Hilbcrt space 
Hxx'- Both approaches will lead to the evaluation of the first moment of the local 
purity at the order C (1/L). 

8.1 Divide and conquer 

The approach we will discuss in this section for the computation of the first moment 
M.i{x, 0) of the local purity TTyi (■?/') is based on the possibility to perform separately 
the integral over the unitary group and the integral over the spectrum. Indeed, let 
us consider the expression (7.16) of the partition function. It is worth introducing 
the spectral decomposition of Ax, i.e. 



where Pi = \i)x{'i-\ with 1 < i < L are the projections over the reference basis 
C Hx, that is Pj = {Sk,A,i), I < k,l < L. By defining the reduced density 




(8.1) 



l<i<L 
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matrix p^i = PAiiUx) = TrB(^7x-Pi^x^), associated to each eigenstate Ux \i)x 
p, the partition function becomes 

ZA{x,P)^C,j dpH{Ux) j da{Ax)s(x- A2)e-^^--<-^'^^'^-(^-) , (8.2) 



l<i<L 



where HijiUx) = Tryi(p^jp^j). As aheady observed the exphcit computation of 
Za{x^ 0) of Eq. (8.2) involves a nontrivial integration over the imitary group (Itzyk- 
son and Zuber [1980]). We will consider here the first order moment (7.23). We 
get 



da{Ax)s(^x- J2 ^l)^i^J dpH{Ux)n, 



= E 

i<i,j<L 

= J2 {\Xj){7rij{Ux)) , 



l<k<L 



^ij{Ux 



(8.3) 



l<i,j<L 



where {XiXj) and {"KijiUx)) stands for two independent averages (the first over 
the eigenvalues and the second over the unitary matrices Ux)- We will separately 
compute these quantities in Sec. 8.1.1 and 8.1.2, respectively, and then compute 
the first order moment of the purity A4^(a;, 0) up to order O (1/L). 



8.1.1 Average over the eigenvalues 

In order to evaluate the average over the eigenvalues (AjAj), 1 < i, j < L, we will in- 
troduce the empirical density function of the eigenvalues p{\) and the joint density 
function p(A, A'). Given a typical realization of eigenvalues, the empirical density 
function is defined as 

P^^^^l E Ae[0,l]. (8.4) 

l<i<L 

The above definition enables to encode the information on the eigenvalues Aj, 
1 < i < L, in the probability distribution p{\) on the interval [0, 1]. The empirical 
density function is thus correctly normalized 

[ d\p{\) = 1. (8.5) 
Jo 

We also have 

(A«)^ fd\X^p{\) = \ (8-6) 
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where, in fact, in order to compute (A") we should take the average of (8.4) over 
the distribution of the eigenvalues (in particular we get (A) = ^/L). We will show 
that the contribution of fluctuations is of higher order in 1/L. If the global system 
is in a generic mixed state, with purity Trp^ = x e [1/L, 1], the distribution of the 
eigenvalues must satisfy the constraint 

E A,2 = X. (8.7) 

l<i<L 

Since the purity is of O (1), the above condition can be satisfied only if one or more 
eigenvalues are of (9 (1). In what follows we will first consider the simplest scenario 
with only one dominating eigenvalue, and then we will generalize the analysis to 
the case of more than one eigenvalue of order C (1), proving the equivalence of the 
two approaches. 

One eigenvalue of O (1) 

Let assume to have one eigenvalue of O (1), detached from the sea of the first L — 1 
eigenvalues of O (1/L), namely 

Al = v^, Y1 Ai = l-v^, (8.8) 

l<i<L-l 

with Xi — 0{1/L) for i < L. The constraint (8.7) is satisfied up to order O (1/L) 



1 T. N / 



(8.9) 



l<i<L 

In this hypothesis, the distribution of the eigenvalues (8.4) splits as 

P(A) = lH^-^) + (^^-l) Pse^W, (8-10) 

where Psea(A) is the normalized probability density of the first L — 1 eigenvalues 

1 /"^ 

Psea(A) = -^ J2 '^(A-AO, / ciApsea(A) = 1. (8.11) 

l<i<L-l 

Its first two moments are given by 

(A)3ea = jy^^Pse.{X)=^-J^, (8-12) 
•^0 l<i<L-l ^ ^ 
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As before, (A^)ggg^ has been computed by neglecting the average over the distribution 
of the eigenvalues. This choice will be justified at the end of this subsection. In 
order to evaluate the average (XiXj), 1 < i, j < L we need to introduce the joint 
empirical density function, which similarly to Eq. (8.4) can be written as 

p{X, A') = — ^ Y.' '^(^ - '^(^ - ^^O' (8-14) 

^ ' l<i,j<L 

the prime standing for i ^ j. Under the above hypothesis the joint density proba- 
bility splits as 

p(A, A') = -^5(A-VS)psea(A') + ^5(A'- VS)Psea(A)+ (^1 - |^ Psea(A, A'), (8.15) 

where we have introduced the joint probability density of the sea: 

Psea(A, AO ^ ^1 _ (5(A-A,)5(A'-A,). (8.16) 

It easy to show that the above definition can be written as 

Psea(A,AO = p,,,(A)psea(A')--^Psea(A)[5(A-A')-Psea(A')]. (8.17) 

Notice that Psea(A) is the marginal probability of the above joint density function, 
that is ^ 

Psea(A)= / (iA'psea(A,A'). (8.18) 
JO 

Furthermore, given A and A' belonging to the sea of the first L — 1 eigenvalues, we 
have 

(AA').ea = (A)L-^((A^La-WL) (8-19) 

and from Eqs. (8.12) and (8.13) we get 

We are now ready to compute the explicit expression for (A^) and (AA') (with 
A^A'): 

i2\ I / 1 ^ \ /\2\ 



(AA') = |x/^(A)sea+ (AA')sea. (8.22) 
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In the limit of large values of L, we get 

X 



{^') = y + ^ly^)' (8-23) 



L \L 

Notice that the average values computed above show fluctuations of O (l/L"^) and 
0{1/L^'). We will soon show that this result docs not depend on our hypothesis 
on the distribution of the eigenvalues, and can be generalized to the case of more 
eigenvalues of O {!). 

Before considering this point, let us compute the density distribution of the sea 
of eigenvalues in the thermodynamic limit. Similarly to what seen in Sec. 3.1 we 
can introduce the natural scaling for the sea of eigenvalues 

A, = (l-V^)-^, 0<t,= -^<l, Vz = l,...,L-l (8.25) 
and define, in the thermodynamic limit, the density function 

Psea(A) = ^ dtS(X- A(i)) . (8.26) 
The distribution Psea(A) is normalized and satisfies the constraint (8.8): 

(iApsea(A) = 1 (8.27) 

dAAp,ea(A) = (8.28) 



being [0, a] the domain of the scaled eigenvalues and a > 0. In analogy to the 
Wishart distribution found at ^ = for the eigenvalues pA in Chap. 3, we set 



PsUX) = Cxl^, (8.29) 



where a and C are determined by conditions (8.27) and (8.28): 

aC 



2^ = 1 (8.30) 
= 1-^/x. (8.31) 
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Therefore, expression (8.29) becomes: 



^-°-t^) = Mr^ f ""f ^^P-^d-vJ)] (8.32) 

being a = 4(1 — ^/x) and C = l/27r(l — y/x). Notice that when x = 0, Psea 
reduces to the unbiased average density at (3 ~ found for the distribution of the 
eigenvalues of the reduce density matrix in Chap. 3. By comparing (8.11) and 
(8.26) we get 

L-1 ~, L-1 / L- 1 \ 

Psea(A) ^ Psea(A) = ^ Psea A ^ , (8.33) 



from which all moments can be easily evaluated 

(nea-/.AAWA)^(^)/ 

(1 - v^)'" 4-r [n + 1) 
(L-1)" ^{n + l)l' 



dyy^'pses.iy) 

(8.34) 



In particular, (A^)^^^ = 2 (1 - y^) /{L - if, in agreement with (8.13), that is 

(A^)sea = 0(lA^)- 



k eigenvalues of O (1) 

In general, for large L, a density matrix p with purity x can have finite number k 
of eigenvalues of order O (1), say Aj = y/xl ior L — k < i < L, with 



= a;, (8.35) 



L-k<i<L 



and the first L — A; eigenvalues Aj = 0(l/L)forl<i<L — /c. Thus the constraint 
on the purity of the global system is satisfied up to order (9 (1/L), that is 

+ (8-36) 



l<i<L 

The empirical density (8.10) generahzes to 



P(^)-l E <^(A-V^)+(l-|)psea(A), (8.37) 
L-k<i<L ^ ^ 
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with the empirical density for the sea of the first L — k eigenvalues of C (1/L) given 

by 

Psea(A) = -^ J2 S{X-Xi). (8.38) 
l<i<L-k 

Thus in this set up the joint probability density shows three main contributions 



p(A, A') 



L(L-l) 
L-k 



L—k<i,j<L 



_ -^x 5^ (A - \/x'i) Psea(A') + 5 (A' - ^) Psea(A)] 

Psea(A,A'), (8.39) 



L-k<i<L 

{L-k){L-k-l) 



L(L-l) 

being Psea(A, A') the joint probability density for the sea, analog to (8.16), 



Psea(A, A') 



{L-k){L 



i-^-- ^'5(A-A,)5(A'-A,). 



(8.40) 



l<i,j<L-k 



Wc can now compute the average product of eigenvalues (AjAj), 1 < i, j < L. For 
Aj — Xj, from Eq. (8.37) we have 



(A^ 



L-k<i<L ^ 



l<i<L-k 

On the other hand for Aj 7^ A^, from Eq. (8.39) we have 
1 2 



.41) 



(AA') 



L(L-l) ^ ^ ' ' L(L-l) 



I). 5: 



Xi 1 



L-k<j<L 



{L-k- 1] 



L{L-l){L-k) 
1 



L(L-l) 



1 - E 

Z,-fc<i<L 

+2 E 



L—k<i,j<L 



L-k<i<L 
2 ■ 



1- E 

L-k<j<L 



E 



+ I 1 - > VX,: 
L-fc<i<L 



L(L-l) 



(8.42) 



L-k<i<L 
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that recalling condition (8.35) becomes: 

(AA'> = (8.43) 

By comparing Eqs. (8.23) and (8.24) with Eqs. (8.41) and (8.43), we conclude that 
the average over the eigenvalues does not depend on our assumption on the (finite) 
number of eigenvalues of O (1). 

8.1.2 Unitary group average 

The next step, in order to compute the first moment of the purity (8.3), is the 
evaluation of the average over the unitary group 

{TTijiUx)) = J dfiH{Ux)7^ij{Ux), (8.44) 

where nij{Ux) = TY{pAiPAj), and pAk = ^t^b{Ux \k)x {k\ Ux^), \k)x being the k-th 
vector of the reference basis, with 1 < i, j, k < L. 

Let us, first of all, choose a convenient representation of the tensor product 
structure of Hx = 'Ha®'Hb: i.e. = ® C^, with L = NM. Given ua = 
{un)i<n<N £ and vb — {vm)i<m<M G C^, wc define the tensor product of the 
two vectors as 

wx = ua^vb = {'We)i<e<L e (8.45) 

where 

'e = M(n - 1) +m, 
Wi = UnVm, {n^ [^] + 1, (8.46) 
m = £-M(n-l), 

[x] being the integer part of x. Therefore, for any w € C^, the representation of 
the partial trace 

PAiw) = Ttb {\w)x{w\) = iPA,ij)l<i,j<N, (8.47) 

reads 

{Pa{w))i^^ = WM{i-i)+k WM{m-i)+k , V Z, m G {1, . . . , N}. (8.48) 

l<k<M 
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In order to write the explicit expression for T^ijiUx), we have to choose the com- 
putational orthonormal basis for the system. Let U — {uki)\<k,i<L be an element 
of the unitary group U{L). Its columns u,i = {uu)i<k<Li with 1 < / < form 
an orthonormal basis of C^, i.e. = Ylii<k<L'^kiUkm — ^i,m, as follows from 

WU — 1. Thus, we can write 

TTijiU) = IV^ (paMpaM) , V i,i e {1, . . . , L} (8.49) 

where the elements of the reduced density matrix ^^('W.i), 1 < i < iV, are given by 

{PA{u,i))mn = X] UM(m-l)+k,i UM(n-l)+k,i V m, 77, G {1, . . . , N} . (8.50) 
l<k<M 

From the last two relations we can finally write the explicit expression for i^ijiUx) 

= 5^ 5^ UM{m-l)+k,i UM{n-l)+k,iUM{n-l)+l,j UM(m-l)+l,j , 

l<k,l<M l<m,n<N 

Vi,je{l,...7V}. (8.51) 

Let us compute the average (8.44) over the unitary group. If G U{L) is sam- 
pled according to the Haar measure, we can parametrize its columns as u,j — 
{rke^'^'')i<k<L for 1 < j < L: the phases (f)^ are independent identically distributed 
(i.i.d.) random variables with a uniform probability density on [— vr, vr], and the 
squared moduli r| arc uniformly distributed on the simplex = {(r|)i<fc<i | < 
rl < l,X]^fc — 1} (Zyczkowski and Sommers [2007]). Since the volume of the 
simplex is [A^,! = 1/(1/ — 1)!, the joint probability density of the square moduli, on 



i^, reads 



and therefore 



PL{rlrl...,rl) = {L-l)\6(l- ^ , (8.52) 

V l<fe<L / 



PL(ri,r2,...,ri)=2^(L-l)!<5(l- r^J J] r^. (8.53) 

In appendix B we show how L random variables can be generated with a uniform 
distribution on the simplex A^. 

If we integrate out the last L — K variables, 1 < X < L, by defining 

Vi^rl^i, yie{l,...,L-K} (8.54) 
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we find 



PLin,r2,...,rK) = 2^(L-1)! J] r, /' d^'^y <5(l - ^ r,^ - ^ y,) 

l<j<i<: "'0 l<i<K l<j<L-K 

nl r^-yi r^-yi-y2 

= 2^(L-1)! TT rj / dy, / d?/2 / dya ■ ■ ■ 

T^^^i^ Jo Jo Jo 



l<j<K 
-l-yi--yL-K-3 

X ' 



r^—yi-.—yL-K-i I 

j dyL-K-2 1 1 - XI ~ Yl 



l<i<K l<j<L-K-2 

(8.55) 



and finally 

Pdn, r,,..., tk) = 2^7^3^31)1 ( 1 - E ) n (8-56) 



L-_ft:-i 



l<i</i" / l<j<K 



In particular for K = 1 we get the probability density of the amplitude of an 
arbitrary element of U 

PL{r) = 2{L-l)r{l-r^)^-^. (8.57) 
The n-th moments of this density function are given by: 

(r") = {L- 1)B (- + 1, L - l) = (L - ^y-Y^J^^ (8-58) 

where B{x,y) is the Beta function, 

B{x, y) = dt t^-\l - t)y-' = Y^^- (8.59) 

For large L (L ^ oo) the probability density piir) becomes 

PL{r) ~ 2Lre-^''\ (8.60) 

It follows that in the thermodynamic limit we can approximate Ux = (i^ij) with 
Uij = rjj-e*'^'-', i.i.d. random variables with probability density functions 

p{nj) = 2Ln,e-^4 , p(0, .) = i-, (8.61) 

that is ^ 

p{uij) = -exp{-L\uij\'^), Vi,j e {1,...L}. (8.62) 

TT 
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We then have that Ui/s are i.i.d. complex Gaussian random variables with zero 
mean and variance = 1/L. Notice that, in the above approximation, uni- 

tarity is satisfied at order 1/VL. Indeed, 

^ UklUkm ) = (UklUkm) = ^ ^ = ^ (8.63) 

\l<k<L / l<k<L l<k<L 



while 

UklUkm — ^l,m 

l<k<L 



2\ ^ 

/ ^ XI {^klUkmUjlUjm)-5i^rn = V Z, m G {1, . . . L}. 

/ l<k,j<L 

(8.64) 

The above quantity measures the additional fluctuations introduced by the Gaus- 
sian approximation, which enables us to easily evaluate the average 

("""ij) = X {uM{m-l)+k,i UM{n-l)+k,iUM(n-l)+l,j UM{m-l)+l,j) , 

l<k,l<M l<m,n<N 

Vi,je{l,...L}. (8.65) 

Indeed, by {\uij\'^)x — 1/L and {\uij\'^)x — 2/L^ it follows that, when 1 < i j < 



L, 



{uM{m-l)+k,i UM(n-l)+k,iUM{n-l)+l,j UM{m-l)+l,j) 

1? 



= 7f (8.66) 



while, for 1 < i = J < I/, 

{uM{m-l)+k,i UM(n-l)+k,iUM{n-l)+l,i UM{m-l)+l,i) 

= iSm,n + Sk,l) {\uM[rn-l)+k,if) {\uM{n-l)+l,if) — Sm,nSk,l {\uM{m-l)+k,i\^) 
_ ^m,n + Sk,l — '2-5m,n5k,l 

Vm,n e {1,...A^} and V A;, Z e {1, . . . M}. 

(8.67) 

Incidentally, notice that (8.67) can be computed exactly by using (8.58) and reads 
{5m,n + 5k,i)/L'^ - 25m,Jk,i/L{L + 1). Therefore, 

, , , NM^ + N'^M -2NM , . NM'^ 
K^iAUx)) = Sij — + (1-^^^)-^- 

When N = M = VL, one gets (iry) = {H{)>a,PAj)) ~ (1 + ki)/^- 
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8.1.3 Average purity 

From the eigenvalue averages (8.23) and (8.24), and from the average (8.68) over 
the unitary group we can evaluate the average purity at order O (1/L) 



1 T 1 T \ 



M 



i<i,j<L i<i,j<L 

'i-<i,j<L ^ ^ l<i<L 

1 , , f 1 1\ 1 X , , 



which ioT N — M — \fL yields 

Mi^(x,0)~^(l + x). (8.70) 

Notice that for a; = 1 Eq. (8.69) is in agreement with the thermodynamic behavior 
of the pure states distribution computed in (Scott and Caves [2003]; Facchi, Mar- 
zohno, Parisi, Pascazio and Scardicchio [2008]) for the general case of unbalanced 
bipartions. Furthermore, Eq. (8.70) yields A1^(l,0) = {t^ab) = 2/iV, according to 
the results of Chap. 3 for N = M (see Fig. 3.8 where u — N {ttab}), which can be 
also inferred from Eq. (3.53) for n = 1. 



8.2 Multiply and conquer 

In this section we introduce an alternative approach for the computation of the 
cumulants of the partition function Za{x,/3) of Eq. (8.2) that will allow us to de- 
termine the leading terms of the average purity Aif{x, 0). Wc will then compute 
the first moment of the distribution of the local purity and compare the result 
with the one obtained in the previous section. In particular, we go back to expres- 
sion (7.14) of the partition function which involves the integration over the larger 
unitary group U{L^) acting on the enlarged purification space T-Lxx' — 'Hx <8) Tix' 

ZA{x,f3) = C,JdfiH{Uxx')H^-T^iT^x'mxx'W) 

xe-l^ i^(^sCiVx'(l*>xx'(*l)))'. (8.71) 

Here, \'^)xx' ~ ^xx' \^o)xx'^ where |^o)xx' ^ reference state and Uxx' an 
element of the unitary group U^L"^). If we choose the vector (5i,i)i<i<L2 € 
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as reference state, the partition function (8.71), according to the tensor product 
representation introduced in Sec. 8.1.2, reads 

Za{x,P) ^C,J dl,H{u.l) S {X - TT{TTx'\u.l)xX'{u.l\f) e-/^^(^-Tr,,|n.i>,,,(n.,|)^^ 

(8.72) 

where 

1 2 2 

dnH{u,i) = ,^ sr2 PL^{ri,r2, ■ ■ ■■,rL2) (p r (8.73) 
(27rj^ 

is the exphcit expression of the Haar measure over the unitary group for the repre- 
sentation {uii) = (rjC**^'), 1 < i < L^, with the joint probabihty density Pl2 given 
by Eq. (8.53). By introducing the simphfied notation, Zi = Un, we rewrite the 
partition in the more compact form 

Za{x, P)^ J dfxniz) S {x - nx{z)) e-^-^(^), (8.74) 

and the Haar measure over the unitary group reads 

df,H{z)=6(i- \zA n (8-75) 

y i<fe<L2 J i<j<L2 

With the same technique introduced in the first part of this chapter, let us de- 
termine the exphcit expression for the purity for 7Tx{z) and nA{z). The matrix 
elements of the density matrix p according to the notation introduced in Sec. 8.1.2 
are 

Pmn = {T^X'\u,i){u,i\)mn = ZL(m-l)+k ZLin-l)+k, W 171, U E {1, . . . , L} . 

l<k<L 

(8.76) 

We then get that the purity of the global system, whose value is constrained to be 
X, is 

T^x{z) ^ ^ ZL{m-l)+k ZL{n-l)+k ZL{n-l)+l ZL{m-l)+l ^ X, (8.77) 

l<k,l,m,n<L 

whereas for the purity of subsystem A we get 

T^Aiz) = ^ ^ ^ ZL(M{i-l)+l-l)+k ZL{Mij-l)+l-l)+k 
l<k,k'<L l<i,j<N 1<1,1'<M 

X ZL{M{j-l)+l'-l)+k' ZL{M{i-l)+l'-l)+k'- 

(8.78) 
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Analogously to what seen in Sec. 8.1.2, for large L we will consider the ZiS as i.i.d. 
complex Gaussian random variables, i.e. 

d/iniz) ^ TT — ex.p{-L'^\zi\'^)dzidzi. (8.79) 

J. J. 7^ 

l<i<Z,2 

The fluctuations due to this approximation scale as 0(1/L^) (see Eqs. (8.62) - 
(8.64)). However, will show that the fluctuations of the average purity are of order 
0{1/L). We now have all the elements for obtaining an alternative - and more 
direct - derivation of (8.69). 



8.2.1 Average purity 

The average purity (8.3) reads 

Mtix, 0) = y" d^iHiz) 5 (x - nxiz)) 7r^(z), (8.80) 
and involves 4-point correlation functions the form 

{zL{M{i-l)+l-l)+k ZL{M{j-l)+l-l)+k ZL{M{j-l)+l'-l)+k' ^L(M(i-l)+Z'-l)+fe') , (8.81) 

with 1 < i,j < N, 1 < k,k' < L, and 1 < 1,1' < M. See Eq. (8.78). The 
two dominant contributions come from the contractions 5ij and 5i^i' that involve 
sums over or M terms in (8.78). The remaining contractions yield higher order 
corrections. The first contraction reads 

{\zL{M(i-i)+i-i)+k\^) {\zL(M{i-i)+i'-i)+k'\^) = -j^, ^ i,j ^ {1, - N} (8.82) 



hence the first leading term is given by 

(''.^>.= E E E 

l<k,k'<L l<i,j<N 1<1,1'<M 



L4 



1 

TV' 



The second contraction reads 

^1,1' {zL{m-l)+k ZL{n-l)+k ZL{n-l)+k' ^L(m-l)+fc') , V /, /' G {1, . . . M} 

with 

m = M{i-l) + l, n = M{j - 1) + I' . 
Note that the condition I — I' is equivalent to thb — tib, where 



niB = m — M 



m 



M 



ub = n — M 



n — 1 
M 



(8.83) 

(8.84) 
(8.85) 

(8.86) 
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Thus the other dominant contribution reads 



M2 



^mB,nB {^L{m-l)+k ^L(n-l)+fe ^L(n-l)+fe' ^L(m-l)+fe')x ' (8.87) 

l<k,k' ,m,n<L 



which, by virtue of (8.77), yields 



X 



88) 



M 



On the other hand, if we require both the conditions 5ij and 5;^// for 1 < z,,; < A?" 
and 1 < Z, < M to be satisfied, we get: 



Other correlations between the above indices lead to higher order corrections. Sum- 
ming up the results we obtain 



in agreement with (8.69), and with the first cumulant computed in (Facchi, Mar- 
zohno, Parisi, Pascazio and Scardicchio [2008]) for both balanced and unbalanced 
bipartitions, at x — 1. 

The techniques discussed in this chapter can be generalized to the computation 
of the leading terms of the higher order moments of the local purity. In the next 
chapter we will exactly compute the first two moments of the local purity for a 
generic set of states in &x{'Hx), confirming the results found in the previous section 
thru the Gaussian approximation. 




(8.89) 




(8.90) 



Chapter 9 

Exact moments and twirling maps 



In this chapter we will determine the high temperature expansion of the first mo- 
ment of the local purity, from the exact computation of its first two moments for 
/3 = 0. This generahzes previous results achieved for the case of pure states both 
in the thermodynamic limit (Lubkin [1978]; Lloyd and Pagels [1988]; Zyczkowski 
and Sommers [2007]; Facchi, Marzolino, Parisi, Pascazio and Scardicchio [2008]) 
and for systems with finite dimension (Scott and Caves [2003]; Giraud [2007a]). 
In Sec. 9.1 we will establish a formal connection between our problem and the 
theory of quantum channels. We will exploit the symmetry properties of the so 
called twirling transformations in order to compute the exact expression of the 
first moment of the purity of one subsystem for the case of pure states and then, 
through a purification scheme apply the same strategy to the case of arbitrary 
mixed states. However the generalization of this procedure to the computation of 
higher moments is very hard to be performed, since the purification procedure in- 
troduces a higher number of copies of the system making the structure of the total 
Hilbert space more involved. In Sec. 9.2 we will introduce an alternative approach 
based on the solution of some basic integrals over the unitary group. We will then 
compute the exact expression of the first two cumulants, and determine in Sec. 9.3 
the high temperature expansion of the first moment of the local purity. 



9.1 Exact computation of the first moment by a 
twirling map 

In this section we will determine the exact expression of the first moment of the 
local purity of an arbitrary mixed state. In particular we will exploit the proper- 
ties of the so called twirling transformations (VoUebrech and Werner [2001]; Lee, 
Chi, Oh and Kim [2003]; Hayashi, Markham, Murao, Owari and Virmani [2008]; 
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Bennett, Di Vincenzo, Smolin and Wootters [1996]; Bennett, Brassard, Popescu, 
Schumacher, Smohn and Wootters [1996]). This approach will also establish a 
formal connection between our problem and the theory of quantum channels. We 
start in Sec. 9.1.1 by deriving a general expression for the moments A^^(1,0) in 
the case of pure p (i.e. the quantities Mn of Eq. (3.52)) and verify that it yields the 
exact value given by Page (Page [1993]; Giraud [2007a]; Scott and Caves [2003]) for 
n — 1. We will then generahze this technique to the case of mixed p in Sec. 9.1.2, 
and compute the exact expression of the first moment of the local purity. 



9.1.1 Pure initial states 

Let us consider a fixed (normalized) pure state {'4') ab the global system X — AB 
and parametrize the pure states of T-Lx through the mapping \iPu)ab ~ ^ab 
with Uab e U{'Hx) — U{L) distributed according to the unique invariant Haar 
measure dfiH{U), as done in Chap. 2. The reduced density matrix of subsystem A 
associated to \iPu)ab given by 

PA = TrB[UAB\ll^)AB{ll^\Uis]^ Yl BmAB\^)AB{^\Uis\£)B (9.1) 

i<e<M 



where {|^)b}i<£<m is an orthonormal basis of T-Lb, with dim'Hs = M. Prom Eq. 
(9.1) the local purity of |'0?7)ab is given by 



T^AbW = Yl T^^[BmAB\i^)AB{i^\UU^)BB{i'\UAB\^)AB{^\UiB\nB 

i<e,e'<M 

= J2 AB{i^\{uUi)B{i'\UAB)\^P) AB AB 



i<e,e'<M 



(9.2) 



where in the last term we used the property of the trace Trdip) {ipl O) — 

operator O on T-Lx, and write the resulting expression as a 
product of two expectation values on \iP)ab- We can recast Eq. (9.2) into a more 
compact form by doubling the Hilbert space, i.e. adding the auxiliary copies A' 
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that is 



and B' of A and B, respectively. Namely, 

i<e,i'<M 

XA'B'{lp\ (uls,\i')B'{i\UA'B') Wa'B' 

i<ei'<M 



X 



(Uab C/a' 



B' 



lA'B' 



IV') as), 



(9.3) 



7rAB(V') ^ Tr (UAB^UA'B')(^\ij)AB{ij\^\ij)A'B'{'l/j\) (U^B^UI,^,) 



X (5b|B' ® Iaa') 



(9.4) 



Here, the trace is performed over all degrees of freedom (i.e. AA'BB'), Iaa' is the 
identity operator on AA , and 



^i<e,e'<M 



(9.5) 



is the swap operator on BB'. The swap operator is a unitary, self-adjoint trans- 
formation which, for all operators Qb and T^/ on T-Lb and T-Lb', gives 



SBlB'i^B ® '^B')Sb\B' = (g) ©B/ 



(9.6) 



Recall that by definition of the Haar measure on U{L), the first moment Aii{x = 
1, 0) = M.1 of Eq. (3.52) is obtained by averaging over all possible Ux G W(L) (see 
also Sec. 2.1), that is 



Ml 



Tr 



r(2) {\iP)x{ip\ ® \ip)x'{ip\) [Sb\b' (8) Iaa') 



(9.7) 



where X = AB, X' = A'B', and T*-^-* is the Completely Positive Trace Preserving 
(CPTP) twirling channel (Vollebrech and Werner [2001]; Lee, Chi, Oh and Kim 
[2003]; Hayashi, Markham, Murao, Owari and Virmani [2008]; Bennett, Di Vin- 
cenzo, Mor, Shor, Smohn and Terhal [1999]) which transforms the operators Qxx' 
of XX' into 



r'-'\Qxx') = J diiH{U) {Ux ® Ux') Qxx' {U^x ® ^1') , 



(9.8) 
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duniU) being the usual Haar measure. This map plays an important role in quan- 
tum information theory where it was first introduced as a tool for characterizing the 
distillability of bipartite entanglement (Bennett, Brassard, Popescu, Schumacher, 
Smolin and Wootters [1996]; Bennett, Di Vincenzo, Smolin and Wootters [1996]). 
It has several properties which allow us to simplify the calculation. For instance, 
it is known that T^^^ maps all the states of the system into (generalized) Werner 
states (VoUebrech and Werner [2001]; Lee, Chi, Oh and Kim [2003]). Furthermore 
it is self-adjoint - i.e. its description in Heisenberg picture, when we extend the 
action of quantum channels from density matrices to linear operators on Hxx', 
coincides with T*^^^. In particular, this last property can be used to rewrite (9.7) 
as 

Ml = (^x{^\(^x'{^\)T^^^iSB\B'(^lAA') {\^)x(^\^)x') . (9.9) 

One of the main reasons that make the twirling channels very useful for our scopes 
is that the explicit expression for the action of T*-^-* can be obtained by exploiting 
the symmetry of duniU). In particular it is possible to show that T'^'^\Qxx') can 
be decomposed as a linear combination of projectors on the symmetric and antisym- 
metric subspaces of XX' — ABA'B' (with respect to the bipartition AB\A'B'). 
More precisely, introducing the swap operator which exchanges X with X' (or 
equivalently A with A' and B with B') 

Sx\X' — Sab\A'B' — Sb\B' <H) Sa\A': (9.10) 

it can be shown that the action of T^^^ is given by 

r(^)(e..o = ^^g^^Tr(e..o + ^^f^iv(^x,x^e..o 

LTr(exx') - Tr{Sx\x'exx') , 



L(L2 - 1) 



I XX' 



LTr{Sx\x'Qxx') - T^j^xx') ^ .„^^. 

+ Kl^) "^^1^' ' ^^-^^^ 

where L = NM is the dimension of Hx and Ixx' is the identity operator on Hxx'- 
This is the central formula of our approach with the twirling transformations, and 
will lead to the exact computation of the first moment of the local purity both for 
pure and mixed global states. This can be achieved for x = 1, either using Eq. (9.7) 
or Eq. (9.9). We have first to compute the quantities Tr(Gxx') and Tt{Sx\x'Qxx') 
with Qxx' = \'ip)x{'ip\ ® \ip)x'{'4'\- 

Tr(|V')x(V'l®lV')x'(^|) = l, 

TtiiSBlB' SAlA'MABilpl \lp)A'B'{lp\] = 1 , (9.12) 
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where in the second expression we have used the fact that \iI')ab 
invariant uder -S'xix', i-e. 

{Sb\B' ® Sa\A'){\'4')ab ® \ij)A'B') = \i^)AB ® |V')a'S'- 

Replacing all this in Eq. (9.11) we get 



r^'\\ip)x{ip\^\ip)x'm 



Ixx' + Sx\x' 
L{L + 1) ' 



and thus 



Ml 



1 



-Tr 



L(L + 1^ 
1 

N{N + 1) 



{Iaba'B' + Sb\b' ® 'S'^iA') {Sb\b' ® Iaai) j 



Tr(5'B|B' /aa') + Tr(/BB/ (g) 5'a|a') 



(Here we exploited the fact that S'^^^, = Ibb')- Notice that 

T^{Sb\B') = BW)B{i'\^BB'm')B'mB' = M, 

l<i,j,l,e<M 

while 

Tr (7bbO = = dim(?{B 7^^/). 
Thus, by using the condition L — NM and the identities 



Tx{Sb\b' ® Iaa') = MN^ , Tr{lBB' ® Sa\a') = NM' 



we get 



Ml 



N + M 
NM + 1 



)A'B' IS 



(9.13) 



(9.14) 



(9.15) 

(9.16) 

(9.17) 
(9.18) 

(9.19) 



which coincides with the exact value found in (Giraud [2007a]), see Eq. (3.58). 

We mention that the same technique can also be applied to higher moments 
Mn- The extension of Eq. (9.7) for n ^ 2 is obtained by introducing 2n copies of 
AB organized in n couples, i.e. AiBi,A[B[, ^2^2, ^2^2' ^nBn,A'^B'^. We 
then introduce the following generalized twirling transformation acting on XX' — 
AiBiA[B[A2B2A',B', ■ ■ ■ A^B^A^B'^: 



r^'^^'KOxx') = J duniU) { U®U®---®U) Qxx' { U^®U^®---® U\ ), 

2n 



2n 



(9.20) 
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where X — AiBi . . . A^B^ and Qxx' is a generic operator on l-Lxx' = ^as"- 
This channel is a proper generahzation of the map T^^^ whose properties can be 
estabhshed along the lines of (Zuber [2008]). With this choice we can express the 
n-th moment of the purity (3.52) as 



A^;^(l,0) = ^„ 



Tr 



02 |®n 



(2n) 
B 



(9.21) 



,B, <8 \^)a'.b'), I a is the identity on the 2n copies 

' 3 3^ 

of A, i.e. A — AiA^ ■ ■ ■ A^A'^, and is the swap operator which exchanges 



where |*®2^®" 

of A, i.e. A = 
subsystems B1B2 



Bn with B[B2 ■ ■ ■ B'^ pairwisely, i.e. 



(2n) 



Bj\B'.- 



9.1.2 Mixed initial states 

Consider now the case x <1. Following the parameterization introduced in Sec. 7.2 
we split the average over the set GxijH-x) of the density matrices of global purity 
X, as an average over the unitaries acting on l-Lx followed by an average over the 
space of the eigenvalues, see Eq. (7.16). Specifically this is accomplished by writing 

p{U) = Ux Ax , (9.22) 

with X — AB, Ux a generic unitary transformation on Ux, and Ax a diagonal 
density matrix which can be chosen as 

= Yl ^^i\^)A{^\®\i)B{i\ 

l<i<Nl<e<M 

= diag(AM(i-i)+j)i<i<Ar,i<€<jw- (9.23) 

Here, by keeping the same notation of Chap. 8, we have set Xu — XM{i-i)+e, where 
{\'i)A}i<i<N and B}i<e<M are orthonormal sets of Ha and Hb, respectively. 
Introducing then the ancillary systems a and b isomorphic to A and B respectively, 
we define the following purification of p, 

|*)ASa6= Yl ® ® |^)a <H) K)6 ■ (9.24) 

l<i<Ar l<e<M 

Let us fix the spectrum A^ of the global density matrix p of the system. The set 
of vectors in Hxx — Hx ® 'Hx with the same Schmidt coefficients is given by 

l*c/)x. = Uxx mxx ' Uxx ^Ux^U^e U{L^), (9.25) 
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with Ux,Ux e ^^{L) and yields the set of density matrices (9.22) with the same 
spectrum Ax- Here, we have set X — AB and x — ah. By partial tracing over 
subsystem B we obtain the set of reduced density matrices 

Pa{U) = TlBipABiU)) = Tl BabiUAB\'^)ABab (mis) 

ABab {'^\UU<l)Bab , (9.26) 

l<5<ArM2 

where {\q)Bab}i<q<NM'^ is an orthonormal basis of T-Lb ® 'Hab, and Uab G U{L). 
Notice that Pa{U) does not depend on U^- The local A-purity of Pa(U) becomes 



7rA(p) = IV(p/(t/)) = T^{Bab{q\UAB\'if)ABab{^\UiB\(l)Bab 

l<q,q'<NM^ 

XBab{q'\UAB\^) ABab q')Bab^ 



J2 ABab{'^\UUq)Bab{q'\UAB\^) ABab 
l<q,q'<NM^ 

X ABab ( * I C/is I g') Bab I C/aB I ABab ■ (9.27) 

In analogy to what seen for the pure case, it can be casted as an expectation value 
on l^f)®^ by doubling the Hilbert space 

7^a{p) = [xx{^\ ® X'x'i^l) (u\b ® ^A'B') {SsablB'a'b' IaA') 

X (Uab ® Ua'B') (I^)xx <E) , (9.28) 

where 

SBab\B'a'b' = \q)Bab{q \ ® \q)Bab{q\ (9.29) 

l<q,q'<NM^ 

is the swap operator between Bab and B'a'b'. Here, we have introduced the an- 
cillary bipartite systems X' = A'B', and x' = a'h' (here A',a',B',b' are the aux- 
iliary copies of A. a. B. and b respectively), thus |^)x.t = \'^)ABab and — 
\'^)A'B'a'b'- Therefore, by integrating over Uab we get, 

MtiAx) = Tr[r(2)(|*)x.(*| ® |*)XV(*|) [SBab\B'a'b'<^lAA')], (9.30) 

where T*-^-* is the twirling transformation on Hxx', Eq- (9.8). As mentioned in the 
previous subsection, from the properties of the twirling maps T^'^\ Aii{Ax) can 
be written as 

Mt{Ax) = [xx{^\^X'x'{^\)r^'^[SBab\B'a'b'^lAA'){\'^)xx^\'^)x'x'). (9.31) 
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See Eq. (9.9) for comparison. We will thus use the last identity. According to 
Eq. (9.11) we have to compute Trxx'(©xx') and Trxx'{Sx\x'Qxx') where Qxx' 
is 

SsablB'a'b' ® IaA' = Sb\B' <8) S^^jy <8) Sa\a' <H) IaA'- (9.32) 

We get 

T^'^ABA'B'iSsablB'a'b' ® ^AA') = T^T^ABA'B'{Sb\B' ® ® Sa\a' ® ^AA') 

= TybB'{Sb\B') TtaA'{IaA') Sb\b> (g) Sa\a' = N'^M Sb\b> (8) Sa\a' , (9.33) 

T^ABA'B'[SAB\A'B'{SBab\B'a'b' <H) ^Ayl')] = TtABA'B' (-^BB' <8) 'S'aIA' <H) 'S'bl;,/ (g) Sa\a') 
= TTBB'ilBB') TrAA'{SA\A') Sblb' Sa\a' = M^AT ^^|^, ® . (9.34) 

Thus, from Eq. (9.11) we obtain 

T*-^-* (^SBab\B'a'b' ® -^ylA' j = 



A^2M2 - 


1^^ 


Ar(M2 . 


-1) 


' mM^ 


- 1 


M{N^ - 




Ar2M2 - 






-1) 


' Ar2M2 


- 1 



<E) Sa\A' <8) '5'b|b' <8) 'S'ala' 

-S^ix' ® 'S'a;|a;', (9.35) 
with L = NM. Replacing this expression into Eq. (9.31) and using the identities 

l<i<N l<e<M 

(x.(*|<»xv(*|)(^x|x'0^.|.')(l*)x.0|*)xv) = 1, (9-36) 
we finally get 

A4^rA ^ - ^^^'-^V fAM I ^(^'-1) rq37l 
- iV2M2 - 1 + iV2M2-l' (^-^^^ 

which depends on the spectrum Ax only through the purity of p, i.e. Tr(A^) = 
Tr(p2). This immediately tells us that averaging upon Ax, while keeping fixed x, 
in (7.9), will give 

. .4/ ^^ M(N^ - 1) N(M^ - 1) 
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Some special ceises:- It is worth noticing that for a balanced bipartition N 
M = VI > 1 Eq. (9.38) yields 



iAt \/L{l+x) 1 + X 



Mt(..0) = ^^^^-^, (9.39) 

as found by using the Gaussian approximation in Chap. 8, see Eqs. (8.70) and 
(8.90). On the other hand for x = 1 (i.e. pure global states), Eq. (9.39) coincides 
with that obtained in (Giraud [2007a]; Facchi, Marzolino, Parisi, Pascazio and 
Scardicchio [2008]). 

Finally, consider the case in which p is maximally mixed, i.e. is the density 
matrix Iab/L. In this case x — 1/L and Eq. (9.38) correctly gives 

Mt{l/L,0)^j^, (9.40) 

in agreement with the general result (7.25). In analogy to what seen for pure states, 
the above analysis can in principle be extended to the case of higher moments of 

A1;t(Ax) = Tr[r(^")(|^)x.(^r(^"))(4'ab®^rO] ' (9-41) 

where = <^]=i{\'ip)AjBjajbj ® \'^) Ars'^arb'.) , I A is the identity on the 2n copies 

of A, i.e. A — AiA\ ■ ■ ■ A^A^, and S^^^^ is the swap operator which exchanges 
BiQibi with -B-a-6- for 1 < i < n, S^^^^ = <S>j=iSBjajbj\B'.arbr- However, forn > 1 Eq. 
(9.11) is no more valid and the application of the twirling transformation is far 
more complicated than for n = 1, since it cannot be reduced to a combination of 
projection operators on the symmetric and antisymmetric subspaces. In the next 
section we will introduce an alternative approach for the computation of moments 
of the local purity of an arbitrary mixed state, that will lead to a formula for the 
n-th moment of 7r^('0) and to the exact expression for its first two cumulants. 



9.2 Moments of the local purity at f3 = 

The technique we will introduce in this section in order to exactly compute the 
moments of the purity is based on the explicit solution due to Zuber of some basic 
integrals over the unitary group (Zuber [2008]). Notice that, as mentioned above, 
also the symmetry properties of the twirling maps can be determined along the lines 
of these solutions. We will exactly compute the first two moments and determine 
the dependence on the spectrum of the higher order ones. 
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9.2.1 First moment 

Let us introduce a purification for p on the doubled Hilbert space Tixx' = "Hx ® 
T-Lx'- According to what seen in the previous section, after integrating over all 
states in Hx with a given spectrum Ax, we will show that the only dependence on 
the latter is in terms of the purity x, fixed by the delta function in the partition 
function (7.9). If we fix the spectrum of the global density matrix p of the system 
as Ax = diag(AM(a~i)+/3) (see Eq. (9.23)), the purification of Ax in the space 
Hxx', with X = AB and X' = A'B' is 

|*)xx'= Yl E V^V \(^P)ab®Ma'B': (9-42) 

l<a<iV 1</3<M 

where = ® l/?)^, {|a)^}i<a<iv and {|/3)s}i</3<m A'}i<a<N and 

{|/^)b'}i</3<m) are the reference bases in T-La and T-Lb {'Ha' and ?^s'), respectively. 
We have seen that those vectors in Hxx' with the same Schmidt coefficients are 
given by \'^u)xx' = Uxx' |^)xxm Uxx' = Ux ® Ux' e U{L'^), with Ux,Ux' G 
hl{L) and yields the set of density matrices with the same spectrum Ax, namely 
p = Ux-^xUx- By partial tracing over subsystem B one obtains the set of reduced 
density matrices pAiU) = TrBTrx'(|^c/)xx' {"^u])- 

l<a<N 1<(3<M 

= E E ^c.0B{j\UAB\aP)^^{aP\U],^\j)^. (9.43) 

l<a<N l<l3,j<M 

Thus we can write the local purity as 

7rA([/xAxt/l) = IV(p/) = E E E V/3iAa,/3, 

l<oi,a2<A^ l</3i,/32<Af l<ji,j2<N 

XAB(a2/32|f^AB li2)B (ill ^AB 

X AB{aiPi\U[s O2I Uab ab > (9-44) 

which, by inserting the completeness relation for subsystem A, becomes 

nA{UxAxU]i) = ^ E E E ''^aiftAaafe 

l<ai,a2<N l</3i,/32<M l<ii,i2<N l<ji J2<M 

XAB{<y2P2\UAB KiJ2)ab (^l^ll ^AB {(^if^l) AB 

X AB{aiPi\U]^B N2Ji)ab (^2j2| Uab \0i2P2) ab ■ (9-45) 
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We can now compute the first moment of the purity (7.23) at P — 0. More generally, 
by recalling that 

Za{x,0) = 1, (9.46) 

dfJ,xfi{p) = dfix{p), (9.47) 

dl^Ap) = C,diiH{Ux)da{Ax)s(^x-Tr{A.j,)'^, (9.48) 

we get 

M^{x,Q) = j di,,{p)7rA''{p) 

^ C,j da{Ax)S{x - Tr(A^)) M^{Ax) = {M^i^x))., 

(9.49) 

with 

M^{Ax) = j di^niUx) 7r^"(C/xAxt/l). (9.50) 

Prom Eq. (9.45) we have that the average over the unitary group of the first moment 
particularizes to 

Ml (Ax) = XI XI -^ai/3iAa2/32 

l<ai,a2<A?" l</3i,/32<M 

X ^ / '^^^^^'^)^^i-?i'"i/^i^^2j2,a2/32^l2/32,iU2^cii/3i,i^^^^ 



X 

l<ii,i2<N l<ji,j2<M ' 

(9.51) 



where 



Uij,a(i = AB{tj\UAB\aP)^s V a, i G {1, . . . , TV}, V /3, j e {1, . . . , M} 

(9.52) 

are the matrix elements of Ux G U{L). This integral can be exphcitly done by 
applying the solution for the basic integral over the unitary group U{L) given by 
Zuber (Zuber [2008]): 

/n 
d^lH{U)Ui,j^ . . Mi^j^Uli^ . . .Ulj^ = X C[a]Yl5{iaJr{a))S{ja,kra(a)), 

(9.53) 

with 

^['^^ = ^ ^!^^;(7^ ' ^ ^ ^ 

|y|=n ^ 
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where C[a] is the sum over the Young diagrams Y where x^^\W]) is the character 
of the symmetric group Sn associated to Y, depending on the conjugacy class of 
the permutation a, Sk{I) is the dimension of the representation of Y in terms of 
the hnear group GL(L), and S{a,b) is the Kroneker delta function Sa^j,. Applying 
this solution to (9.51) we get 

Ml (Ax) = XI XI '^ai/SiAazft 

l<ai,a2<N l<0i,02<M 

X X ^M/l('^)<^('^l^l''^TO-(2)^rCT(2))(^(a2^2,aro-(l)^rCT(l)) 

T,<7eS2 

= E E E E C[a]Mr)S{[ras],c) 

l<ai,a2<N l<0i,02<M T,a&S2 ceC{S'2) 



(9.55) 



where /i (p) depends on the permutation p & S2 



flip) = E ^(^i'Mi))^(^2, V2)) E ^(ii'ip(2))^(i2,ip(i)), (9.56) 

i<n,i2<Ar i<jij2<M 

s is the transposition (swapping) of pairs of nearby indices (s = [2] G 5*2) 

is{i) = «2 and ^^(2) = ii (9.57) 

and C{S2) is the set of the conjugacy classes of the symmetric group 5'2, €(82) — 
{[1^], [2]}. From (9.55) it can be easily inferred that the only possible contributions 
of the spectrum are related to the conjugacy classes of the symmetric group 5*2 : 

[ras] = [1^] ^ f E E ^"/^) = 1' 

\1<Q<7V 1</3<M / 

[ras] ^[2] ^ f E E ^lp]-T^i^x)- (9-58) 

\l<a<Ar 1</3<M / 

By summing and by using the explicit expressions of the coefficients (9.54) (Zuber 
[2008]) 

^t^'l = (L-1)(L + 1) - = - (L-l)i(L+l) - f^-^^) 



we get 
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The first moment of the purity of subsystem A, is the average (9.60) over the 
spectrum of the system. By plugging (9.60) into (9.49), we finally get 

in total agreement with Eq. (9.38), derived by using the properties of the twirling 
transformations. 



9.2.2 k-th moment 

The technique shown in the previous section can be easily generalized in order to 
compute from (9.50) higher moments at ^ = 0. More precisely, we can write 



Mt(Kx) = I diJ.„(Ux) TA^iUx'^xUl, 



E Z n 

l<ai,...,a2fc<iV l<Pi,...Ak<M l<i<2fe 



X 



l<ii,...,i2k<N l<h,...j2k<M l<i<2k 

X TT r/^ 77^ 

J. J. CX2ml32 

l<m<fe 

(9.62) 

Thus, Eq. (9.53) for n ^ 2k gives 

MtiAx) = E E E E C[a]Mr)6{[ras],c) 

l<ai,...,a2k<N l<Pi,...,hk<M T,aeS2k ceC{S2k) 



X 



l<i<2k 



n 

(9.63) 

where fk{p) depends on the permutation p G S2k 

fk{p)= E n ^(^^'VW) E n ^ihm-lJp{2m))S{j2m:jp{2m-l)) 

l<ii,:.,i2k<N l<e<2k l<ji,---,j2k<M l<m<k 

(9.64) 

and, analogously to Eq. (9.57), s is the swapping of pairs of nearby indices, 1 < 
i < or 1 < i < M: 

is(2e-i) = i2e and is(2e) = i2e-i V£ e {1, . . . , k}. (9.65) 
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Observe that when A; = 1 we retrieve Mf{A.x) (see Eq. (9.55)). The different 
contributions from the spectrum can be classified in terms of the conjugacy classes 
of the symmetric group, as shown in Eq. (9.58). However, for k > 1, they do not 
depend only on the identity and the purity Tr(A^). They show a more complex 
dependence on the spectrum, thru its higher order invariants Tr(A^), with k > 2. 
Thus the integral on the spectrum (9.49) is in general non trivial. 



9.2.3 Second moment 



Let us now compute the second moment of the purity. In this section we will 
consider the case of arbitrary bipartite states, with purity x e [1/L, 1], generalizing 

some results found for the case of pure states x = 1, (Lubkin [1978]; Lloyd and 
Pagels [1988]; Scott and Caves [2003]; Giraud [2007a]; Facchi, Marzolino, Parisi, 
Pascazio and Scardicchio [2008]). The second moment of the purity can be directly 
computed by setting k — 2 in Eq. (9.62). The expression for the coefficients C[p] 
in (9.53), when p e €{84) is (Zuber [2008]): 




- 8L2 + 6 



(L - 3)(L - 2)(L - 1)L2(L + 1)(L + 2)(L + 3) ' 



C [2, 1^] 



1 




(L-3)(L-l)L(L + l)(L + 3)' 
L2 + 6 



C [3, 1] 



(L - 3)(L - 2)(L - 1)L2(L + 1)(L + 2)(L + 3) ' 
2L2-3 



C[4] 



{L - 3)(L - 2)(L - 1)L2(L + 1)(L + 2)(L + 3) ' 

5 



(9.66) 



(L - 3)(L - 2)(L - 1)L(L + 1)(L + 2)(L + 3) ' 
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The symmetric group 5*4 consists of five conjugacy classes, giving the following 
contributions to the integral (9.63) in terms of the spectrum of p: 



[ras] = [1'] 
[ras] = [2, 1^] 

[21 



irasi 



[ras] = [3, 1] 
[ras] = [4] 



Y^K.) =(Tr(A^))^ 



Tr(A^), 



(9.67) 



with r, 0" G 5*4 and ,s G 5*2 defined in (9.65), where, in this case, we set k = 2. 
We can now compute the expression for (Ax): 



M^{Ax) = CN., 



M 



{M' - 1){N^M\M^ - 1) - 2N\6M^ - 7) + 22) 

+Tr(A^) {2NM{N'^ - 1){M^ - 1){N^M'^ - 14)) 

+(Tr(A^))2 (A^^ - l)(M^Ar^ + M^A^^ _ u]y2j^2 ^ g^2 ^ gg^ 

+Tr(A^) 40(Ar^ - 1){M^ - 1) 



+Tr(A^) (-10ArM)(Ar2 - 1){M^ - 1) 



(9.68) 



where 



cn,m — 



N^M^{N^M^ - 7)2 - 36 



(9.69) 
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In particular ii M — N — \fL wc get 



(L^ - 2L^ - IIL^ + 26L2 + 8L - 22) 

+Tr(A^) (2L^ - 4L^ - 26L-^ + SGL^ - 28L) 
+ (Tr(A^))2 {1} - ISL^ + 20L2 + 24L - 30) 
+Tr(A^) 40(L - 1)=^ 



+IY(A^^) (-10L)(L-1)^ 



(9.70) 



with 



L2(L2- 7)2 -36' 



(9.71) 



Prom the exact expression for J\A^{x,^) and A1^(x,0), Eqs. (9.61) and (9.68) 
respectively, we can now compute the second cumulant of the purity at /3 = 0: 



/C^(x,0) = M^{x,^)- M^{x,^)' 



(iV2M2 - l)2(iV4M4 - 13iV2M2 + 36) 
4iVM(iV2 - 1)(M2 - l)(Ar2M2 + 11) 



— X 



(iV2M2 - l)2(iV4M4 - 13iV2M2 + 36) 
, 2(iV2 - 1)(M2 - l)(Ar4M4 - 4A^2^2 ^ 
(7V2M2 - l)2(Ar4M4 - 13A^2^2 ^ 36) 
40(A^2 _ i)(^2 _ 1) 



+ (Tr(Al)), 
-(Tr(A^)), 



7V2M2(Ar2M2 - 7)2 - 36 
10A^M(A^2 _ i)(^2 _ 1) 

7V2M2(A^27^4-2_7)2_36- 



(9.72) 
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For the case of balanced bipartitions, N — M — vT, we get 

2(L-1)2(L2 + 11) 



/C^(a;,0) 



(L2 - 1)2(L4 - 13L2 + 36) 

4L(L- 1)2(^2 + 11) 



—X 



(L2 - 1)2(L4 - 13L2 + 36) 
2 2(L- 1)2(L4-4L2 + 15) 
(L2 - 1)2(L4 - 13L2 + 36) 



L2(L2 - 7)2 - 36 
. 10L(L-1)2 



Notice that for the pure case, corresponding to = 1 for all k, we retrieve 

the results of Eqs. (3.59) and (3.57), as found in (Scott and Caves [2003]; Giraud 
[2007a]). Furthermore in the large limit we get 

^^(a^,0) ^, L^M^y/Z (9.74) 

as found in (Facchi, Marzolino, Parisi, Pascazio and Scardicchio [2008]), see Eq. 
(3.55). 



9.3 High temperature expansion of the first mo- 
ment of the purity 

In Chap. 7 we have seen that for x 7^ 1 due to the nonfactorization of the integral 
over the unitary group, the computation of the partition function is more involved 
than for pure states, with the only exception given by completely mixed states 
p — Ix/L. However, for sufficiently small /?, that is in the high temperature 
regime, we can study the statistical distribution of the local purity by expanding 
its n-th moment around /? = 0, up to order /J^, as shown in Eq. (7.26): 

Mt{x, /3) ^ Mt{x, 0) -/3 [Mt+,{x, 0) -M^{x, 0) Mt{x, 0)] , /3 ^ 0. (9.75) 
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In particular, by plugging in the above expression for n = 1 Eqs. (9.60) and (9.72) 
for the first moment and the second cumulant in ^ = we get 



Mt{x,P) ~ Mf{x,0) + P[{Mf{x,0)f-M^{x,Q)] 



—X 



—X 



1 " ' 

2{N^ - 1)(M2 - 1){N^M'^ + 11) 
(7V2M2 - lyimM"^ - IZN'^M^ + 36) 

2(iV2 - 1)(M2 - l)(-2iVM)(iV2M2 + 11) 
(Ar2M2 - 1)2(A^47\^4 _ 137V2M2 + 36) 

2 2(iV2 - 1)(M2 - 1)(A^4^4 _ 47V2M2 + 15) 



-(Tr(Ai)) 



■(Tr(A^)). 



(7V2M2 - l)2(Ar47i^4 _ i3Ar2M2 + 36) 
40(A^2 - 1)(M2 - 1) 



7V2M2(Ar2M2 - 7)2 - 36 

{-imM){N'^ - 1)(M2 - 1) 
N^M^iN^M"^ - 7)2 - 36 



(9.76) 



which for the case of balance bipartitions, N = M = \/L, particularizes to 



Mtix, 0)+/3 [{Mtix, 0))' - M^ix, 0)] 
y/L{l + x) 



1 + L 



2(12 + 11) 



+x 



(L + l)2(L4-13L2 + 36) 
4L(L2 + 11) 



-X 



(L + l)2(L4-13L2 + 36) 

2 2(L4 - 4L2 + 15) 
(L + l)2(L4-13L2 + 36) 
40(L- 1)2 



-(1V(A^)). 



L2(L2 -7)2 -36 
10L(L- 1)2 



L2(L2 - 7)2 - 36 



(9.77) 
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Eq. (9.77) for X = 1 can be compared with the results we have found in Chap. 3, 
where (3 was scaled as (3' = (3L^^'^ = (3N^. With this choice our expression yields 

M^{1, L'/^) ~ >lf^(l,0) + /3'L=^/2 [(7Wf(l,0))2-A^^(l,0)] 




(9.78) 



in exact agreement with the behavior of the average purity studied in Chap. 3. 
Indeed, by expanding the inverse temperature in Eq. (3.38) around 5 = 2 (corre- 
sponding to /3 = 0), we get 

P = --^ + l{5-2f + ((5 - 2f) (9.79) 

that is 

(5 = 2-4/3 + 24/3^ + O (/J^) . (9.80) 

Therefore, from Eq. (3.48) for the internal energy u = N {ttab) = \/LM.f{l, jSL'^^"^) 
we get 

u = 2(1 - /3) + O . (9.81) 

In Fig. 9.1 we show the behavior of the first moment of the local purity vr^ as a 
function of /3, for different values of the global purity ttx = a; G [1/L, 1] and for 
balanced bipartions N = M = a/Z. 




Figure 9.1: First moment of the purity as a function of /3. The red line refers to 
the set of totally mixed states, such that Ai^{l/L,(3) = 1/N V (3. The black one 
refers to pure states for /3 > 0, see Fig. 3.8. Finally the blue line corresponds to 
the first moment of the local purity at /3 = 0, for arbitrary mixed states, while the 
light blue region refers to the hight temperature region. 

Summarizing we have determined the exact analytical expression for Aif^{x, (3) for 
the set e^CHx) (7.2) of: 
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Exact moments and twirling maps 



• maximally mixed states, x = 1/L (see Eq. (7.25)) 

Mf{l/L,P)^l/N, V/3 (9.82) 

• arbitrarily mixed states x & [1/L, 1] at /3 = (see Eq. (9.39)) 

Mt{x,0)^y^^^, (9.83) 

• pure states x = 

Mt{x,P) = 
Mtix,p) = 

for < /? < /?+ 

Finally we have also 
seeEq. (9.77). 



1 and for every /3's (see Eqs. (3.48), (3.38) and (3.42)) 

(9.85) 

= 2 and /3+ < /3, respectively. 

found the high temperature expansion of the first moment. 



Conclusions and outlook 



In this thesis, we have discussed the behavior of bipartite entanglement of a large 
quantum system. We have seen that given a bipartite system, the purity of one part 
can be considered a measure of quantum correlations between the two subsystems 
only for the case of pure states (that is for isolated quantum systems), while for the 
more general case of mixed states, this quantity is just a lower bound, and should 
be substituted with its convex roof. In particular, we have developed a canonical 
approach for the study of the distribution of the Schmidt coefficients for a fixed 
value of the average entanglement. We have thus introduced a partition function 
depending on a fictitious temperature, which localizes the measure on the set of 
states with higher and lower entanglement with respect to typical (random) states, 
with respect to the Haar measure. The role of the energy in the partition function 
is played by our entanglement measure/indicator, the purity of one subsystem. 

In the first part of the thesis, we have obtained a complete characterization of 
the distribution of the purity and of the eigenvalues for the case of pure states. 
The global picture is interesting as several locally stable solutions exchange sta- 
bilities. On the stable branch (solutions of minimal/maximal free energy for pos- 
itive/negative temperatures) we have unveiled the presence of three main ther- 
modynamic phases for the system, namely the maximally entangled, the typical 
and the separable phases, separated by critical points. In particular we have en- 
lightened a second order phase transition, associated to a Z2 symmetry breaking 
and related to the vanishing of some Schmidt coefficients, followed by a first order 
phase transition, associated to the evaporation of the largest eigenvalue from the 
sea of the others. In the different phases the distribution of the Schmidt coefficients 
has very dissimilar profiles. While for large /9 (small purity) the eigenvalues, all 
0{1/N)^ follow the Wigner semicircle law, they become distributed according to 
the Wishart law for smaller /3 and larger purity, across the second order transition. 
For even smaller (and eventually negative) values of /3, when purity becomes finite, 
across the first order phase transition, one eigenvalue evaporates, leaving the sea 
of the other eigenvalues 0(1/ N) and becoming 0{1). This is the signature of 
separability, this eigenvalue being associated with the emergence of factorization 
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in the wave function (given the bipartition) . This interpretation is suggestive and 
hints at a profound modification of the distribution of the eigenvalues as and 
therefore the purity, are changed. Summarizing, for the stable solution the purity 
goes asymptotically from 1/N to 1. This analysis has been achieved by properly 
scaling /3, in order to yield the correct thermodynamic limit in the partition func- 
tion. By keeping the same scaling of positive temperatures also for negative /3's, 
we have determined a metastable solution for the system. It shows two phase tran- 
sitions that are both of second order. The first one, called gravity branch, takes 
place when the Wishart distribution, characterizing typical states, evolves into an 
asymmetric arcsine distribution. Then, with a second order phase transition, the 
distribution of the Schmidt coefficients returns symmetric, leading for /3 — )■ — oo to 
maximally entangled states. At the gravity branch we can interpret the expansion 
of the free energy around the critical temperature as the partition function of ran- 
dom 2D surfaces, yielding a theory of pure gravity (Morris [1991]; Di Francesco, 
Ginsparg and Zinn- Justin [1995]). 

In the second part of the thesis we have considered the more involved case of 
mixed states and through the same statistical approach, we have determined the 
exact expression of the first two cumulants at /3 = and the high temperature 
expansion, for the first moment of the local purity. We have thus generalized the 
results of (Giraud [2007a]; Scott and Caves [2003]), valid for the case of global pure 
states. Furthermore, we have also introduced an interesting connection between 
our problem and the theory of quantum channels, more precisely we have exploited 
the symmetry properties of the twirling transformations (VoUebrech and Werner 
[2001]) in order to compute the exact expression for the first moment of the local 
purity. 

Our characterization of bipartite entanglement of Haar-distributed states, where 
the least set of assumptions is made on their generation, can be used in an exper- 
iment, mutatis mutandis, as a check of the lack of correlations. If one observes 
that the moments of the purity deviate from the expected values, one could argue 
for non-randomness (or additional available information) of the states. In turn, 
our fictitious inverse temperature (3 acquires physical meaning, in that it measures 
deviations from typicality. This analysis also answers to another, more practical, 
question: since during the study of some features of entanglement one often re- 
lies on numerical simulations, it is important to know to what extent entangled 
quantum states may be considered typical (Zyczkowski, Horodecki, Sanpera and 
Lewenstein [1998]). 

The analysis performed in this thesis for the study of bipartite entanglement of 
large quantum systems, offers several possible perspectives for future investigation. 
From the point of view of quantum information, it would be of great interest to 
determine the evolution of the entanglement phase transitions we have studied for 
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the case of pure states, for the more general case of mixed states, when a new 
degree of freedom, given by the purity of the global state, is introduced. In this 
context, another outstanding result would be to understand whether these phase 
transitions survive even in the multipartite entanglement scenario (Facchi, Florio, 
Marzolino, Parisi and Pascazio [2010]), if one views the distribution of purity (over 
all balanced bipartitions) as a characterization of the global entanglement of the 
many-body wave function of quantum systems (Facchi, Florio, Parisi and Pascazio 
[2008]). 

Finally, we observe that the phase transitions investigated in this thesis, through a 
classical statistical mechanics approach arc not quantum phase transitions (Sachdev 
[1999]). Nevertheless, entanglement is known to be a good indicator of quantum 
phase transitions (Osterloh, Amico, Falci and Fazio [2002]), not only for what con- 
cerns bipartite entanglement (Vidal, Latorre, Rico and Kitaev [2003]; Verstraete, 
Popp and Cirac [2004]) but also for the case of multipartite entaglement (Costan- 
tini, Facchi, Florio and S. Pascazio [2007]; De Pasquale, Costantini, Facchi, Florio, 
Pascazio and Yuasa [2008]). Thus it would be interesting to investigate the link, if 
any, between these different phase transitions, and therefore apply the techniques 
shown in this work in order to explore new appealing features of many-body sys- 
tems, such as spin chains and systems which show a critical behavior. 

We believe that, since the approach followed in this thesis is rooted on the 
intersection of many aspects of theoretical physics (from classical statistical me- 
chanics, to quantum information, to random matrix theory) it can be regarded 
as a useful methodological framework for further and even more general topics of 
investigation. 



Appendix A 

Steepest descent method 



The steepest descent method is an approximation scheme to evaluate a class of 
contours integrals in a complex domain of the form 

In^ [ dz e^^(^) (A.l) 
Jc 

in the large N limit, being the contour C contained within the domain of ana- 
lyticity D of A. The steepest descent method is an application of the Cauchy's 
integral theorem according to which the contour of integration can be deformed in 
D without changing 1^. If we assume that the global maximum of A{z) is not on 
the boundary of D, it corresponds either to a singularity of the function or to a 
regular point where the derivative of A vanishes, 

A'{z) = 0. (A.2) 

We will focus on the second case. A point Zc where condition (A.2) is satisfied 
is, in general, a saddle point for the curves of constant ReA{z). In order to avoid 

cancellations due to the oscillatory behavior of the phase of the integrand we need 
to choose as contour of integration a path of steepest descent, that is a directed 
curve whose tangent at each point has a direction in which the rate of descent is 
maximal {direction of steepest descent). Indeed it can be proved that such curves 
are defined by 

lmA{z)^ImA{zc). (A.3) 

Thus along this path the phase of the integrand remains constant and the leading 
contribution to the integral comes from a neighbourhood of Zc- If we make the 
change of variables 

t = A{z,)-A{z) (A.4) 



153 



154 



Steepest descent method 



the path of steepest descent C can be mapped onto the positive real axis: 



p+oo 

/ dt A{t)e^\ A{t) 
Jo 



dz 



dA{z) 



(A.5) 

z=A-^A(zc)-t) 



We have assumed that the upper hmit of integration can be extended up to +00 
due to the fast exponential decay away from A{zc)- By applying the Watson's 
lemma for i — >■ 0+, in the limit N ^ 00, we find 



" 2N\A"(z. 



e^^(^=)+^^i-^ (A.6) 



where 9i and 62 arc the directions of steepest descent for the case A"{zc) 7^ 0, that 
is {z — Zc) = \z — Zc\e^^* for z & C and 1 < i < 2. For a more general form of the 
integral (A.l) 



I dzg{z) e^^(^), (A.7) 
Jc 

provided lng{z) is an analytic function in a region of the complex plane that 
includes C, we have to multiply expression (A.6) by the constant factor g{zc)- 
We remark that in most cases the leading term of this general expression can be 
retrieved by expanding A around the maximum and approximate expression (A.l) 
with a Gaussian integral. 

Another generalization of the above results is given by the n dimensional inte- 
gral on R" 

In^ J d^'z e^^(^i'^2,...^„)_ (A.8) 
In the limit — )■ cxd it can be shown that 

[detA(^)]-^/%^^(-), (A.9) 
being z*c = {^ci, ■ ■ ■ , -^cn} ^iid A^^^ the Hessian matrix 

= Vi,ie{l,...,n}. (A.10) 



Appendix B 



How to generate uniform random 
variables on the simplex 



In this appendix we discuss a useful trick in order to generate random numbers 
with the probabihty density function (8.52). 

Consider L i.i.d. random variables {yj)i<j<L with an exponential density function 

p{yj)=eM-yj), Vje{l,...,L} (B.l) 

Let us define Xj — 5^i<fc<£, yfc, for 1 < j < L. In what follows, in order to 
simplify the notation, if not specified the indices will run between 1 and L. The 
joint density function of the normalized random variables Xj is 

Pl{xu...:Xl) = / d'^ye-^^y^llsfxj-^^] 

J^^ J V z^kykJ 

= JJ^ dt e-* Jd^yS{t-J2 Vk) n ^ (^^ - f ) 

= ^ dt e-* J d^z t^6 (l - -Sfc) j n ^ " -^i) 

= y di t^-^e-^ j d^zs(i-j2 ^k) n ^ - 

= (L-l)!5(l-5^x,), (B.2) 
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where we have used the definition of the Gamma function 

poo 

r(L) ^ dt t^-^e-* (B.3) 
Jo 

and since L is an integer number, we get r(L) — [L — 1)!. Therefore, the random 
variables 

r/ = ^, Vje{l,...,L} (B.4) 
are distributed with joint probabihty density (8.52). 
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